B. Zwiebach 
February 9, 2016 


Chapter 1: Key Features of Quantum Mechanics 


Quantum mechanics is now almost one-hundred years old, but we are still discovering some of its 
surprising features and it remains the subject of much investigation and speculation. The framework 
of quantum mechanics is a rich and elegant extension of the framework of classical physics. It is also 
counterintuitive and almost paradoxical. 


Quantum physics has replaced classical physics as the correct fundamental description of our phys- 
ical universe. It is used routinely to describe most phenomena that occur at short distances. Quantum 
physics is the result of applying the framework of quantum mechanics to different physical phenomena. 
We thus have Quantum Electrodynamics, when quantum mechanics is applied to electromagnetism, 
Quantum Optics, when it is applied to light and optical devices, or Quantum Gravity, when it is 
applied to gravitation. Quantum mechanics indeed provides a remarkably coherent and elegant frame- 
work. The era of quantum physics begins in 1925, with the discoveries of Schrodinger and Heisenberg. 
The seeds for these discoveries were planted by Planck, Einstein, Bohr, de Broglie, and others. It is 
a tribute to human imagination that we have been able to discover the counterintuitive and abstract 
set of rules that define quantum mechanics. Here we aim to explain and provide some perspective on 
the main features of this framework. 

We will begin by discussing the property of linearity, which quantum mechanics shares with elec- 
tromagnetic theory. This property tells us what kind of theory quantum mechanics is and why, it 
could be argued, it is simpler than classical mechanics. We then turn to photons, the particles of light. 
We use photons and polarizers to explain why quantum physics is not deterministic and, in contrast 
with classical physics, the results of some experiments cannot be predicted. Quantum mechanics is a 
framework in which we can only predict the probabilities for the various outcomes of any given exper- 
iment. Our next subject is quantum superpositions, in which a quantum object somehow manages to 
exist simultaneously in two mutually incompatible states. A quantum light-bulb, for example, could 
be in a state in which it is both on and off at the same time! 


1 Linearity of the equations of motion 


In physics a theory is usually described by a set of equations for some quantities called the dynamical 
variables of the theory. After writing a theory, the most important task is finding solutions of the 
equations. A solution of the equations describes a possible reality, according to the theory. Because 
an expanding universe is a solution of Albert Einstein’s gravitational equations, for example, it follows 
that an expanding universe is possible, according to this theory. A single theory may have many 
solutions, each describing a possible reality. 

There are linear theories and nonlinear theories. Nonlinear theories are more complex than linear 
theories. In a linear theory a remarkable fact takes place: if you have two solutions you obtain a 
third solution of the theory simply by adding the two solutions. An example of a beautiful linear 
theory is Maxwell’s theory of electromagnetism, a theory that governs the behavior of electric and 
magnetic fields. A field, as you probably know, is a quantity whose values may depend on position 
and on time. A simple solution of this theory describes an electromagnetic wave propagating in a 
given direction. Another simple solution could describe an electromagnetic wave propagating in a 


different direction. Because the theory is linear, having the two waves propagating simultaneously, 
each in its own direction and without affecting each other, is a new and consistent solution. The sum 
is a solution in the sense that the electric field in the new solution is the sum of the electric field 
in the first solution plus the electric field in the second solution. The same goes for the magnetic 
field: the magnetic field in the new solution is the sum of the magnetic field in the first solution plus 
the magnetic field in the second solution. In fact you can add any number of solutions to still find 
a solution. Even if this sounds esoteric, you are totally familiar with it. The air around you is full 
of electromagnetic waves, each one propagating oblivious to the other ones. There are the waves of 
thousands of cell phones, the waves carrying hundreds of wireless internet messages, the waves from 
a plethora of radio-stations, TV stations, and many, many more. Today, a single transatlantic cable 
can carry simultaneously millions of telephone calls, together with huge amounts video and internet 
data. All of that courtesy of linearity. 

More concretely, we say that Maxwell’s equations are linear equations. A solution of Maxwell’s 
equation is described by an electric field E a magnetic field B, a charge density p and a current density 
J, all collectively denoted as (E, B,p, J). This collection of fields and sources satisfy Maxwell’s 
equations. Linearity implies that if (E, B,p, J) is a solution so is (@E, aB, ap, aJ), where all fields 
and sources have been multiplied by the constant a. Given two solutions 


(E1, Bi, 91,51), and (Eo, Bo, p2, Ja), (1.1) 
linearity also implies that we can obtain a new solution by adding them 
(E, + Ey, Bi + Bo, pi +p2, Jit+Jz). (1.2) 
The new solution may be called the superposition of the two original solutions. 
It is not hard to explain what is, in general, a linear equation or a linear set of equations. Consider 


the equation 
ba S00; (1.3) 


where, schematically, u denotes the unknown. The unknown may be a number, or a function of time, 
a function of space, a function of time and space, essentially anything unknown! In fact, u could 
represent a collection of unknowns, in which case we would replace u above by wi, u2,.... The symbol 
L denotes a linear operator, an object that satisfies the following two properties 


L(uy + u2) = Ilu,+ Lu, L(a u) = aLu, (1.4) 
where a is a number. Note that these conditions imply that 
L(au, + Bug) = alu, + BLue2, (1.5) 


showing that if u; is a solution ( Lu; = 0) and uz is a solution (Lug = 0) then au; + Buz is also a 
solution. We call au; + Suz the general superposition of the solutions u; and ug. An example may 
help. Consider the equation 


eae eS (1.6) 


where 7 is a constant with units of time. This is, in fact, a linear differential equation, and takes the 
form Lu = 0 if we define 


Lu = —+-u (1.7) 


Exercise 1. Verify that (1.7) satisfies the conditions for a linear operator. 


Einstein’s theory of general relativity is a nonlinear theory whose dynamical variable is a gravita- 
tional field, the field that describes, for example, how planets move around a star. Being a nonlinear 
theory, you simply cannot add the gravitational fields of different solutions to find a new solution. This 
makes Einstein’s theory rather complicated, by all accounts much more complicated than Maxwell the- 
ory. In fact, classical mechanics, as invented mostly by Isaac Newton, is also a nonlinear theory! In 
classical mechanics the dynamical variables are positions and velocities of particles, acted by forces. 
There is no general way to use two solutions to build a third. 

Indeed, consider the equation of motion for a particle on a line under the influence of a time- 
independent potential V(x), which is in general an arbitrary function of x. The dynamical variable 
in this problem is x(t), the position as a function of time. Letting V’ denote the derivative of V with 
respect to its argument, Newton’s second law takes the form 


Pat) _ vn 
aa = —V(a(t)). (1.8) 


The left-hand side is the mass times acceleration and the right hand side is the force experienced by 
the particle in the potential. It is probably worth to emphasize that the right hand side is the function 
V’(x) evaluated for x set equal to x(t): 


(1.9) 


While we could have used here an ordinary derivative, we wrote a partial derivative as is commonly 
done for the general case of time dependent potentials. The reason equation (1.8) is not a linear 
equation is that the function V’(z) is not linear. In general, for arbitrary functions u and v we expect 


V'(au) 4 aV'(u), and V'(utv) AV'(u)+V(v). (1.10) 


As a result given a solution x(t), the scaled solution ax(t) is not expected to be a solution. Given two 
solutions x;(t) and 2(t) then 2(t) + x2(t) is not guaranteed to be a solution either. 


Exercise. What is the most general potential V (a) for which the equation of motion for x(t) is linear? 


Quantum mechanics is a linear theory. The signature equation in this theory, the so-called 
Schrodinger equation is a linear equation for a quantity called the wavefunction and it determines 
its time evolution. The wavefunction is the dynamical variable in quantum mechanics but, curiously, 
its physical interpretation was not clear to Erwin Schrodinger when he wrote the equation in 1925. It 
was Max Born, who months later suggested that the wavefunction encodes probabilities. This was the 
correct physical interpretation, but it was thoroughly disliked by many, including Schrédinger, who 
remained unhappy about it for the rest of his life. The linearity of quantum mechanics implies a pro- 
found simplicity. In some sense quantum mechanics is simpler than classical mechanics. In quantum 
mechanics solutions can be added to form new solutions. 

The wavefunction W depends on time and may also depend on space. The Schrodinger equation 
(SE) is a partial differential equation that takes the form 


ih— = HV, (1.11) 


w 


where the Hamiltonian (or energy operator) H is a linear operator that can act on wavefunctions: 


H(av) =aHv, A(¥,+%.) = A(W,) +A), (1.12) 


with a a constant that in fact need not be real; it can be a complex number. Of course, A itself does 
not depend on the wavefunction! To check that the Schrédinger equation is linear we cast it in the 


form LW = 0 with L defined as 
woh he (1.13) 
ot 
It is now a simple matter to verify that L is a linear operator. Physically this means that if VW; and 
Wp» are solutions to the Schrddinger equation, then so is the superposition aW + GW.2, where a and 8 


are both complex numbers, i.e. (a, 6 € C) 


2 Complex Numbers are Essential 


Quantum mechanics is the first physics theory that truly makes use of complex numbers. The numbers 
most of us use for daily life (integers, fractions, decimals) are real numbers. The set of complex 
numbers is denoted by C and the set of real numbers is denoted by R. Complex numbers appear 
when we combine real numbers with the imaginary unit 7, defined to be equal to the square root of 
minus one: i = —1. Being the square root of minus one, it means that i squared must give minus 
one: i? = —1. Complex numbers are fundamental in mathematics. An equation like 2? = —4, for an 


unknown x cannot be solved if x has to be real. No real number squared gives you minus one. But if 


we allow for complex numbers, we have the solutions x = +27. Mathematicians have shown that all 
polynomial equations can be solved in terms of complex numbers. 
A complex number z, in all generality, is a number of the form 


z=a+ibeC, a,beER. (2.1) 


Here a and 0 are real numbers, and ib denotes the product of 7 with b. The number a is called the real 
part of z and 6 is called the imaginary part of z: 


Rez =a, Imz=b. (2.2) 
The complex conjugate z* of z is defined by 
SS 1b (2.3) 


You can quickly verify that a complex number z is real if z* = z and it is purely imaginary if z* = —z. 
For any complex number z = a+ib one can define the norm |z| of the complex number to be a positive, 


real number given by 
\z] = Va? +)?. (2.4) 


You can quickly check that 
|z/? = zz", (2.5) 


where z* = a — ib is called the complex conjugate of z = a+ 7b. Complex numbers are represented as 
vectors in a two dimensional “complex plane”. The real part of the complex number is the x component 
of the vector and the imaginary part of the complex number is the y component. If you consider the 
unit length vector in the complex plane making an angle 0 with the x axis has x component cos @ and 


y component sin 9. The vector is therefore the complex number cos @+7sin@. Euler’s identity relates 
this to the exponential of 70: 


e” — cosO+ ising. (2.6) 


A complex number of the form e’%, with x real is called a pure phase. 


While complex numbers are sometimes useful in classical mechanics or Maxwell theory, they are 
not strictly needed. None of the dynamical variables, which correspond to measurable quantities, 
is a complex number. In fact, complex numbers can’t be measured at all: all measurements in 
physics result in real numbers. In quantum mechanics, however, complex numbers are fundamental. 
The Schrodinger equation involves complex numbers. Even more, the wavefunction, the dynamical 
variable of quantum mechanics it itself a complex number: 


WeEC. (2.7) 


Since complex numbers cannot be measured the relation between the wavefunction and a measurable 
quantity must be somewhat indirect. Born’s idea to identify probabilities, which are always positive 
real numbers, with the square of the norm of the wavefuntion was very natural. If we write the 
wavefunction of our quantum system as W, the probabilities for possible events are computed from 
||?. The mathematical framework required to express the laws of quantum mechanics consists of 
complex vector spaces. In any vector space we have objects called vectors that can be added together. 
In a complex vector space a vector multiplied by a complex number is still a vector. As we will see in 
our study of quantum mechanics it is many times useful to think of the wavefunction V as a vector in 
some complex vector space. 


3 Loss of Determinism 


Maxwell’s crowning achievement was the realization that his equations of electromagnetism allowed 
for the existence of propagating waves. In particular, in 1865 he conjectured that light was an elec- 
tromagnetic wave, a propagating fluctuation of electric and magnetic fields. He was proven right in 
subsequent experiments. Towards the end of the nineteenth century physicists were convinced that 
light was a wave. The certainty, however, did not last too long. Experiments on blackbody radiation 
and on the photo-emission of electrons suggested that the behavior of light had to be more complicated 
than that of a simple wave. Max Planck and Albert Einstein were the most prominent contributors 
to the resolution of the puzzles raised by those experiments. 

In order to explain the features of the photoelectric effect, Einstein postulated (1905) that in a light 
beam the energy comes in quanta — the beam is composed of packets of energy. Einstein essentially 
implied that light was made up of particles, each carrying a fixed amount of energy. He himself found 
this idea disturbing, convinced like most other contemporaries that, as Maxwell had shown, light was 
a wave. He anticipated that a physical entity, like light, that could behave both as a particle and as 
a wave could bring about the demise of classical physics and would require a completely new physical 
theory. He was in fact right. Though he never quite liked quantum mechanics, his ideas about particles 
of light, later given the name photons, helped construct this theory. 

It took physicists until 1925 to accept that light could behave like a particle. The experiments of 
Arthur Compton (1923) eventually convinced most skeptics. Nowadays, particles of light, or photons, 
are routinely manipulated in laboratories around the world. Even if mysterious, we have grown 
accustomed to them. Each photon of visible light carries very little energy — a small laser pulse can 


contain many billions of photons. Our eye, however, is a very good photon detector: in total darkness, 
we are able to see light when as little as ten photons hit upon our retina. When we say that light 
behaves like a particle we mean a quantum mechanical particle: a packet of energy and momentum 
that is not composed of smaller packets. We do not mean a classical point particle or Newtonian 
corpuscle, which is a zero-size object with definite position and velocity. 

As it turns out, the energy of a photon depends only on the color of the light. As Einstein 
discovered the energy FE and frequency v for a photon are related by 


E=hv (3.1) 


The frequency of a photon determines the wavelength of the light through the relation vA = c, where 
c is the speed of light. All green photons, for example, have the same energy. To increase the energy 
in a light beam while keeping the same color, one simply needs more photons. 


As we now explain, the existence of photons implies that Quantum Mechanics is not deterministic. 
By this we mean that the result of an experiment cannot be determined, as it would in classical 
physics, by the conditions that are under the control of the experimenter. 

Consider a polarizer whose preferential direction is aligned along the x direction, as shown in 
Figure 1. Light that is linearly polarized along the x direction namely, light whose electric field points 
in this direction, goes through the polarizer. If the incident light polarization is orthogonal to the x 
direction the light will not go through at all. Thus light linearly polarized in the y direction will be 
totally absorbed by the polarizer. Now consider light polarized along a direction forming an angle a 
with the z-axis, as shown in Figure 2. What happens? 


Figure 1: A polarizer that transmits light linearly polarized along the x direction. 


| A 


Figure 2: Light linearly polarized along the direction at an angle a hitting the polarizer. 


Thinking of the light as a propagating wave, the incident electric field E, makes an angle a with 
the x-axis and therefore takes the form 


Ey = Egcosax+Eosinay. (3.2) 


This is an electric field of magnitude Eo. In here we are ignoring the time and space dependence 
of the wave; they are not relevant to our discussion. When this electric field hits the polarizer, the 
component along x goes through and the component along y is absorbed. Thus 


Beyond the polarizer: E = Eocosa x. (3.3) 


You probably recall that the energy in an electromagnetic wave is proportional to the square of the 
magnitude of the electric field. This means that the fraction of the beam’s energy that goes through 
the polarizer is (cosa). It is also well known that the light emerging from the polarizer has the same 
frequency as the incident light. 

So far so good. But now, let us try to understand this result by thinking about the photons that 
make up the incident light. The premise here is that all photons in the incident beam are identical. 
Moreover the photons do not interact with each other. We could even imagine sending the whole 
energy of the incident light beam one photon at a time. Since all the light that emerges from the 
polarizer has the same frequency as the incident light, and thus the same frequency, we must conclude 
that each individual photon either goes through or is absorbed. If a fraction of a photon went through 
it would be a photon of lower energy and thus lower frequency, which is something that does not 
happen. 

But now we have a problem. As we know from the wave analysis, roughly a fraction (cos a)? of the 
photons must go through, since that is the fraction of the energy that is transmitted. Consequently a 
fraction 1 — (cosa)? of the photons must be absorbed. But if all the photons are identical, why is it 
that what happens to one photon does not happen to all of them? 

The answer in quantum mechanics is that there is indeed a loss of determinism. No one can predict 
if a photon will go through or will get absorbed. The best anyone can do is to predict probabilities. 
In this case there would be a probability (cosa)? of going through and a probability 1 — (cosa)? of 
failing to go through. 

Two escape routes suggest themselves. Perhaps the polarizer is not really a homogeneous object 
and depending exactly on where the photon his it either gets absorbed or goes through. Experiments 
show this is not the case. A more intriguing possibility was suggested by Einstein and others. A 
possible way out, they claimed, was the existence of hidden variables. The photons, while apparently 
identical, would have other hidden properties, not currently understood, that would determine with 
certainty which photon goes through and which photon gets absorbed. Hidden variable theories would 
seem to be untestable, but surprisingly they can be tested. Through the work of John Bell and others, 
physicists have devised clever experiments that rule out most versions of hidden variable theories. No 
one has figured out how to restore determinism to quantum mechanics. It seems to be an impossible 
task. 

When we try to describe photons quantum mechanically we could use wavefunctions, or equivalently 
the language of states. A photon polarized along the x direction is not represented using an electric 
field, but rather we just give a name for its state: 


|photon; ay (3.4) 


We will learn the rules needed to manipulate such objects, but for the time being you could think of 
it like a vector in some space yet to be defined. Another state of a photon, or vector is 


|photon; Y) (3.5) 


representing a photon polarized along y. These states are the wavefunctions that represent the photon. 
We now claim that the photons in the beam that is polarized along the direction a are in a state 
|photon; a) that can be written as a superposition of the above two states: 


| photon; a) = cosa |photon; x) + sina |photon; y) ; (3.6) 


This equation should be compared with (3.2). While there are some similarities both are superpositions— 
one refers to electric fields and the other to “states” of a single photon. Any photon that emerges 
from the polarizer will necessarily be polarized in the x direction and therefore it will be in the state 


Beyond the polarizer: |photon; i (3:7) 


This can be compared with (3.3) which with the factor cosa carries information about the amplitude 
of the wave. Here, for a single photon, there is no room for such a factor. 


In the famous Fifth Solvay International Conference of 1927 the world’s most notable physicists 
gathered to discuss the newly formulated quantum theory. Seventeen out of the twenty nine attendees 
were or became Nobel Prize winners. Einstein, unhappy with the uncertainty in quantum mechanics 
stated the nowadays famous quote: “God does not play dice”, to which Niels Bohr is said to have 
answered: “Einstein, stop telling God what to do.” Bohr was willing to accept the loss of determinism, 
Einstein was not. 


4 Quantum Superpositions 


We have already discussed the concept of linearity; the idea that the sum of two solutions representing 
physical realities represents a new, allowed, physical reality. This superposition of solutions has a 
straightforward meaning in classical physics. In the case of electromagnetism, for example, if we have 
two solutions, each with its own electric and magnetic field, the “sum” solution is simply understood: 
its electric field is the sum of the electric fields of the two solutions and its magnetic field is the sum of 
the magnetic fields of the two solutions. In quantum mechanics, as we have explained, linearity holds. 
The interpretation of a superposition, however, is very surprising. 

One interesting example is provided by a Mach-Zehnder interferometer; an arrangement of beam 
splitters, mirrors, and detectors used by Ernst Mach and Ludwig Zehnder in the 1890’s to study 
interference between two beams of light. 

A beam splitter, as its name indicates, splits an incident beam into two beams, one that is reflected 
from the splitter and one that goes through the splitter. Our beam-splitters will be balanced: they 
split a given beam into two beams of equal intensity (Figure 3). The light that bounces off is called 
the reflected beam, the light that goes through is called the transmitted beam. The incident beam 
can hit the beam splitter from the top or from the bottom. 

The Mach-Zehnder configuration, shown in Figure 4, has a left beam splitter (BS1) and a right 
beam splitter (B52). In between we have the two mirrors, M1 on the top and M2 on the bottom. 
An incoming beam from the left is split by BS1 into two beams, each of which hits a mirror and is 
then sent into BS2. At BS2 the beams are recombined and sent into two outgoing beams that go into 
photon detectors DO and D1. 

It is relatively simple to arrange the beam-splitters so that the incoming beam, upon splitting at 
BS1 and recombination at BS2 emerges in the top beam which goes into DO. In this arrangement 
no light at all goes into D1. This requires a precise interference effect at BS2. Note that we have 
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Figure 3: An incident beam hitting a beam-splitter results in a reflected beam and a transmitted beam. Left: 
incident beam coming from the top. Right: incident beam coming from the bottom. 
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Figure 4: A Mach-Zehnder interferometer consists of two beam splitters BS1 and BS2, two mirrors M1 and 
M2, and two detectors DO and D1. An incident beam will be split into two beams by BS1. One beam goes 
through the upper branch, which contains M1, the other beam goes through the lower branch, which contains 
M2. The beams on the two branches recombine at BS2 and are then sent into the detectors. The configuration 
is prepared to produce an interference so that all incident photons end at the detector DO, with none at D1. 


two beams incident upon BS2; the top beam is called ‘a’ and the lower beam is called ‘b’. Two 
contributions go towards D0: the reflection of ‘a’ at BS2 and the transmission from ‘b’ at BS2. These 
two contributions interfere constructively to give a beam going into D0. Two contributions also go 
towards D1: the transmission from ‘a’ at BS2 and the reflection from ‘b’ at BS2. These two can indeed 
be arranged to interfere destructively to give no beam going into D1. 


It is instructive to think of the incoming beam as a sequence of photons that we send into the 
interferometer, one photon at a time. This shows that, at the level of photons, the interference is not 
interference of one photon with another photon. Each photon must interfere with itself to give the 
result. Indeed interference between two photons is not possible: destructive interference, for example, 
would require that two photons end up giving no photon, which is impossible by energy conservation. 

Therefore, each photon does the very strange thing of going through both branches of the inter- 
ferometer! Each photon is in a superposition of two states: a state in which the photon is in the top 
beam or upper branch, added to a state in which the photon is in the bottom beam or lower branch. 
Thus the state of the photon in the interferometer is a funny state in which the photon seems to be 
doing two incompatible things at the same time. 


Equation (3.6) is another example of a quantum superposition. The photon state has a component 
along an x-polarized photon and a component along a y-polarized photon. 

When we speak of a wavefunction, we also sometimes call it a state, because the wavefunction 
specifies the “state” of our quantum system. We also sometimes refer to states as vectors. A quantum 
state may not be a vector like the familiar vectors in three-dimensional space but it is a vector 
nonetheless because it makes sense to add states and to multiply states by numbers. Just like vectors 
can be added, linearity guarantees that adding wavefunctions or states is a sensible thing to do. Just 
like any vector can be written as a sum of other vectors in many different ways, we will do the same 
with our states. By writing our physical state as sums of other states we can learn about the properties 
of our state. 


Consider now two states | A) and |B : Assume, in addition, that when measuring some property 
Q in the state | A) the answer is always a, and when measuring the same property Q in the state |B ) 
the answer is always b. Suppose now that our physical state |W) is the superposition 


|v) = alA)+ |B), a: p.eC. (4.1) 


What happens now if we measure property Q in the system described by the state |W)? It may seem 
reasonable that one gets some intermediate value between a and b, but this is not what happens. A 
measurement of Q will yield either a or b. There is no certain answer, classical determinism is lost, but 
the answer is always one of these two values and not an intermediate one. The coefficients a and 6 in 
the above superposition affect the probabilities with which we may obtain the two possible values. In 
fact, the probabilities to obtain a or b 


Probability(a) ~ |a|?, | Probability(b) ~ |6|?. (4.2) 


Since the only two possibilities are to measure a or b, the actual probabilities must sum to one and 
therefore they are given by 


Jor? al? 
Ja|? + |)?” jal? + |B? 


If we obtain the value a, immediate repeated measurements would still give a, so the state after 
the measurement must be |A). The same happens for b, so we have 


Probability(a) = Probability(b) = (4.3) 


After measuring a the state becomes |W) = |A), (4.4) 
After measuring b the state becomes |V) = |B). , 


In quantum mechanics one makes the following assumption: Superposing a state with itself doesn’t 
chance the physics, nor does it change the state in a non-trivial way. Since superimposing a state with 
itself simply changes the overall number multiplying it, we have that VW and aW represent the same 
physics for any complex number a different from zero. Thus, letting = represent physical equivalence 


|A) = 2|A) = a] A) = -|A). (4.5) 


This assumption is necessary to verify that the polarization of a photon state has the expected number 
of degrees of freedom. The polarization of a plane wave, as one studies in electromagnetism, is described 
by two real numbers. For this consider an elliptically polarized wave, as shown in Figure 5. At any 
given point, the electric field vector traces an ellipse whose shape is encoded by the ratio a/b of 


10 


the semi-major axes (the first real parameter) and a tilt encoded by the angle @ (the second real 
parameter). Consider for this a general photon state formed by superposition of the two independent 
polarization states |photon; x) and |photon; y): 


a|photon; x) + 6|photon;y), a,BeEC. (4.6) 


At first sight it looks as if we have two complex parameters a and £, or equivalently, four real 
parameters. But since the overall factor does not matter we can multiply this state by 1/a to get the 
equivalent state that encodes all the physics 


|photon; x) + B |photon; y) , (4.7) 


showing that we really have one complex parameter, the ratio 3/a. This is equivalent to two real 
parameters, as expected. 


Figure 5: Parameters that define an elliptically polarized state. 


Let us do a further example of superposition using electrons. Electrons are particles with spin. 
Classically, we imagine them as tiny balls spinning around an axis that goes through the particle itself. 
Once an axis is fixed, the electron has two and only two options: its rotation may be clockwise or 
counterclockwise about the axis, but in both cases it spins at the same fixed rate. These opposite ways 
of spinning are called spin up and spin down along the axis (see Figure 6). The up and down refer to 
the direction of the angular momentum associated with the rotation, and it is indicated by an arrow. 
According to quantum mechanics, and as verified by multiple experiments, the same possibilities, up 
or down, arise whatever axis we use to measure the spin of the electron. 

Physicists usually set up coordinate systems in space by choosing three orthogonal directions, the 
directions of the x, y, and z axes. Let us choose to describe our spinning electrons using the z axis. 
One possible state of an electron is to be spin up along the z axis. Such a state is described as | +; z), 
with an arrow pointing up, and the label z indicating that the spin arrow points along the increasing z 
direction. Another possible state of an electron is spin down along the z axis. Such a state is described 
as | |; z), with an arrow pointing down, meaning this time that the spin points along the decreasing 
z direction. If these two are possible realities, so it would be the state |W) representing the sum 


[Y) = |t3z) + |L32). 


The state |W) is in a superposition of a spin up and a spin down state. What kind of physics does this 
sum |W) represent? It represents a state in which a measurement of the spin along the z axis would 
result in two possible outcomes with equal probabilities: an electron with spin up or an electron with 
spin down. Since we can only speak of probabilities, any experiment must involve repetition until 
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probabilities can be determined. Suppose we had a large ensemble of such electrons, all of them in 
the above state |W). As we measured their spin along z, one at a time, we would find about half of 
them spinning up along z and the other half spinning down along z. There is no way to predict which 
option will be realized as we measure each electron. It is not easy to imagine superposition, but one 
may try as follows. An electron in the above state is in a different kind of existence in which it is able 
to both be spinning up along z and spinning down along z simultaneously! It is in such a ghostly, 
eerie state, doing incompatible things simultaneously, until its spin is measured. Once measured, the 
electron must immediately choose one of the two options; we always find electrons either spinning up 
or spinning down. 


NZ AN 
| 
| 
| 
| 


Spin up Spin down 


== 


Figure 6: An electron with spin along the z axis. Left: the electron is said to have spin up along z. Right: 
the electron is said to have spin down along z. The up and down arrows represent the direction of the angular 
momentum associated with the spinning electron. 


A critic of quantum mechanics could suggest a simpler explanation for the above observations. He 
or she would claim that the following simpler ensemble results in identical experimental results. In the 
critic’s ensemble we have a large number of electrons with 50% of them in the state | +; z) and 50% of 
them in the state | |; z). He or she would then state, correctly, that such an ensemble would yield the 
same measurements of spins along z as the ensemble of those esoteric |W) states. The new ensemble 
could provide a simpler explanation of the result without having to invoke quantum superpositions. 

Quantum mechanics, however, allows for further experiments that can distinguish between the 
ensemble of our friendly critic and the ensemble of |W) states. While it would take us too far afield 
to explain this, if we measured the spin of the electrons in the x direction, instead of z direction, 
the results would be different in the two ensembles. In the ensemble of our critic we would find 50% 
of the electrons up along x and 50% of the electrons down along x. In our ensemble of |W) states, 
however, we would find a very simple result: all states pointing up along x. The critic’s ensemble 
is not equivalent to our quantum mechanical ensemble. The critic is thus shown wrong in his or her 
attempt to show that quantum mechanical superpositions are not required. 


5 Entanglement 


When we consider superposition of states of two particles we can get the remarkable phenomenon 
called quantum mechanical entanglement. Entangled states of two particles are those in which we 
can’t speak separately of the state of each particle. The particles are bound together in a common 
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state in which they are entangled with each other. 


Let us consider two non-interacting particles. Particle 1 could be in any of the states 


{|u1), |u2),--}, (5.1) 


while particle 2 could be in any of the states 


{|v1), \v2), -.-} (5.2) 


It may seem reasonable to conclude that the state of the full system, including particle 1 and particle 
2 would be specified by stating the state of particle 1 and the state of particle 2. If that would be the 
case the possible states would be written as 


jus) ® ley), 4FEN, (5.3) 


for some specific choice of 7 and 7 that specify the state of particle one and particle two, respectively. 
Here we have used the symbol ®, which means tensor product, to combine the two states into a single 
state for the whole system. We will study ® later, but for the time being we can think of it as a kind 
of product that distributes over addition and obeys simple rules, as follows 


(ay|u1) + a2|u2)) ® (B1]v1) + B2\v2)) = a1Bi\ur) @ |v1) + 0182] u1) ® |v2) 


(5.4) 
ar 0281] U2) ®& v1) 1 02/32 |u2) & |v2) : 


The numbers can be moved across the ® but the order of the states must be preserved. The state on 
the left-hand side -expanded out on the right-hand side— is still of the type where we combine a state 
of the first particle (a1 |u1) + a2|t2)) with a state of the second particle (81 |v1) + B2|v2)). Just like 
any one of the states listed in (5.3) this state is not entangled. 

Using the states in (5.3), however, we can construct more intriguing superpositions. Consider the 
following one 


|u1) ® |v1) + |u2) ® |ve). (5.5) 


A state of two particles is said to be entangled if it cannot be written in the factorized form (--- )@(---) 
which allows us to describe the state by simply stating the state of each particle. We can easily see 
that the state (5.5) cannot be factorized. If it could it would have to be with a product as indicated 
in (5.4). Clearly, involving states like |u3) or |v3) that do not appear in (5.5) would not help. To 
determine the constants a1, a2, 31,2 we compare the right hand side of (5.4) with our state and 
conclude that we need 


afi=1, amB2=0, aehi=0, arfg=1. (5.6) 


It is clear that there is no solution here. The second equation, for example, requires either a1 or 62 to 
be zero. Having a, = 0 contradicts the first equation, and having 32 = 0 contradicts the last equation. 
This confirms that the state (5.5) is indeed an entangled state. There is no way to describe the state 
by specifying a state for each of the particles. 


Let us illustrate the above discussion using electrons and their spin states. Consider a state of two 
electrons denoted as | t) ®| |). As the notation indicates, the first electron, described by the first 
arrow, is up along z while the second electron, described by the second arrow, is down along z (we 
omit the label z on the state for brevity). This is not an entangled state. Another possible state is 
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one where they are doing exactly the opposite: in | |) @ |) the first electron is down and the second 
is up. This second state is also not entangled. It now follows that by superposition we can consider 
the state 


It} @1L) + 148] 1). (5.7) 
This is a entangled state of the pair of electrons. 


Exercise. Show that the above state cannot be factorized and thus is indeed entangled. 


In the state (5.7) the first electron is up along z if the second electron is down along z (first term), 
or the first electron is down along z if the second electron is up along z (second term). There is a 
correlation between the spins of the two particles; they always point in opposite directions. Imagine 
that the two entangled electrons are very far away from each other: Alice has one electron of the pair 
on planet earth and Bob has the other electron on the moon. Nothing we know is connecting these 
particles but nevertheless the states of the electrons are linked. Measurements we do on the separate 
particles exhibit correlations. Suppose Alice measures the spin of the electron on earth. If she finds 
it up along z, it means that the first summand in the above superposition is realized, because in 
that summand the first particle is up. As discussed before, the state of the two particles immediately 
becomes that of the first summand. This means that the electron on the moon will instantaneously go 
into the spin down-along-z configuration, something that could be confirmed by Bob, who is sitting 
in the moon with that particle in his lab. This effect on Bob’s electron happens before a message, 
carried with the speed of light, could reach the moon telling him that a measurement has been done 
by Alice on the earth particle and the result was spin up. Of course, experiments must be done with 
an ensemble that contains many pairs of particles, each pair in the same entangled state above. Half 
of the times the electron on earth will be found up, with the electron on the moon down and the other 
half of the times the electron on earth will be found down, with the electron on the moon up. 

Our friendly critic could now say, correctly, that such correlations between the measurements of 
spins along z could have been produced by preparing a conventional ensemble in which 50% of the 
pairs are in the state | +) ®| {) and the other 50% of the pairs are in the state | |) ®@| t). Such 
objections were dealt with conclusively in 1964 by John Bell, who showed that if Alice and Bob are 
able to measure spin in three arbitrary directions, the correlations predicted by the quantum entangled 
state are different from the classical correlations of any conceivable conventional ensemble. Quantum 
correlations in entangled states are very subtle and it takes sophisticated experiments to show they are 
not reproducible as classical correlations. Indeed, experiments with entangled states have confirmed 
the existence of quantum correlations. The kind of instantaneous action at a distance associated with 
measurements on well-separated entangled particles does not lead to paradoxes nor, as it may seem, 
to contradictions with the ideas of special relativity. You cannot use quantum mechanical entangled 
states to send information faster than the speed of light. 


Sarah Geller transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this docu- 
ment. 
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1 Mach-Zehder Interferometer 


We have discussed before the Mach-Zehnder interferometer, which we show again in Figure 1. It 
contains two beam-splitters BS1 and BS2 and two mirrors. Inside the interferometer we have two 
beams, one going over the upper branch and one going over the lower branch. This extends beyond 
BS2: the upper branch continues to DO while the lower branch continues to D1. 


Figure 1: The Mach-Zehnder Interferometer 


Vertical cuts in the above figure intersect the two beams and we can ask what is the probability to 
find a photon in each of the two beams at that cut. For this we need two probability amplitudes, or 
two complex numbers, whose norm-squared would give probabilities. We can encode this information 
in a two component vector as 

a 
. ed: 
(5) an 


Here a is the probability amplitude to be in the upper beam and £ the probability amplitude to be in 
the lower beam. Therefore, |a|? would be the probability to find the photon in the upper beam and 
|G|? the probability to find the photon in the lower beam. Since the photon must be found in either 


one of the beams we must have 
lol? + |e)? =1. (1.2) 


Following this notation, we would have for the cases when the photon is definitely in one or the other 
beam: 


7) , photon on bottom beam : (7) : (1.3) 


photon on upper beam: G 


We can view the state (1.1) as a superposition of these two simpler states using the rules of vector 


(5) = (0) +(s) = (a) +2 (0) aa 


In the interferometer shown in Figure 1 we included in the lower branch a ‘phase shifter’, a piece of 


addition and multiplication: 


Figure 2: A phase shifter of phase factor e’®. The amplitude gets multiplied by the phase. 


material whose only effect is to multiply the probability amplitude by a fixed phase e” with 6 € R. 
As shown in Figure 2, the probability amplitude a to the left of the device becomes e’°a to the right 
of the device. Since the norm of a phase is one, the phase-shifter does not change the probability to 
find the photon. When the phase 6 is equal to 7 the effect of the phase shifter is to change the sign 
of the wavefunction since e’” = —1. 


Let us now consider the effect of beam splitters in detail. If the incoming photon hits a beam- 
splitter from the top, we consider this photon to belong to the upper branch and represent it by G 
If the incoming photon hits the beam-splitter from the bottom, we consider this photon to belong to 
4 . We show the two cases in Figure 3. The effect of the beam 


1 
splitter is to give an output wavefunction for each of the two cases: 


Left BS: a = @ , Right BS: (1) = (‘) (1.5) 


As you can see from the diagram, for the photon hitting from above, s may be thought as a reflection 


the lower branch, and represent it by ( 


amplitude and t as a transmission coefficient. Similarly, for the photon hitting from below, v may be 
thought as a reflection amplitude and u as a transmission coefficient. The four numbers s,t,u,v, by 
linearity, characterize completely the beam splitter. They can be used to predict the output given 
any incident photon, which may have amplitudes to hit both from above and from below. Indeed, an 


Qa 
would give 
;) 


()-«()+0Q) + «Ge -GE- 69). o9 


In summary, we see that the BS produces the following effect 


(5) > GG): an 


2 


incident photon state ( 


Figure 3: Left: A photon incident from the top; s and t are the reflected and transmitted amplitudes, re- 
spectively. Right: A photon incident from the bottom; v and u are the reflected and transmitted amplitudes, 
respectively. 


We can represent the action of the beam splitter as matrix multiplication on the incoming wavefunction, 


(; ‘) (1.8) 


We must now figure out the constraints on s,t, u,v. Because probabilities must add up to one, equation 
(1.5) implies that 


with the two-by-two matrix 


Isl + lel? 


ful” + fol? 


i (1.9) 
if (1.10) 


The kind of beam splitters we use are called balanced, which means that reflection and transmission 
probabilities are the same. So all four constants must have equal norm-squared: 


(1.11) 


Nie 


i. 4 
S U\) [Wf Wf 
(; y-@ #) (1.12) 
v2 V2 
This fails if acting on normalized wavefunctions (or column vectors) does not yield normalized wave- 
functions. So we try with a couple of wavefunctions 


4. i a Js a. 
= H@Q-&) &UE)-O ~ 
1 1}i\g om Ah fh ae 1) 
v2 V2 v2 V2 Va) \VB 


While the first example works out, the second does not, as |1|? + |1|? = 2 4 1. An easy fix is achieved 


by changing the sign of v: 
+ + 1/1 1 
S ") _ [| v2 | _ ( ) (1.14) 
= \). Lo) = ; : 
(; ' (2 V2 ya\iool 


a 


B 


ial) ))-de29) vs 


Indeed the resulting state is well normalized. The total probability is what we expect 


Let’s check that this matrix works in general. Thus acting on a state ( ) with |a|? +|8|? = 1 we find 


sla + 6] + 3/a— 8)” =3(lal” + |B)? + 06% + 08) + 3(\a|” + |B)’ 
=|o|? + |e)’ = 1, 


ee (1.16) 


The minus sign in the bottom right entry of (1.14) means that a photon incident from below, as it is 
reflected, will have its amplitude changed by a sign or equivalently a phase shift by a (check this!). 
This effect, of course, is realized in practice. A typical beam splitter consists of a glass plate with a 
reflective dielectric coating on one side. The refractive index of the coating is chosen to be intermediate 
between that of glass and that of air. A reflection causes a phase shift only when light encounters a 
material of higher refractive index. This is the case in the transition of air to coating, but not in the 
transition from glass to coating. Thus the beam splitter represented by (1.14) would have its coating 
on the bottom side. Transmitted waves have no phase shift. 
Another possibility for a beam splitter matrix is 


1 /-1 1 
et), aay 
which would be realized by a dielectric coating on the top side. You can quickly check that, like the 


previous matrix, its action also conserves probability. We will call the left beam-splitter BS1 and the 
right beam splitter BS2 and their respective matrices will be 


BS1: a4 iN BS2: = (; Aye (1.18) 


The two beam splitters are combined to form the interferometer shown in Figure 4. If we now assume 
a 


B 
is obtained by acting first with the BS1 matrix and then with the BS2 matrix: 


won (3) om St DEC YO AG IO-C)- 0» 


With the help of this result, for any input photon state we can write immediately the output photon 
state that goes into the detectors. 


an input photon wavefunction from the left, the output wavefunction that goes into the detectors 


, the output from the interferometer is , and therefore a 


1 
1 0 
photon will be detected at DO. This is shown in Figure 5. We can make a very simple table with the 
possible outcomes and their respective probabilities P: 


If the input photon beam is c 


Outcome P 
| photonat DO | 1 (1.20) 
| photon at D1 


mM 


Figure 4: The Mach-Zehnder interferometer with input and output wavefunctions indicated. 


a 
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Figure 5: Incident photon from below will go into DO. 


Now, block the lower path, as indicated in Figure 6. What happens then? It is best to track down 
things systematically. The input beam, acted by BS1 gives 


i f/-1. 1\70\.. + 71 
Al 1 , () 32 () | ie 
This is indicated in the figure, to the right of BS1. Then the lower branch is stopped, while the upper 
1 


branch continues. The upper branch reaches BS2, and here the input is c , because nothing is 


coming from the lower branch. We therefore get an output 


=(; - (#) 2 ()) (1.22) 


In this experiment there are three possible outcomes: the photon can be absorbed by the block, or 


Figure 6: The probability to detect the photon at D1 can be changed by blocking one of the paths. 


can go into any of the two detectors. As we see in the diagram, the probabilities are: 


Outcome P 

photon at block 5 
(1.23) 

photon at DO : 

I 

photonatD1 | 7 


It is noteworthy that before blocking the lower path we could not get a photon to D1. The probability 
to reach D1 is now 1/4 and was increased by blocking a path. 


2 Elitzur-Vaidman Bombs 


To see that allowing the photon to reach D1 by blocking a path is very strange, we consider an imag- 
inary situation proposed by physicists Avshalom Elitzur and Lev Vaidman, from Tel-Aviv University, 
in Israel. They imagined bombs with a special type of trigger: a photon detector. A narrow tube goes 
across each bomb and in the middle of the tube there is a photon detector. To detonate the bomb 
one sends a photon into the tube. The photon is then detected by the photon detector and the bomb 
explodes. If the photon detector is defective, however, the photon is not detected at all. It propagates 
freely through the tube and comes out of the bomb. The bomb does not explode. 

Here is the situation we want to address. Suppose we have a number of Elitzur-Vaidman (EV) 
bombs, but we know that some of them have become defective. How could we tell if a bomb is 
operational without detonating it? Assume, for the sake of the problem, that we are unable to 
examine the detector without destroying the bomb. 

We seem to be facing an impossible situation. If we send a photon into the detector tube and 
nothing happens we know the bomb is defective, but if the bomb is operational it would simply 
explode. It seems impossible to confirm that the photon detector in the bomb is working without 
testing it. Indeed, it is impossible in classical physics. It is not impossible in quantum mechanics, 
however. As we will see, we can perform what can be called an interaction-free measurement! 

We now place an EV bomb on the lower path of the interferometer, with the detector tube properly 
aligned. Suppose we send in a photon as pictured. If the bomb is defective it is as if there is no detector, 
the lower branch of the interferometer is free and all the photons that we send in will end up in DO, 


Figure 7: A Mach-Zehnder interferometer and an Elitzur-Vaidman bomb inserted on the lower branch, with 
the detector tube properly aligned. If the bomb is faulty all incident photons will end up at DO. If a photon ends 
up at D1 we know that the bomb is operational, even though the photon never went into the bomb detector! 


just as they did in Figure 5. 


Outcome P 
photonatD0O | 1 
no explosion 
photonatD1 | 0 
no explosion 

bomb explodes | 0 


(2.24) 


If the bomb is working, on the other hand, we have the situation we had in Figure 6, where we placed 
a block in the lower branch of the interferometer: 


Outcome P 
bomb explodes $ 
photon at DO i (2.25) 
no explosion 
photon at D1 i 
no explosion 


Assume the bomb is working. Then 50% of the times the photon will hit it and it will explode, 25% 
of the time the photon will end in DO and we can’t tell if it is defective or not. But 25% of the time 
the photon will end in D1, and since this was impossible for a faulty bomb, we have learned that 
the bomb is operational! We have learned that even though the photon never made it through the 
bomb; it ended on D1. If you think about this you will surely realize it is extremely surprising and 
counterintuitive. But it is true, and experiments (without using bombs!) have confirmed that this 
kind of interaction-free measurement is indeed possible. 


Sarah Geller transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this docu- 
ment. 
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1 Photoelectric Effect 


The photoelectric effect was first observed by Heinrich Hertz in 1887. When polished metal plates are 
irradiated, he observed, they may emit electrons, then called “photo-electrons”. The emitted electrons 
thus produce a photoelectric current. The key observations were: 


e There is a threshold frequency vg. Only for frequencies v > vg is there a photoelectric current. 
The frequency vo depends on the metal and the configuration of the atoms at the surface. It is 
also affected by inhomogeneities. 


e The magnitude of the photoelectric current is proportional to the intensity of the light source. 
e Energy of the photoelectrons is independent of the intensity of the light source. 


A natural explanation for the features in this effect didn’t come until 1905, when Einstein explained 
the above features by postulating that the energy in light is carried by discrete quanta (later called 
photons) with energy hy. Here h is Planck’s constant, the constant used by Planck to to produce a 
fit for the blackbody energy as a function of frequency. 


AES ts 


Figure 1: Electrons in a metal are bound. If the photon energy is greater than the work function W 
an electron may be ejected. 


A given material has a characteristic energy W, called the work function, which is the minimum 
energy required to eject an electron. This is not easily calculated because it is the result of an 


interaction of many electrons with the background of atoms. It is easily measured, however. When 
the surface of the material is irradiated, electrons in the material absorb the energy of the incoming 
photons. If the energy imparted on an electron by the absorption of a single photon is greater than 
the work function W, then the electron is ejected with kinetic energy E,- equal to the difference of 
the photon energy and the work function: 


E.- = mv = hv-W = E,-W. (1.1) 


This equation, written by Einstein explains the experimental features noted above, once we assume 
that the quanta act on individual electrons to eject them. The threshold frequency is defined by 


hvo = W, (1.2) 


as it leads to a photoelectron with zero energy. For v > vp the electrons will be ejected. Increasing the 
intensity of the light source increases the rate that photons arrive, which will increase the magnitude 
of the current, but will not change the energy of the photoelectrons because it does not change the 
energy of each incoming quanta. 

Equation (1.2) allowed Einstein to make a prediction: The kinetic energy of the photo-electrons 
increases linearly with the frequency of light. Einstein’s prediction was confirmed experimentally 
by Millikan (1915) who measured carefully the photoelectron energies and confirmed their linear 
dependence on the energy. Millikan’s careful work allowed him to determine the value of Planck’s 
constant h to better than 1% accuracy! Still, skepticism remained and physicists were not yet convinced 
about the particle nature of these light quanta. 


Example: Consider UV light with wavelength \ = 290nm incident on a metal with work function 
W =4.05eV What is the energy of the photo-electron and what is its speed? 


Solution: It is useful to solve these problems without having to look up constants. For this try 
recalling this useful relation 


he = 197.33 MeV.fm, h= (1.3) 


where MeV = 10°%eV and fm= 10~!°m. Let us use this to compute the photon energy. In this case, 
Cc 2m - 197.33 MeV.fm 27 - 197.33 


Ey =hv= aah = OS 10a - 500 eV & 4.28eV, (1.4) 
and thus 
E.- = Ey, — W = 0.23eV. (1.5) 
To compute the energy we set 
0.23eV = gmev? = 1 (mec?) (2)? (1.6) 
Recalling that mec? ~ 511,000eV one finds 
46 = (2)? a} 2 = 0,0000488., (1.7) 


With and c = 300,000 Km/s we finally get v ~ 284.4 Km/s. 


This is a good point to consider units, in particular the units of h. We can ask: Is there a physical 
quantity that has the units of h. The answer is yes, as we will see now. From the equation EF = hv, 
we have 


yp Ma, (1.8) 


n= ) _ MI?/T? | c 


V 


2 


where [-] gives the units of a quantity, and M,L,T are units of mass, length, and time, respectively. 
We have written the right-most expression as a product of units of length and momentum. Therefore 


[A] = [rx p] = [L]. (1.9) 


We see that A has units of angular momentum! Indeed for a spin one-half particle, the magnitude of 
the spin angular momentum is sh. 


With [h] = [r][p] we also see that one has a canonical way to associate a length to any particle 
of a given mass m. Indeed, using the speed of light, we can construct the momentum p = mc, and 
then the length @ is obtained from the ratio h/p. This actually is the Compton wavelength (c of a 
particle: 


pene (1.10) 

mec 
then has units of length; this is called the Compton wavelength of a particle of mass m. Note that this 
length is independent of the velocity of the particle. The de Broglie wavelength of the particle uses 
the true momentum of the particle, not mc! Thus, Compton and de Broglie wavelengths should not 


be confused! 


It is possible to get some physical intuition for the Compton wavelength Ac of a particle. We claim 
that Ac is the wavelength of a photon whose energy is equal to the rest energy of the particle. Indeed 
we would have 
me = hy = ne oe (1.11) 
mc 
confirming the claim. Suppose you are trying to localize a point particle of mass m. If you use light, 
the possible accuracy in the position of the particle is roughly the wavelength of the light. Once we 
use light with A < Ac the photons carry more energy than the rest energy of the particle. It is possible 
then that the energy of the photons go into creating more particles of mass m, making it difficult, if 
not impossible to localize the particle. The Compton wavelength is the length scale at which we need 
relativistic quantum field theory to take into account the possible processes of particle creation and 
annihilation. 
Let us calculate the Compton wavelength of the electron: 
h 27he 27 - 197.33 MeV.fm 


Wie = = = 2426fm = 2.426 pm. 1.12 
ole) ae 0.511 MeV ca pm (1.12) 


This length is about 20 times smaller than the Bohr radius (53 pm.) and about two-thousand times 
the size of a proton (1 fm.). The Compton wavelength of the electron appears in the formula for the 
change of photon wavelength in the process called Compton scattering. 


2 Compton Scattering 


Originally Einstein did not make clear that the light quantum meant a particle of light. In 1916, 
however, he posited that the quantum would carry momentum as well as energy, making the case for 
a particle much clearer. In relativity, the energy, momentum, and rest mass of a particle are related 
by 

BS oe Se (2.13) 


(Compare this with the classical equation E = p*/2m.) Of course, one can also express the energy 
and momentum of the particle in terms of the velocity: 


E = 5 p ——————— (2.14) 


You should use these expressions to confirm that (2.13) holds (|p| = p). A particle that moves with 
the speed of light, like the photon, must have zero rest mass, otherwise its energy and momentum 
would be infinite due to the vanishing denominators. With the rest mass set to zero, equation (2.13) 
gives the relation between the photon energy E, and the photon momentum py: 


Ey = Dy, (2.15) 
Then, using Av = c, we reach 
_E, _ bw _h (2.16) 
re iY 


We will see this relation again later when we discuss matter waves. 


Figure 2: Unpolarized light incident on an electron scatters into an angle 6. Classically, this is described 
by Thomson scattering. The light does not change frequency during this process. 


Compton carried out experiments (1923-1924) scattering X-rays off a carbon target. X-rays corre- 
spond to photon energies in the range from 100 eV to 100 KeV. The goal was scattering X-ray photons 
off free electrons, and with some qualification, the electrons in the atoms behave this way. 

The classical counterpart of the Compton experiment is the scattering of electromagnetic waves 
off free electrons, called Thompson scattering. Here an electromagnetic wave is incident on a electron. 
The electric field of the wave shakes the electron which oscillates with the frequency of the incoming 
field. The electron oscillation produces a radiated field, of the same frequency as that of the incoming 
radiation. In classical Thomson scattering the differential scattering cross section is given by 


do e2 \? 1 9 
40 (5) 5 (1+ cos” 6) , (217) 
where @ is the angle between the incident and scattered wave, with the radiated energy at the same 
frequency as the incoming light. This is shown in Figure 2. The cross-section has units of length- 
squared, or area, as it should. It represents the area that would extract from the incoming plane wave 
the amount of energy that is scattered by the electron. Indeed the quantity e?/(mc7) is called the 
classical electron radius and it is about 2.8 fm! not much bigger than a proton! 

If we treat the light as photons, the elementary process going on is a collision between two particles; 
an incoming photon and a roughly stationary electron. Two facts can be quickly demonstrated: 


e The photon cannot be absorbed by the electron. It is inconsistent with energy and momentum 
conservation (exercise) 


e The photon must lose some energy and thus the final photon wavelength A; must be larger 
than the initial photon wavelength ;. This is clear in the laboratory frame, where the initially 
stationary electron must recoil and thus acquire some kinetic energy. 
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Figure 3: The results of Compton’s scattering experiment. The incident photon wavelength is ;, and 
the scattered photon wavelength is Af ~ A; + €c, corresponding to @ = 90°. 


Indeed, Compton’s observations did not agree with the predictions of Thompson scattering: the 
X-rays changed frequency after scattering. A calculation using energy and momentum conservation 
shows that the change of wavelength is correlated with the angle between the scattered photon and 
the original photon: 


h 
At = Ai + — (1 — cos 8) = A4,+ fc (1 —cos8). (2.18) 


Note that appearance of the Compton wavelength of the electron, the particle the photon scatters off 
from. The maximum energy loss for the photon occurs at 6 = 7, where 


Ap(@ = 180°) = A3 + 2AC. (2.19) 
The maximum possible change in wavelength is 2A¢. For 6 = 5 the change of wavelength is exactly lc 
Ap(9 = 90°) = AV + Ac. (2.20) 

Compton’s experiment used molybdenum X-rays with energy and wavelength 
E,17.5keV, ; = 0.0709nm, (2.21) 


incident on a carbon target. Placing the detector at an angle 6 = 90° the plot of the intensity (or 
number of photons scattered) as a function of wavelength is shown in Figure 2. One finds a peak for 
for Af = 0.0731 nm, but also a second peak at the original wavelength A; = 0.0709 nm. 

The peak at Ay is the expected one: Af — Ay Y 2.2 pm, which is about the Compton wavelength 
of 2.4 pm. Given that the photons have energies of 17 KeV and the bound state energies of carbon 


are about 300 eV, the expected peak represents instances where the atom is ionized by the collision 
and it is a fine approximation to consider the ejected electrons. The peak at A; represents a process 
in which an electron receives some momentum from the photon but still remains bound. This is not 
very unlikely: the typical momentum of a bound electron is actually comparable to the momentum of 
the photon. In this case the photon scatters at 90° and the recoil momentum is carried by the whole 
atom. The relevant Compton wavelength is therefore that of the atom. Since the mass of the carbon 
atom is several thousands of times larger than the mass of the electron, the Compton wavelength of 
the atom is much smaller than the electron Compton wavelength and there should be no detectable 
change in the wavelength of the photon+ 


3 Matter Waves 


As we have seen, light behaves as both a particle and a wave. This kind of behavior is usually said to 
be a duality: the complete reality of the object is captured using both the wave and particle features 
of the object. The photon is a particle of energy E,, but has frequency v which is a wave attribute, 
with E = hv. It is a particle with momentum p, but it also has a wavelength A, a wave attribute, 
given by (2.16) 


A= (3.22) 


h 
Py 
In 1924, Louis de Broglie proposed that the wave/particle duality of the photon was universal, and 
thus valid for matter particles too. In this way he conjectured the wave nature of matter. Inspired by 
(3.22) de Broglie postulated that associated to a matter particle with momentum p there is a plane 
wave of wavelength given by 
h 
\=-. 
Pp 
This is a fully quantum property: if h > 0, then ’ — 0, and the particles have no wave properties. 


(3.23) 


And exciting consequence of this is that matter particles can diffract or interfere! In the famous 
Davisson-Germer experiment (1927) electrons are strike a metal surface and one finds that at certain 
angles there are peaks in the intensity of the scattered electrons. The peaks showed the effect of 
constructive interference from scattering off the lattice of atoms in the metal, demonstrating the wave 
nature of the electrons. One can also do two-slit interference with electrons, and the experiment can 
be done shooting one electron at a time. A recent experiment [arXiv:1310.8343] by Eibenberger et.al 
reports interference using molecules with 810 atoms and mass exceeding 10 000 amu (that’s 20 million 
times the mass of the electron!) 

The de Broglie wavelength can be calculated to estimate if quantum effects are important. Consider 
for this purpose a particle of mass m and momentum p incident upon an object of size x, as illustrated 
in Figure 3. Let A = h/p denote the de Broglie wavelength of the particle. The wave nature of the 
particle is not important if is much smaller than x. Thus, the “classical approximation,” in which 
wave effects are negligible, requires 


A 
Wave effects negligible: - ale (3.24) 


Using \ = h/p, this yields 
Wave effects negligible: zrp>h, (3.25) 
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Figure 4: A particle of momentum p incident on an obstacle of size x. 


A 


a relation in which both sides have units of angular momentum. 

Classical behavior is a subtle limit of quantum mechanics: a classical electromagnetic field requires 
a large number of photons. Any state with an exact, fixed number of photons, even if large, is not 
classical, however. Classical electromagnetic states are so-called coherent states, in which the number 
of photons fluctuates. 


Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this 
document. 
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1 de Broglie wavelength and Galilean transformations 


We have seen that to any free particle with momentum p, we can associate a plane wave, or a “matter 
wave”, with de Broglie wavelength \ = h/p, with p = |p|. The question is, waves of what? Well, this 
wave is eventually recognized as an example of what one calls the wavefunction. The wavefunction, 
as we will see is governed by the Schrédinger equation. As we have hinted, the wavefunction gives us 
information about probabilities, and we will develop this idea in detail. 

Does the wave have directional or polarization properties like electric and magnetic fields in an 
electromagnetic wave? Yes, there is an analog of this, although we will not delve into it now. The 
analog of polarization corresponds to spin! The effects of spin are negligible in many cases (small 
velocities, no magnetic fields, for example) and for this reason, we just use a scalar wave, a complex 
number 


U(x,t)EeCc (1.1) 


that depends on space and time. A couple of obvious questions come to mind. Is the wavefunction 
measurable? What kind of object is it? What does it describe? In trying to get intuition about this, 
let’s consider how different observers perceive the de Broglie wavelength of a particle, which should 
help us understand what kind of waves we are talking about. Recall that 
h h 20 
baa See (1.2) 
where & is the wavenumber. How would this wave behave under a change of frame? 

We therefore consider two frames S$ and S’ with the x and 2’ axes aligned and with S’ moving to 
the right along the +z direction of S with constant velocity v. At time equal zero, the origins of the 
two reference frames coincide. 

The time and spatial coordinates of the two frames are related by a Galilean transformation, which 
states that 

g=z2—vt, t=t. (1.3) 


Indeed time runs at the same speed in all Galilean frames and the relation between x and 2’ is manifest 
from the arrangement shown in Fig. 1. 

Now assume both observers focus on a particle of mass m moving with nonrelativistic speed. Call 
the speed and momentum in the S frame v and p = mv, respectively. It follows by differentiation with 
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Figure 1: The S’ frame moves at speed v along the x-direction of the S frame. A particle of mass m 
moves with speed v, and thus momentum p = mv, in the S$ frame. 


respect to t = ¢’ of the first equation in (1.3) that 


dx! dz 
_ ae 1.4 
dt! a or 


which means that the particle velocity v’ in the S’ frame is given by 
tifa, (1.5) 
Multiplying by the mass m we find the relation between the momenta in the two frames 
p = p-—mv. (1.6) 


The momentum p’ in the S$’ frame can be appreciably different from the momentum p in the S' frame. 
Thus the observers in S” and in S$ will obtain rather different de Broglie wavelengths \’ and X\! Indeed, 
h h 


ye — 
AS = 
ph p—mv 


TaN (1.7) 


This is very strange! As we review now, for ordinary waves that propagate in the rest frame of 
a medium (like sound waves or water waves) Galilean observers will find frequency changes but no 
change in wavelength. This is intuitively clear: to find the wavelength one need only take a picture 
of the wave at some given time, and both observers looking at the picture will agree on the value of 
the wavelength. On the other hand to measure frequency, each observers must wait some time to see 
a full period of the wave go through them. This will take different time for the different observers. 


Let us demonstrate these claims quantitatively. We begin with the statement that the phase 
ob = kx — wt of such a wave is a Galilean invariant. The wave itself may be cos@ or sing or some 
combination, but the fact is that the physical value of the wave at any point and time must be agreed 
by the two observers. The wave is an observable. Since all the features of the wave (peaks, zeroes, 
etc, etc) are controlled by the phase, the two observers must agree on the value of the phase. 

In the S frame the phase can be written as follows 

20 27x 2nV 


where V = @ is the wave velocity. Note that the wavelength is read from the coefficient of x and w is 


minus the coefficient of t The two observers should agree on the value of ¢. That is, we should have 


¢(2',t') = $(2,t) (1.9) 


where the coordinates and times are related by a Galilean transformation. Therefore 


20 


¢' (2', t’) = +z (2 - Ve) = Te + ot! Ve) = F Qnr(V —v),, 


: 1.1 
12 i t (1.10) 


Since the right-hand side is expressed in terms of primed variables, we can read X’ from the coefficient 
of x’ and w’ as minus the coefficient of t’: 


NX (1.11) 
uf = =v v) = at (1 =) =w (1 = =) (1.12) 


This confirms that, as we claimed, for a physical wave propagating in a medium, the wavelength is a 
Galilean invariant and the frequency transforms. 

So what does it mean that the wavelength of matter waves change under a Galilean transforma- 
tion? It means that the UY waves are not directly measurable! Their value does not correspond to a 
measurable quantity for which all Galilean observers must agree. Thus, the wavefunction need not be 
invariant under Galilean transformations: 


W(ar,t) 4 W(2',t’), (1.13) 


where (a,t) and (2’,t’) are related by Galilean transformations and thus represent the same point and 
time. You will figure out in Homework the correct relation between U(x,t) and W/(a’,t’). 


What is the frequency w of the de Broglie wave for a particle with momentum p? We had 
p=hk (1.14) 


which fixes the wavelength in terms of the momentum. The frequency w of the wave is determined by 
the relation 
E = hw, (1.15) 


which was also postulated by de Broglie and fixes w in terms of the energy F of the particle. Note 
that for our focus on non-relativistic particles the energy E is determined by the momentum through 


the relation 
Pr 
R= —., 1.16 
om (1.16) 


We can give three pieces of evidence that (1.15) is a reasonable relation. 


1. If we superpose matter waves to form a wave-packet that represents the particle, the packet 
will move with the so called group velocity vg, which in fact coincides with the velocity of the 
particle. The group velocity is found by differentiation of w with respect to k, as we will review 
soon: 


dw dE d (2) Dp 


= = vy. 1.1 
dk dp dp \2m m i“ kt) 


2. The relation is also suggested by special relativity. The energy and the momentum components 
of a particle form a four-vector: 
E 
=p (1.18) 


Similarly, for waves whose phases are relativistic invariant we have another four-vector 
w 
(= ; k) (1.19) 


Setting two four-vectors equal to each other is a consistent choice: it would be valid in all Lorentz 
frames. As you can see, both de Broglie relations follow from 


(=.») =o; (Ss k) (1.20) 


3. For photons, (1.15) is consistent with Einstein’s quanta of energy, because E = hy = hw. 


In summary we have 


pH=hk, H=he. (1.21) 


These are called the de Broglie relations, and they are valid for all particles. 


2 Phase and Group Velocities 


To understand group velocity we form wave packets and investigate how fast they move. For this we 
will simply assume that w(k) is some arbitrary function of k. Consider a superposition of plane waves 
ei(kt—w(k)t) given by 


(oa= / dk &(k)eithe-W()0), (2.22) 


We assume that the function ®(k) is peaked around some wavenumber k = ko, as shown in Fig. 2. 


Figure 2: The function ®(k) is assumed to peak around k = ko. 


In order to motivate the following discussion consider the case when ®(k) not only peaks around 
ko but it also is real (we will drop this assumption later). In this case the phase y of the integrand 
comes only from the exponential: 

p(k) = kx —wu(k)t. (2:23) 


We wish to understand what are the values of x and t for which the packet ~(z,t) takes large values. 

We use the stationary phase principle: since only for k ~ ko the integral over k has a chance to give a 

non-zero contribution, the phase factor must be stationary at k = ko. The idea is simple: if a function 

is multiplied by a rapidly varying phase, the integral washes out. Thus the phase must have zero 

derivative at kg. Applying this idea to our phase we find the derivative and set it equal to zero at ko: 
dp dw 


This means that (x,t) is appreciable when x and ¢ are related by 


dw 
= _ 2.25 
x dk ko d ( ) 
showing that the packet moves with group velocity 
dw 
Sli i 2.26 
"9 dk lig ee) 


Exercise. If ®(kg) is not real write ®(k) = |®(k)|e’*). Find the new version of (2.25) and show that 
the velocity of the wave is not changed. 


Let us now do a more detailed calculation that confirms the above analysis and gives some extra 
insight. Notice first that 


Hedy = i, dk ®t, (2.27) 
We expand w(k) in a Taylor expansion around k = ko 
dw 2 
w(k) = w(ko) + (k — ko) 7 ee O ((k — ko)*) . (2.28) 
ko 
Then we find, neglecting the O((k — ko)”) terms 
: : —~i(k— dw 
wot) = ja B(k) eth ew (ho)te ve) di lao (2.29) 


It is convenient to take out of the integral all the factors that do not depend on k: 
w(x t) = ae Caer tho tho! f ak &(k)ctte th ala" 
p) 


iw(ko)ttiky | t ik(x—4#| +4 20) 
= gee 0) +t 0 Fe lig fa @(k)e (= 41,5 ) 


Comparing with (2.27) we realize that the integral in the above expression can be written in terms of 
the wavefunction at zero time: 


UG ee A ee © 


x“- 


—iw(ko )t+iko % kot v( “ 


. t). (2.31) 


The phase factors in front of the expression are not important in tracking where the wave packet is. 
In particular we can take the norm of both sides of the equation to find 


dw 


Iv(@,t)| =|v(@- = 


t, 0) . (2.32) 


If ¢(a,0) peaks at some value xo it is clear from the above equation that |w(x,t)| peaks for 


dw 
t=% 7 «= 2+ —| Ft, (2.33) 
ie dk 


ko 


showing that the peak of the packet moves with velocity vg, = a evaluated at ko. 


3 Choosing the wavefunction for a free particle 


What is the mathematical form of the wave associated with a particle a particle with energy EF and 
momentum p? We know that w and k are determined from EF = homega and p = hk. Let’s suppose 
that we want our wave to be propagating in the +% direction. All the following are examples of waves 
that could be candidates for the particle wavefunction. 


1. sin (kx — wt) 
2. cos (kx — wt) 
3. el(ku—wt) — eikx e—iwt _ time dependence x e~i* 


4, et(ka—wt) — e—ike piwt _ time dependence x ett 


In the third and fourth options we have indicated that the time dependence could come with either 
sign. We will use superposition to decide which is the right one! We are looking for a wave-function 
which is non-zero for all values of x. 


Let’s take them one by one: 


1. Starting from (1), we build a superposition in which the particle has equal probability to be 
found moving in the +z and the —x directions. 


W(a,t) = sin (ka — wt) + sin (kx + wt) (3.1) 
Expanding the trigonometric functions this can be simplified to 
W(x,t) = 2sin(kx) cos(wt). (3:2) 


But this result is not sensible. The wave function vanishes identically for all x at some special 


times 
wt = (5, F, 3...) (3.3) 
A wavefunction that is zero cannot represent a particle. 
2. Constructing a wave function from (2) with a superposition of left and right going cos waves, 
W(x,t) = cos(kx — wt) + cos(kx + wt) = 2cos(kx) cos(wt) . (3.4) 
This choice is no good, it also vanishes identically when wt = (G, ar ss) 


3. Let’s try a similar superposition of exponentials from (3), with both having the same time 


dependence 
U(x, t) = eilka—wt) + el(-ke—wt) (3.5) 
ae (ae ae ete) e wt 
= 2coskre™. (3.7) 
This wavefunction meets our criteria! It is never zero for all values of 2 because e~*” is never 


Zero. 


4. A superposition of exponentials from (4) also meets our criteria 


U(x, t) Ler Aig i(ka—wt) ie i(—ka—wt) (3.8 
= (en a en) eivt (3 9 
= 2coskre™. (3.10) 


This is never zero for all values of x 


Since both options (3) and (4) seem to work we ask: Can we use both (3) and (4) to represent a 
particle moving to the right (in the +%# direction)? Let’s assume that we can. Then, since adding a 
state to itself should not change the state, we could represent the right moving particle by using the 
sum of (3) and (4) 

U(x, t) = ete) 4 e-Uke—wt) — 9 cos(kx — wt). (3.11) 


This, however, is the same as (2), which we already showed leads to difficulties. Therefore we must 
choose between (3) and (4). 


The choice is a matter of convention, and all physicists use the same convention. We take the free 
particle wavefunction to be 


Free particle wavefunction: (x,t) = e' : (3.12) 


representing a particle with 
p=hk, and B=hw. (3.13) 


In three dimensions the corresponding wavefunction would be 


Free particle wavefunction: W(x,t) = eer ery, (3.14) 


representing a particle with 
p=hk, and E=hw. (3.15) 
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1 Equations for a wavefunction 


We determined that the wavefunction or de Broglie wave for a particle with momentum p and 
energy E is given by 
VE =e (1.1) 


where w and k are determined from 
p=hk, E=hw, E=— (1.2) 


The wavefunction (1.1) represents a state of definite momentum. It is then of interest to find an 
operator that extracts that information from the wavefunction. The operator must be, roughly, 
a derivative with respect to x, as this would bring down a factor of k. In fact more precisely, 
we take 


#2 W(2,1) = (aK) U(e 0) 
= AkW(z,t) (1.3) 
= pW (a, t) 


where the p factor in the last right-hand side is just the momentum. We thus identify the 


fie) - 
operator 75> as the momentum operator p 


(1.4) 


and we have verified that acting on the wavefunction WV(z,t) for a particle of momentum p it 
gives p times the wavefunction: 
pv = pv. (1.5) 


The momentum operator it acts on wavefunctions, which are functions of space and time to 
give another function of x and t. Since p on WV gives a number (p, in fact) times V we say 


that V is an eigenstate of p. The matrix algebra analogy is useful: matrices are the operators 
and column vectors are the states. Matrices act by multiplication on column vectors. An 
eigenvector of a matrix is a special vector. The matrix acting on an eigenvector gives a number 
times the eigenvector. After the action of the matrix the direction of the vector is unchanged 
but its magnitude can be scaled. The same for eigenstates of operators: an operator acting on 
an eigenstate gives the eigenstate up to a multiplicative constant. We also say that W is a state 
of definite momentum. 


Let us now consider extracting the energy information from the free particle wavefunction. This 
time we must avail ourselves of the time derivative: 


ue. ; 

thy U(a, t) = iA(—iw)W(2,t) = hwWV(az,t) = EWV(z,t). (1.6) 
It would seem plausible to say that the time derivative ind is an energy operator but, for a 
free particle the energy is given in terms of the momentum, so we can construct the relevant 
energy operator by working on the above right-hand side 


EV = vy = py = —-—v, (1.7) 


where we used equation (1.5) to write pW as the momentum operator acting on WV. Since p is 
a constant we can move the p factor on the last right-hand side close to the wavefuntion and 
then replace it by the momentum operator: 


1 ho 1hoOhO 
EW = —-——pV = —-—-—-—— VW. 1.8 
am idx? 2m 1 Ox i Ox ee) 
This can be written as ; 
1 F 
Ev = —ppv = 2 vw, (1.9) 
2m 2m 
which suggests the following definition of the energy operator BE: 
. ~) K2 @? 
Se ce eee ca (1.10) 
2m 2m Ox? 


Indeed, for our free particle wavefunction, (1.9) shows that EW = EW. 
Our work also allows us to find an differential equation for which our de Broglie wavefunction 
is a solution. Consider (1.6) and replace the right hand side EW by E'W giving us 


(1.11) 


This is the free-particle Schrodinger equation. More schematically, using the energy operator, 
it can be written as 


(1.12) 


It is worth re-checking that our de Broglie wavefunction satisfies the Schrodinger equation 
(1.11). Indeed for UV = e*(**-“4) we find 


ih(—iw)W = ——(ik)?U (1.13) 


which is a solution since the WV factors cancel and all that is needed is the equality 


hehe 


2m 


hw = 


(1.14) 


which is recognized as the familiar relation E = fap 


Note that the Schrodinger equation admits more general solutions than the de Broglie wave- 
function for a particle of definite momentum and definite energy. Since the equation is linear, 
any superposition of plane wave solutions with different values of k is a solution. Take for 
example 

U(x, t) = el(kiz—wit) ua ei(k2a—wat) (1.15) 


This is a solution, and note that while each summand corresponds to a state of definite mo- 
mentum, the total solution is not a state of definite momentum. Indeed 


pVU(z,t) = Ak fo + fiky 2-9 | (1.16) 


and the right hand side cannot be written as a number times V. The full state is not a state of 
definite energy either. The general solution of the free Schrodinger equation is the most general 
superposition of plane waves: 


U(x,t) = i dk ®(k) et*e—¥ (6) (1.17) 
where ®(k) is an arbitrary function of k that controls the superposition and we have written 
w(k) to emphasize that w is a function of the momentum, as in (1.14). 


Exercise. Verify that WV in (1.17) solves the free Schrédinger equation. 


We now have the tools to time-evolve any initial wavefunction. Namely, given the initial 
wavefunction V(x,0) of any packet at time equal zero, we can obtain W(x,t). Indeed, using 
Fourier transformation one can write 


W(x,0) = [came (1.18) 


where ®(k) is the Fourier transform of U(z,0). But then, the time evolution simply consists 
in adding the exponential e~)! to the integral, so that the answer for the time evolution is 
indeed given by (1.17). 

As we have discussed before, the velocity of a wave packet described by (1.17) is given by the 
group velocity evaluated for the dominant value of k. We confirm that this is indeed reasonable 


Ow Ohw OE O/(Pr D 
fa == 1.1 
Ok Ohk Op Op ( ) ey) 


Ug = 


2m 


3 


which is the expected velocity for a free non-relativistic particle with momentum p and mass m. 


The Schrodinger equation has an explicit 7 on the left-hand side. This 7 shows that it is 
impossible to find a solution for real UV. If V were real the right-hand side of the equation would 
be real but the left-hand side would be imaginary. Thus, the Schrodinger equation forces us to 
work with complex wavefunctions. 

Note also that the Schrodinger equation does not take the form of a conventional wave 
equation. A conventional wave equation for a variable ¢ takes the form 

2 2 
OD we LOD as. fy: (1.20) 
Ogee Veo 
The general solutions of this linear equation are f(a + Vt). This would certainly allow for 
real solutions, which are not acceptable in quantum theory. The Schrodinger equation has no 
second-order time derivatives. It is first-order in time! 


2 Schrodinger Equation for particle in a potential 


Suppose now that our quantum particle is not free but rather is moving in some external 
potential V(x, t). In this case, the total energy of the particle is no longer simply kinetic, it is 
the sum of kinetic and potential energies: 


P 
E = —+4V(z,t), (2.1) 
2m 
This naturally suggests that the energy operator should take the form 


A p 
B= —+4V(qg,t). 2.2 
P+ V (0,1) (2.2) 
The first term, as we already know, involves second derivatives with respect to x. The second 
term acts multiplicatively: acting on any wavefunction W(z, t) it simply multiplies it by V(z, t). 
We now postulate that the Schrodinger equation for a particle in a potential takes the form 
(1.12) with EF’ replaced by the above energy operator: 


O h? 0? 
7 pe ey nae 2, 
ee!) ( On 23) 


(2.4) 


(2.5) 


Let us reconsider the way in which the potential V(z,t) is an operator. We can do this by 
introducing a position operator < that acting on functions of x gives another function of x as 
follows: 


Ef(x) = xf(x). (2.6) 
Note that it follows from this equation and successive applications of it that 
a f(a) Se" fla). (2.7) 


If the potential V(z,t) can be written as some series expansion in terms of x it then follows 
that 
V(4,t)W(2,t) = V(a,t)V(a,t). (2.8) 
The operators we are dealing with (momentum, position, Hamiltonian) are all declared to 
be linear operators. A linear operator A satisfies 


A(ad) =aAd, A(di+¢2) = Adit Ago, (2.9) 


where a is a constant. Two linear operators A and B that act on the same set of objects can 
always be added (A + Bo = Ad + Bo. They can also be multiplied, the product AB is a 
linear operator defined by ABd = A(Be), meaning that you act first with B, which is closest 
to @ and then act on the result with A. The order of multiplication matters and thus AB and 
BA may not be the same operators. To quantify this possible difference one introduces the 
commutator |A, B] of two operators, defined to be the linear operator 


(2.10) 


If the commutator vanishes, the two operators are said to commute. It is also clear that 
[A, A] = 0 for any operator A. 

We have operators ¢ and p that are clearly somewhat related. We would like to know 
their commutator [@, p|. For this we let |%, 6] act on some arbitrary function ¢(a) and then 


attempt simplification. Let’s do it. 
[2 , p\o(x) = (&p — p&)O(x) = Fp d(x) — pk O(z) 
= £(po(x)) — p(#o(z)) 


ho 
= (2D) _ iao(ay) 
we ho pa 
Fae ag eH) ati 
_ fh Ogle) fh Oe(a) _ fh 
i" Oc i Ox Go 
= ~24(2) = thg(2), 
so that, all in all, we have shown that for arbitrary $(x) one has 
[&, p]o(e) = ihe(e). (2.12) 


Since this equation holds for any ¢ it really represents the equality of two operators. Whenever 
we have Ad = Bo for arbitrary ¢ we simply say that A = B. The operators are the same 
because they give the same result acting on anything! We have therefore discovered the most 
fundamental commutation relation in quantum mechanics: 


[,p] = ih. (2.13) 


The right hand side is a number, but should be viewed as an operator (acting on any function it 
multiplies by the number). This commutation relation can be used to prove Heisenberg’s uncer- 
tainty principle, which states that the product of the position uncertainty and the momentum 
uncertainty must always exceed h/2. 

The idea that operators can fail to commute may remind you of matrix multiplication, which 
is also non-commutative. We thus have the following correspondences: 


operators ++ matrices 
wavefunctions <> vectors (2.14) 


eigenstates << eigenvectors 


One can in fact formulate Quantum Mechanics using matrices, so these correspondences are 
actually concrete and workable. 

As an example of useful matrices that do not commute, consider the Pauli matrices, three 
two-by-two matrices given by 


o1= ( 7 , 02= ¢ a 03 = (; a (2.15) 


Actually these matrices are exactly what is needed to consider spin one-half particles. The spin 
operator S has three components S; = Loy. Let us now see if o; and co» commute. 


ma) G a)-G 4) 
ie c a (? : _ & . (2.16) 


2% 0 fl. 0 
[o1, 72] = & >») = 20 6 ay = 2ios (2.17) 


In fact, one also has [02,03] = 2io1 and [o3, 01] = 2ioo. 

Matrix mechanics, was worked out in 1925 by Werner Heisenberg and clarified by Max 
Born and Pascual Jordan. Note that, if we were to write £ and p operators in matrix form, 
they would require infinite dimensional matrices. One can show that there are no finite size 
matrices that commute to give a number times the identity matrix, as is required from (2.13). 
This shouldn’t surprise us: on the real line there are an infinite number of linearly independent 
wavefunctions, and in view of the correspondences in (2.14) it would suggest an infinite number 
of basis vectors. The relevant matrices must therefore be infinite dimensional. 


We then see that 


Two basic properties of the Schrodinger equation 


1. The differential equation is first order in time. This means that for an initial condition it 
suffices to know the wavefunction completely at some initial time tp and the Schrodinger 
equation then determines the wave function for all times. This can be understood very 
explicitly. If we know WV(z, tg) for all x then the right-hand side of the Schrédinger equa- 
tion, which just involves x derivatives and multiplication, can be evaluated at any point x. 
This means that at any point xz we know the time-derivative of the wavefunction (left- 
hand side of the Schrédinger equation) and this allows us to calculate the wavefunction 
a little time later. 


2. Linearity and superposition. The Schrodinger equation is a linear equation for complex 
wavefunctions. Therefore, given two solutions V; and V2, we can form new solutions as 
linear combinations aW,; + 6WV»2 with complex coefficients a and (. 


We have written the Schrodinger equation for a particle on a one-dimensional potential. 
How about for the case of a particle in a three-dimensional potential? As we will see now, this 
is easily done once we realize that in three dimensions the position and momentum operators 
have several components! Recall that the de Broglie wavefunction 


U(x, t) = elk x—wt) = etl kga+kyytkzz —wt) (2.18) 


corresponds to a particle carrying momentum p = hk, with k = (k,, ky, kz). Just as we did in 
(1.3) we can try to extract the vector momentum by using a differential operator. The relevant 
operator is the gradient: 


0 0 O 
Ve= Ga (2.19) 
with which we try 
h De Gee daa 2 
7 VU, i= ; (tks; tkyyth,) U(e,t) = AEV(z, 1) = pvG@,d). (2.20) 


We therefore define the momentum operator p as follows: 


p = ay. (2.21) 


If we call the momentum components (p1, p2, 3) = (Px, Py, Pz) and the coordinates as (#1, £2, %3) = 
(x,y,z) then we have that the components of the above equation are 


Pe = ~~, k=1,2,3. (2.22) 


Just like we defined a position operator <, we now have three position operators (#1, #2, £3) 
making up x. With three position and three momentum operators, we now should state the 
nine possible commutation relations. If you recall our derivation of [%,p] = ih you will note 
that the commutator vanishes unless the superscripts on ~ and p are the same. This means 
that we have 


where the Kronecker delta is defined by 
lifj—4 
unl nee (224) 


In order to write now the general Schrodinger equation we need to consider the kinetic energy 
operator, or the Hamiltonian: 


s p 
A = — t 2.25 
F+V(x,0), (2.25) 
This time h h 
p= pp = 2V-2V = -1?V? (2.26) 
where V? is the Laplacian operator 
OF Or On 
ee ee 2.27 
a Ox? = Oy? = Oz ee 
The Schrodinger equation finally takes the form 
(2.28) 


3 Interpreting the Wavefunction 


Schrodinger thought that the wavefunction V represents a particle that could spread out and 
disintegrate. The fraction of the particle to be found at x would be proportional to the mag- 
nitude of |W|?. This was problematic, as noted by Max Born (1882-1970). Born solved the 
Schrodinger equation for scattering of a particle in a potential, finding a wavefunction that 
fell like 1/r, with r the distance to the scattering center. But Born also noticed that in the 
experiment one does not find fractions of particles going in many directions but rather particles 
remain whole. Born suggested a probabilistic interpretation. In his proposal, 


The wavefunction V(x,t) doesn’t tell us how much of the particle is at position x at 
time t but rather the probability that upon measurement taken at time t we would 
find the particle at position 2. 


To make this precise we use an infinitesimal volume element with volume d?x centered 
around some arbitrary point x. The probability dP to find the particle within the volume 
element d°x at time t is 

dP = |W(x,t)|? d’x. (3.1) 
Consistency requires that the total probability to find the particle somewhere in the whole 
space is unity. Thus the integral of dP over all of space must give one: 


i d’x |W(x,t)|? = 1 (3.2) 


ll space 


Next time we will explore the consistency of this equation with time evolution. 
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1 Normalization and time evolution 


The wavefunction V(x, t) that describes the quantum mechanics of a particle of mass m moving 
in a potential V(z,t) satisfies the Schrédinger equation 


SOW(E,t) Re O° 
or more briefly 
gore. HeGHy (1.2) 


ot 
The interpretation of the wavefunction arises by declaring that dP, defined by 


dP = |W(z,t)|?dz , (1.3) 


is the probability to find the particle in the interval dx centered on x at time t. It follows that 
the probabilities of finding the particle at all possible points must add up to one: 


Z Ww" (et) Vat de= 1. (1.4) 


—Co 


We will try to understand how this equation is compatible with the time evolution prescribed by 
the Schrodinger equation. But before that let us examine what kind of conditions are required 
from wavefunctions in order to satisfy (1.4). 


Suppose the wavefunction has well-defined limits as x — too. If those limits are different 
from zero, the integral around infinity would produce an infinite result, which is inconsistent 
with the claim that the total integral is one. Therefore the limits should be zero: 


lim W(a,t)=0. (1.5) 


L—>xrCO 


It is in principle possible to have a wavefunction that has no well-defined limit at infinity but 
is still is square integrable. But such cases do not seem to appear in practice so we will assume 
that (1.5) holds. It would also be natural to assume that the spatial derivative of UV vanishes as 
x —> +00 but, as we will see soon, it suffices to assume that the limit of the spatial derivative 
of V is bounded 


(1.6) 


We have emphasized before that the overall numerical factor multiplying the wavefunction 
is not physical. But equation (1.4) seems to be in conflict with this: if a given W satisfies it, 
the presumed equivalent 2WV will not! To make precise sense of probabilities it is convenient to 
work with normalized wavefunctions, but it is not necessary, as we show now. Since time plays 
no role in the argument, so assume in all that follows that the equations refer to some time to 
arbitrary but fixed. Suppose you have a wavefunction such that 


[eiwe = Nag i (1.7) 

Then I claim that the probability dP to find the particle in the interval dz about x is given by 
1 

dP = a7 NUP de. (1.8) 


This is consistent because 


_ i on 
pray fem = 5 N= 1. (1.9) 


Note that dP is not changed when WV is multiplied by any number. Thus, this picture makes 
it clear that the overall scale of YW contains no physics. As long as the integral f |W|?dx < oo 
the wavefunction is said to be normalizable, or square-integrable. By adjusting the overall 
coefficient of Y we can then make it normalized. Indeed, again assuming (1.7) the new 
wavefunction W’ defined by 


Ww = — 0, (1.10) 


is properly normalized. Indeed 


fave = x | WwPae = aly (1.11) 


We sometimes work with wavefunctions for which the integral (1.4) is infinite. Such wave- 
functions can be very useful. In fact, the de Broglie plane wave V = exp(ikx — iwt) for a free 
particle is a good example: since |V|? = 1 the integral is in fact infinite. What this means is 
that exp(ikx — iwt) does not truly represent a single particle. To construct a square-integrable 
wavefunction we can use a superposition of plane waves. It is indeed a pleasant surprise that 
the superposition of infinitely many non-square integrable waves is square integrable! 


2 The Wavefunction as a Probability Amplitude 


Let’s begin with a normalized wavefunction at initial time to 


i. W* (x, to) V (a, to)dx = 1. (2.1) 


[oe) 


Since W(x, tg) and the Schrédinger equation determine W for all times, do we then have 


[ W*(2,t)V(x,t)de = 1? (2:2) 


[oe) 


Define the probability density p(z, t) 
plat) = U*(2,1)U(2,t) = We, H/. (2.3) 


Define also V(t) as the integral of the probability density throughout space: 


N(t) = / Baca. (2.4) 
The statement in (2.1) that the wavefunction begins well normalized is 
N (to) = 1, (2.5) 


and the condition that it remain normalized for all later times is V(t) = 1. This would be 
guaranteed if we showed that for all times 
dN (t) 
dt 


=0. (2.6) 


We call this conservation of probability. Let’s check if the Schrodinger equation ensures this 
condition will hold: 


dN(t) _ i. Op(x,t) 1, 


dt a OF Q.7) 
eK OU* OW (zx, t) 
= W(2,t) + U* ——_—— : 
[(Feveo+ wen ) ae 
From the Schrodinger equation, and its complex conjugate 
OW Ow Wd: 
_ Ow tee OW train, 
—ih ae (AWD) = ae oe (HW)*. (2.9) 


In complex conjugating the Schrodinger equation we used that the complex conjugate of the 
time derivative of WV is simply the time derivative of the complex conjugate of V. To conjugate 


3 


the right hand side we simply added the star to the whole of HV. We now use (2.8) and (2.9) 
in (2.7) to find 


aN ys [ (nw w - iY (HY) ae 


dt ioe Nell 
a a (2.10) 
7 -( fl (AW)* Udo — | w"(HV)dr) 
To show that the time derivative of N(t) vanishes, it suffices to show that 
(2.11) 


Equation (2.11) is the condition on the Hamiltonian operator H for conservation of probability. 
In fact, if H is a Hermitian operator the condition will be satisfied. The operator H is a 
Hermitian operator if it satisfies 


Hermitian operator: i (HV,)*U, = / U*(AWV). (2:12) 


[oe) 


[o.e) 


Here we have two wavefunctions that are arbitrary, but satisfy the conditions (1.5) and (1.6). 
As you can see, a Hermitian operator can be switched from acting on the first function to acting 
on the second function. When the two functions are the same, we recover condition (2.11). 

It is worth closing this circle of ideas by defining the Hermitian conjugate J" of the linear 
operator T’. This is done as follows: 


ar (TV,) = (PHY, Wy. (2.13) 
j j 


Co [oe) 


The operator TJ’, which is also linear, is calculated by starting from the left-hand side and 
trying to recast the expression with no operator acting on the second function. An operator T 
is said to be Hermitian if it is equal to its Hermitian conjugate: 


T is Hermitian if TP. (2.14) 


Hermitian operators are very important in Quantum Mechanics. They have real eigenvalues 
and one can always find a basis of the state space in terms of orthonormal eigenstates. It turns 
out that observables in Quantum Mechanics are represented by Hermitian operators, and the 
possible measured values of those observables are given by their eigenvalues. Our quest to show 
that normalization is preserved under time evolution in Quantum Mechanics has come down 
to showing that the Hamiltonian operator is Hermitian. 


3 The Probability Current 


Let’s take a closer look at the integrand of equation (2.10). Using the explicit expression for 
the Hamiltonian we have 


OP _ Law) w — W*(ATW) 
ot h (3.1) 
i flow OV 7 
The contributions from the potential cancel and we then get 
oem - h (or ru 
—((HV)*U—wW*(HV)) = — Vv —-v*— }. ; 
(avy v- ww) = = (Fov-w Fs) (3.2) 


The only chance to get to show that the integral of the right-hand side is zero is to show that 
it is a total derivative. Indeed, it is! 


Papa ite yee. 20 fou ow 
(Hee 0Ey)) = slam (ae x) 
OT h (.,0v av 
‘ -slaa(¥ ae oe) | 
OTR. ow 


O|h Ow 
= ——}]—Im( v*—— 
Ox - m( Ox )] 
where we used that z — 2* = 2iIm(z). Recall that the left-hand side we have evaluated is 
actually oe and therefore the result obtained so far is 


Op. OO NN Ow 
4 |] —Im[( w*2— — 0. A 
+o |* m( =) ee 
This equation encodes charge conservation and is of the type 
Op OJ 
—s —— a 0 . 
a Be 2) 


where J(x,t) is the current associated with the charge density p. We have therefore identified 
a probability current 


(3.6) 


There is just one component for this current since the particle moves in one dimension. The 
units of J are one over time, or probability per unit time, as we now verify. 


For one spatial dimension, [V] = L~'/?, which is easily seen from the requirement that 
J dz|W¥/? is unit free. (When working with d spatial dimensions the wavefunction will have 
units of L~4/?). We then have 


OV 1 ML? h ie 
or el=a W=-=F. [E)=5. (3.7) 
1 ee One 
> [J] ri probability per unit time (3.8) 
We can now show that the time derivative of N is zero. Indeed, using (3.5) we have 
dN Con Op oar 
en fee ~~ f Be = Ue.) - J-00,9) (3.9) 
The derivative vanishes if the probability current vanishes at infinity. Recalling that 
h Ow Ow* 
= ——( V*—_ — V—_ 3.10 
a ( Ox Ox ) om) 


we see that the current indeed vanishes because we restrict ourselves to wavefunctions for which 


Titra Oven Timi a ou remains bounded. We therefore have 


(3.11) 


as we wanted to show. 


To illustrate how probability conservation works more generally in one dimension, focus 
on a segment x € |a,b]. Then the probability P,, to find the particle in the segment |a, }], is 
given by 


b 
Pe if p(x, t) dx. (3.12) 
If we now take the time derivative of this and, as before, use current conservation we get 


dPip a OJ (x,t) 
a Pees 


dt = —J(b,t) + J(a,t). (3.13) 


This is the expected result. If the amount of probability in the region [a, b] changes in time, it 
must be due to probability current flowing in or out at the edges of the interval. Assuming the 
currents at x = 6b and at x = a are positive, we note that probability is flowing out at x = b 
and is coming in at x = a. The signs in the above right-hand side correctly reflect the effect of 
these flows on the rate of change of the total probability inside the segment. 


4 Probability current in 3D and current conservation 


The determination of the probability current J for a particle moving in three dimensions follows 
the route taken before, but we use the 3D version of the Schrdédinger equation. After some work 
(homework) the probability density and the current are determined to be 


(4.1) 
(4.2) 
In three spatial dimensions, [V] = L~? and the units of J are quickly determined 
1 h I? 

WV) = — —|=— 4, 

wu=7. [2|=5 (4.3) 
1 a oe : 

> [J] TE probability per unit time per unit area (4.4) 


The conservation equation (4.2) is particularly clear in integral language. Consider a fixed 
region V of space and the probability Qy(t) to find the particle inside the region: 


Quit) = f olxpa'x. (4.5) 
V 
The time derivative of the probability is then calculated using the conservation equation 
dQv Op 33 3 
ay —d’x = —]| V-Jd’x. 4.6 
di vot - - Ge) 


Finally, using Gauss’ law we find 


dQy 


where S is the boundary of the volume V. The interpretation here is clear. The probability 
that the particle is inside V may change in time if there is flux of the probability current across 
the boundary of the region. When the volume extends throughout space, the boundary is at 
infinity, and the conditions on the wavefunction (which we have not discussed in the 3D case) 
imply that the flux across the boundary at infinity vanishes. 


Our probability density, probability current, and current conservation are in perfect analogy 
to electromagnetic charge density, current density, and current conservation. In electromag- 
netism charges flow, in quantum mechanics probability flows. The terms of the correspondence 
are summarized by the table. 


a Electromagnetism Quantum Mechanics 
charge density probability density 


charge in a volume V | probability to find particle in V 
current density probability current density 
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1 Wavepackets and Uncertainty 


ka 


A wavepacket is a superposition of plane waves e’** with various wavelengths. Let us work with 


wavepackets at ¢t = 0. Such a wavepacket is of the form 
1 = . 
W(x, 0) = —— &(k)e"*"dk. Li 
(2,0) == [ #(8) (1.1 


If we know W(z,0) then ®(k) is calculable. In fact, by the Fourier inversion theorem, the function 
®(k) is the Fourier transform of U(zx,0), so we can write 


B(k) = oe / * W(x, 0)e** de. (1.2) 


Note the symmetry in the two equations above. Our goal here is to understand how the uncertainties 
in U(z,0) and ®(k) are related. In the quantum mechanical interpretation of the above equations one 
recalls that a plane wave with momentum hk is of the form e’**. Thus the Fourier representation of 
the wave V(z,0) gives the way to represent the wave as a superposition of plane waves of different 


momenta. 


‘Ap 


Figure 1: A ®(k) that is centered about k = ko and has width Ak. 


Let us consider a positive-definite ®(k) that is real, symmetric about a maximum at k = ko, and has 
a width or uncertainty Ak, as shown in Fig. 1. The resulting wavefunction U(ax,0) is centered around 
x = 0. This follows directly from the stationary phase argument applied to (1.1). The wavefunction 
will have some width Az, as shown in Fig. 2. Note that we are plotting there the absolute value 
|W(x,0)| of the wave packet. Since U(x,0) is complex, the other option would have been to plot the 
real and imaginary parts of U(x, 0). 


“7 
x 


Figure 2: The U(z,0) corresponding to ®(k) shown in Fig. 1, centered around x = 0 with width Az. 


Indeed, in our case V(x, 0) is not real! We can show that 


W(x, 0) is real if and only if ®*(—k) = ®(k). (1.3) 


Begin by complex conjugating the expression (1.1) for U(x, 0): 


U*(z,0) = = / ” O*(k)e**dk = = / ” o*(—k)e* dk. (1.4) 


In the second step we let k + —k in the integral, which is allowable because we are integrating over 
all k, and the two sign flips, one from the measure dk and one from switching the limits of integration, 
cancel each other out. If 6*(—k) = ®(k) then 
1 v ik 
U*(x,0) = — O(k)e" "dk = V(z,0), 1.5 
(2,0) == [aa (0,0) (1.5) 


as we wanted to check. If, on the other hand we know that V(z,0) is real then the equality of U* and 


=f. &*(—k)edk = sf B(k)ei**dk. (1.6) 


i (®*(—k) — ®(k)) dk = 0. (1.7) 


This equation actually means that the object over the brace must vanish. Indeed, the integral is 


W gives 


This is equivalent to 


computing the Fourier transform of the object with the brace, and it tells us that it is zero. But a 
function with zero Fourier transform must be zero itself (by the Fourier theorem). Therefore reality 
implies ®*(—k) = ®(k), as we wanted to show. 

The condition 6*(—k) = ®(k) implies that whenever ® is non-zero for some k it must also be non- 
zero for —k. This is not true for our chosen ®(k): there a bump around ko but is no corresponding 
bump around —ko. Therefore U(x,0) is not real and U(2x,0) will have both a real and an imaginary 
part, both centered on x = 0, as shown in Fig. 3. 

Let’s now get to the issue of width. Consider the integral for (2, 0) 


1 o tka 
a. (ke dk, (1.8) 


and change variable of integration by letting k = kp + k, where the new variable of integration k 


(7,0) = 


parameterizes distance to the peak in the momentum distribution. We then have 


U(x,0) = 7 | ® (ky +k) elke i, (1.9) 


2 


Figure 3: The real and imaginary parts of U(x, 0). 


As we integrate over k, the most relevant region is 
T. Ak Ak 
ke laos | (1.10) 


because this is where ®(k) is large. As we sweep this region, the phase kx in the exponential varies in 
the interval 7 
kee [-422, AR | (for x > 0), (1.11) 


and the total phase excursion is Akz. We will get a substantial contribution to the integral for a 
small total phase excursion; if the excursion is large, the integral will get washed out. Thus, we get a 
significant contribution for Ak|x| < 1, and have cancelling contributions for Ak|z| > 1. 

From this, we conclude that (2,0) will be nonzero for x € (—20,20) where xo is a constant for 
which Akzo + 1. We identify the width of U(x,0) with Ax := 2x9 and therefore we have Ak 5A x1. 
Since factors of two are clearly unreliable in this argument, we simply record 


Az Ak = 1. (1.12) 


This is what we wanted to show. The product of the uncertainty in the momentum distribution and 
the uncertainty in the position is a constant of order one. This uncertainty product is not quantum 
mechanical; as you have seen, it follows from properties of Fourier transforms. 

The quantum mechanical input appears when we identify hk as the momentum p. This identifica- 
tion allows us to relate momentum and k uncertainties: 


Ap = hAk. (1.13) 
As a result, we can multiply equation (1.12) by h to find: 
Ag Ap & fh. (1.14) 


This is the rough version of the Heisenberg uncertainty product. The precise version requires definining 
Az and Ap precisely. One can then show that 


Heisenberg uncertainty product: Az Ap > ; (1.15) 


The product of uncertainties has a lower bound. 


Example: Consider the case where ®(k) is a finite step of width Ak and height 1// Ak, as shown in 
Fig. 4. Find (x,t) and estimate the value of Az. 


Figure 4: A momentum distribution. 


Note that the U(z,0) we aim to calculate should be real because &*(—k) 


®(k). From the integral 
representation, 
Ak 
1 a 1 : 
W(2z,0) = / ett dk 
aM) Van Jak VAk 
1 etka a 
~ 4/ tx 
ann = (1.16) 
1 ee — ei 
V27Ak wx 
- Ak 
Lay ate _ [AE sin 
x 
2rAk 20 a 
We display U(2,0) in Fig. 5. We estimate 


AX Ab Ab AY me 


Figure 5: The Y(x,0) corresponding to ®(k). 


27 
Ate Ty > AcAks 2r. 


(1.17) 
2 Wavepacket Shape Changes 


In order to appreciate general features of the motion of a wave-packet we looked at the general solution 
of the Schrédinger equation 


U(a,t) = -_ 7 B(k)etlAe—W(h)) dk, (2.18) 


4 


and under the assumption that ®(k) peaks around some value k = kg we expanded the frequency w(k) 
in a Taylor expansion around k = ko. Keeping terms up to and including (k — ko)? we have 
dw 1 » d’w 


w(k) = w(ko) + (k — ko) dk, + 5 (k — ko) dk? |, 


(2.19) 


0 


The second term played a role in the determination of the group velocity and the next term, with 
second derivatives of w is responsible for the shape distorsion that occurs as time goes by. The 
derivatives are promptly evaluated, 

dw dE p __hk aw h 


= = ; 2.2 
dk dp m ™m’ dk? m ey) 


Since all higher derivatives vanish, the expansion in (2.19) is actually exact as written. What kind of 
phase contribution are we neglecting when we ignore the last term in (2.19)? We have 


en tw(k)t ps 13 (kk)? At (2.21) 


Assume we start with the packet at t = 0 and evolve in time to t > 0. This phase will be ignorable as 
long as its magnitude is significantly less than one: 


h 
— Deh: 
(k— ko) —t <1. (2.22) 


We can estimate (k — kg)? ~ (Ak)? since the relevant k values must be within the width of the 
momentum distribution. Moreover since Ap = hAk we get 


1, (2.23) 


Thus, the condition for minimal shape change is 


t< at ; (2.24) 


We can express the inequality in terms of position uncertainty using ArAp = h. We then get 


t< (Ae)? (2.25) 
Also from (2.24) we can write 

Apt es h (2.26) 

m Ap’ : 
which gives 

A 

Pt < Az. (2.27) 

m 


This inequality has a clear interpretation. First note that Ap/m represents the uncertainty in the 
velocity of the packet. There will be shape change when this velocity uncertainty through time 
produces position uncertainties comparable to the width Az of the wave packet. 


In all of the above inequalities we use < and this gives the condition for negligible shape change. 
If we replace < by © we are giving an estimate for some measurable change in shape. 


Exercise: Assume we have localized an electron down to Ax = 107!°m. Estimate the maximum time 
t that it may remain localized to that level. 
Using (2.25) we have 

m(Az)? mc? (Az)? 0.5 MeV - 1077°m? 


x = = 10-8, 2.28 
: h he-c 200 MeVim - 3 x 108m /s ‘7 ee) 


If we originally had Av = 10~?m, we would have gotten t = Is! 


3. Time evolution of a free wave packet 


Suppose you know the wavefunction U(x,0) at time equal zero. Your goal is finding U(z,t). This is 
accomplished in a few simple steps. 


1. Use U(x,0) to compute ®(k): 
1 xs . 
®(k) = = da U(x, 0)e**. (3.1) 


2. Use ®(k) to rewrite U(x,0) as a superposition of plane waves: 


(2,0) = om / * o(Re* dk. (3.2) 


This is useful because we know how plane waves evolve in time. The above is called the Fourier 
representation of U(x, 0). 


3. A plane wave ~;(ax,0) = e*** evolves in time into ~,(a,t) = e***-#)) with hw = aay Using 
superposition we have that (3.2) evolves into 


ey po 
U(x,t) = = vi B(k) el he“) dk, (3.3) 


This is in fact the answer for (x,t). One can easily confirm this is the solution because: (i) it 


solves the Schrédinger equation (check that!) and (ii) setting t = 0 in U(z,t) gives us the initial 
wavefunction (3.2) that represented the initial condition. 


4. If possible, do the integral over k to find a closed form expression for U(z,t). If too hard, the 
integral can always be done numerically. 


Example: Evolution of a free Gaussian wave packet. Take 


1 2 [a2 
Pa(z,0) = Ont Ja ° ame (3.4) 


This is a gaussian wave packet at t = 0. The constant a has units of length and Ar = a. The state 
Wa is properly normalized, as you can check that f dx|Wa(x,0)|? = 1. 


We will not do the calculations here, but we can imagine that this packet will change shape as 
time evolves. What is the time scale 7 for shape changes? Equation (2.25) gives us a clue. The right 
hand side represents a time scale for change of shape. So we must have 

HLS 2. 
TS a. 3.5 
; (3.5) 
This is in fact right. You will discover when evolving the Gaussian that the relevant time interval is 


actually just twice the above time: 
2 


2ma 
= 3.6 
r= (3.6) 
If we consider the norm-squared of the wavefunction 
1 1 2/902 
w* 0 2 = Lag oe /2a 3.7 
Pe a, 0) = are le, (3.7) 
you will find that after time evolution one has 
1 1 2Jo,2 
w* it 2 Gs ei ae /2a (t) 3.8 
We@e.OP = aaa (3.8) 


where a(t) is a time-dependent width. The goal of your calculation will be to determine a(t) and to 
see how 7 enters in a(t). 
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1 Uncovering momentum space 


We now begin a series of developments that lead to the idea of momentum space as a counterpoint 
or dual of position space. In this section the time dependence of wavefunctions will play no role. 
Therefore we will simply suppress time dependence. You can imagine all wavefunctions evaluated at 
time equal zero or at some arbitrary time to. 

We begin by recalling the key identities of Fourier’s theorem: 


= = i elke 
U(x) = a: B(k)e* dk, 

= fh. - x e tke xc 
w= fH) ae 


The Fourier transform ®(k) has all the information carried by the wavefunction U(x). This is clear 


(1.1) 


because knowing ®(k) means knowing V(2). The function ®(k) also acts as the weight with which we 
add the plane waves with momentum fk to form U(x). 

We will now see that the consistency of the above equations can be used to derive an integral 
representation for a delta function. Such a representation is a needed tool for our upcoming discussion. 
The idea is to replace ®(k) in the first equation by the value given in the second equation. In order 
to keep the notation clear, we must use x’ as a dummy variable of integration in the second equation. 


We have : xo / a 
U(x) = — dk ce | dx! ek? Ue! 


_ = ! ! =|. ik(a—2’) 
=f a UV (2) oa Dee : 
a 


Look at the type of integral. The factor indicated by the brace happens to reduce the 2’ integral to 
an evaluation at x. We know that 6(x’ — x) is the function such that for general f(x), 


(1.2) 


[a fe) se! - 2) = Fe), (1.3) 
and so we conclude that the factor indicated by the brace is a delta function 


O(a’ —2)= = | deere?) (1.4) 


In this integral one can let k + —k and since f dk is left-invariant under this replacement, we find 
that 6(a’ — x) = 6(a — 2’), or more plainly 6(a) = 6(—x). We will record the integral representation 
of the delta function using the other sign: 


Another useful property of delta functions is 
d(ax) = —d(a). (1.6) 


At this point we ask: How does the normalization condition for V(x) look in terms of ®(k)? We 
must simply calculate. We have 


[_@ve@ua- [wef emema se fame. 07 


Rearranging the integrals to do the x integral first we write 


| dx U* (x) U(x eee af dk! ®*( ef dx ei —k)x 
aes oz 


=e dk dk! ®*(k)®(k')5(k! — k) (1.8) 


= ee dk ®*(k)®(k) , 


where we recognized the presence of a delta function and we did the integral over k’. Our final result 
is therefore 


(1.9) 


This is known as Parseval’s theorem, or more generally, Plancherel’s theorem. This equation relates 
the U(x) normalization to a rather analogous normalization for ®(k). This is a hint that just like for 
|W(a)|?, we may have a probability interpretation for |®(k)|?. 

Since physically we associate our plane waves with eigenstates of momentum, let us rewrite Par- 
seval’s theorem using momentum p = hk. Instead of integrals over k we will have integrals over p. 
Letting ®(p) = ®(k) equations (1.1) become 


1 aie ipx 
Va) = 5 | Boerap, 


(1.10) 


®(p) = 


Tela 


For a more symmetric pair of equations we can redefine the function ®(p). We will let ®(p) > O(p)Vh 


Je Pe! hy . 


in equations (1.10). We then obtain our final form for Fourier’s relations in terms of momentum: 


V(x) -—= | ) etPe/h dy , 


(1.11) 


1 ms —ipx 


Similarly, Parseval’s theorem (1.9) becomes 


(1.12) 


Exercise. Verify that the redefinitions we did to arrive at (1.11) indeed yield (1.12) when starting 
from (1.9). 


Our interpretation of the top equation in (1.11) is that ®(p) denotes the weight with which we add 
the momentum state e‘?*/" in the superposition that represents U(x). This momentum state eipa/h 
is an eigenstate of the momentum operator p with eigenvalue p. Just like we say that W(x) is the 
wavefunction in position space x, we can think of ®(p) as the wavefunction in momentum space p. 
The Parseval identity (1.12) suggests that ®(p) has a probabilistic interpretation as well. Given that 
a properly normalized W(x) leads to a ®(p) that satisfies [ dp|®(p)|? = 1, we postulate that: 


|©(p)|?dp is the probability to find the particle with momentum in the range (p,p + dp). | (1.13) 


This makes the analogy between position and momentum space quite complete. 


Let’s consider the generalization to 3D. Fourier’s theorem in momentum space language (namely, using 
p as opposed to k) takes the form 


1 - ip-x 
VO) = aoe | ap 8p)e len 


(1.14) 


1 e —ip-x 
®(p) = ae d°x W(x)e P [Re 


Just like we did in the 1D case, if we insert the Fourier transform into the expression for U(x), we 
find an integral representation for the 3D 6-function 


1 ; ee, 
U(x) oat Can 5 /* ie U(x’ )eP* /h 
35/ p:(x—x’)/h 
= f ax’ v ‘a sf@ pe'P Ga15} 


a [ex U(x’ 


which leads to the identification 


k ude a 


P(x —x’) = oo / Pers): (1.16) 


It is then straightforward to derive Parseval’s identity (exercise!). We find 


[. d?x |U(x > = | ap |e)? (1.17) 


2 dp is the probability to 


We use in 3D momentum space the same probability interpretation: |®(p) 
find the particle with momentum in the range d°p around p. 


2 Expectation Values of Operators 


Consider a random variable Q that takes values in the set {Qi,...,Qn} with respective probabilities 
{pi,.--;Pn}. The expectation value (Q) (or expected value) of Q is the average value that we expect 
to find after repeated observation of Q, and is given by the formula 


(Q) = S- Qipi. (2.18) 
i=1 


As we have seen, in a quantum system the probability for a particle to be found in [x, x + dz] at 
time t is given by 


W* (x, t)U (a, t) dx. (2.19) 
Thus, the expected value of x, denoted as (%) is given by 
is) = Ui 2 U" (x, t)U(a, t)de. (2.20) 


Note that this expected value depends on t. What does (%) correspond to physically? If we consider 
many copies of the physical system, and measure the position x at a time ¢ in all of them, then the 
average value recorded will converge to (%) as the number of measurements approaches infinity. 


Let’s discuss now the expectation value for the momentum. Since we have stated that 
®*(p, t)®(p, t) dp (2.21) 


is the probability to find the particle with momentum in the range [p, p + dp] at time t, we define the 
expectation (jp) of the momentum operator as 


On i p®*(p,t)®(p,t)dp. (2.22) 


We will now manipulate this expression to see what form it takes in coordinate space. Using (1.11) 
and its complex conjugate version we have 


[oe) [o-e) : [oe / a 7 
= i, app | ae elpe/h w(c.t) | eae e (Pe MY (a! ,t) 
=o —o V2Th 


oo V2Th 


w= | * yp &*(p,t)®(p,t) dp 


1 sd * me = ipxz/h ,—ipa’ /h 
=sy [av (ot) fact veal f apperle pe ie (2.23) 
= 1 ee * + / / iS ho ipx/h ,—ipax’ /h 
=a f/f #¥ (ot) fae wast) fap (Fz e e 


_ * ! ! ig ipx/h ipa! /h 
[ww ey [a wen (F2) ay fae e : 


Letting p = hu in the final integral we have 
1 (oe) : /h . ‘Ih 1 (oe) : ( 1) ; 
d wpxr ee) a A aa = _ : 2.24 
2rh [. ue . 20 i- ue eg, ae 


4 


As a result, we have 


(P) 
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a 
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(2.25) 
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where we changed the order of integration. The z’ integral is now easily done and we find 


= [aon (FZ) wen), (2.26) 


We have thus shown that 


ho 


5 = de W(0,1)pW(0,t), P= Ae. (2.27) 


Notice the position of the p operator: it acts on V(x). This motivates the following definition for the 


expectation value (Q) of any operator Q: 


Q)= [de ve.) Qu(a,2). (2.28) 


Example: Consider the kinetic energy operator T fora particle moving in 1D: 


AD 2 92 
fos ee 
= 2n 2m Ox?’ Ce 
The definition gives 
‘ h2 . 0? 
P) = f de Ve.) g5 Ve). (2.30) 


The kinetic energy is a positive operator (being proportional to the square of the momentum operator). 
It is therefore of interest to make this positivity manifest. Integrating by parts one of the x derivatives 
and ignoring boundary terms that are presumed to vanish, we find 


(t) = 


2 2 
eee ie AO (2.31) 
2m 


x 


This is manifestly positive! The expectation value of T can also be computed in momentum space 
using the probabilistic interpretation that led to (2.22): 


= fo 2 
= | pF lap,0?, (2.32) 


Other examples of operators whose expectation values we can now compute are the momentum 
operator p > Ay in 3D, the potential energy operator, V(x), and the angular momentum operator 


rx p= (Yb. — 2py , ZPy — £Dz , EPpy — YPx) 
a a a a a Fs) 2) (2.33) 
=-|y—-2z = : 

4 Oz 


ay? Oa "Oz Oy "a 


3 Time dependence of expectation values 


The expectation values of operators are in general time dependent because the wavefunctions repre- 
senting the states are time dependent. We will consider here operators that do not have explicit time 
dependence, that is, operators that do no : 


ih—(Q) = ine [. Ba U* (x, t)QU(x,t) 


dt 
=in [ Bar (5 ous woe) 


ns ot 
a (3.34) 
en a ak pent 
= in f dbx (jury ow = 5V Quy) 
= / hie (worry = (T¥")QU) 
We now recall the hermiticity of H, which implies that 
i, dx (HW,)*Vy = / dx Vi HW. (3.35) 


This can be applied to the second term in the last right-hand side of (3.34) to move H into the other 
wavefunction Be 
ine (Q) = | Ber (wony = wHQW) 
Res (3.36) 
=f dev [QA] y, 


where we noted the appearance of the commutator. All in all, we have proven that for operators Q 
that do not explicitly depend on time, 


(3.37) 
Note that the commutator satisfies the following properties (homework): 
[A, A] = 0 (3.39) 
[A,B+C] =[A, B] +[A,C] (3.40) 
[A, BC] = [A, B]C + B[A,C] (3.41) 
[AB,C] = A[B,C]+[A,C]B (3.42) 
0 = [A, [B, C]] + [B, [C, A] + [C, [A, B]. (3.43) 
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1 Observables and Hermitian operators 


Let’s begin by recalling the definition of a Hermitian operator. The operator Q is Hermitian if for the 
class of wavefunctions UW we work with, 


fu Us QU, = fu (QU,)* Vo. (1.1) 


We will sometimes use a briefer notation for the integrals of pairs of functions: 


(1.2) 
Note that for any constant a 
(aW,, V2) a a* (Wy, Wo), (U1, aW2) — a(W,, V2). (Gigs) 
With this notation the condition of Hermiticity is more briefly stated 
Q is Hermitian: (W1,QU2) = (QU, WV). (1.4) 
The expectation value of Q was defined by 
Q)v = fu wQu = (VW, QU). (1.5) 


For this formula to make sense, the state WV must be normalized. 


Claim 1. The expectation value of a Hermitian operator is real. To prove this we complex conjugate 
the above definition. The complex conjugate of the integral is the integral of the complex conjugate 
of the integrand, therefore 


((Q)w) "= f aol (u* Qu)" = fa U(Qu)* = f a (QU)* (1.6) 


Note that Qu is a wavefunction, so it makes sense to take its complex conjugate (we never have to 
think of conjugating Q). Using the Hermiticity of the operator, we move it into W to get 


(ix) = / dex VOW = (Q)u, (1.7) 


1 


thus showing that the expectation value is indeed real. 


Claim 2. The eigenvalues of a Hermitian operator are real. Assume the operator Q has an eigenvalue 
qi associated with a normalized eigenfunction ~1(z): 


Qvi(x) = aii(z). (1.8) 


Now compute the expectation value of Q in the state yj: 


(Q)u. = (41, Qv1) = (i avi) = ati, v1) =4- (1.9) 


By claim 1, the expectation value is real, and so is the eigenvalue q,, as we wanted to show. Note the 
interesting fact that the expectation value of Q on an eigenstate is precisely given by the corresponding 
eigenvalue. 


Consider now the collection of eigenfunctions and eigenvalues of the Hermitian operator Q: 


Qvi(x) = nvi(2), (1.10) 
Q v2(x) = q2v2(a), (1.11) 


The list may be finite or infinite. 


Claim 3. The eigenfunctions can be organized to satisfy orthonormality: 


Cone / dar Wt (a)abj(a) = by. (1.12) 


For 7 = j, this is just a matter of normalizing properly each eigenfunction, which we can easily do. 
The equation also states that different eigenfunctions are orthogonal, or have zero overlap. We now 
explain why this is so for i 4 j with q; # q;. Indeed, for this we evaluate (wi, Qv;) in two different 
ways. First 


(bi,Qv;) = (Wiads) = Gis dy), (1.13) 


and second, using Hermiticity of Q, and the reality of the eigenvalues 
(bi, Qv;) = (Qviss) = Gidity) = Gis). (1.14) 
Equating the final right-hand sides in the two evaluations we get 


(a5 — 4) (di, 3) = 0. (1.15) 


Since the eigenvalues were assumed different, this proves that (qj, ~);) = 0, as claimed. This is not yet 
a full proof of (1.12) because it is possible to have degeneracies in the spectrum, namely, different 
eigenfunctions with the same eigenvalue. In that case the above argument does not work. One must 
then show that it is possible to choose linear combinations of the degenerate eigenfunctions that 
are mutually orthogonal (the orthogonality with the eigenfunctions beyond the degenerate space is 
automatic). This is done in 8.05! 


Claim 4. The eigenfunctions of Q for a complete set of basis functions. Any reasonable UV can be 
written as a superposition of Q eigenfunctions. (This is the so-called spectral theorem, which is proven 
in 8.05 in the finite-dimensional case.) This means that 


U(x) = ayy (x) + agye(x (1.16) 


with calculable coefficients a;. Indeed, if we know the eigenfunctions we have that a; is calculated 
doing the integral of w; against WU: 


We quickly prove this by doing the integral 
fu w;(x)U(e2) = fu wr (x) S > ayh;(2) = Lei fu wy (x); (x = 20h ig Se (118) 
j 
The condition that W is normalized implies a condition on the coefficients a;. We have 
fev me = [de Davia aio 
a j 
= Dafa; fu bj (a); (2) (1.19) 


II 
M 
a 
< 
= 
&. 
| 
i 
2 
£ 


so that the normalization of UV implies that 


So lei? = 1 (1.20) 


a 


We are finally in the position to state the measurement postulate. This is the way in which we 
understand that Hermitian operators represent observables and learn the rules that they follow. 


Measurement Postulate: If we measure the Hermitian operator Q in the state VU, the possible 
outcomes for the measurement are the eigenvalues qi, q2,.... The probability p; to measure q; is given by 


pe = lal’, (1.21) 
where U(x) = >>, aiy;(x). After the outcome q;, the state of the system becomes 
Wa) Sa (Be): (1.22) 


This is called the collapse of the wavefunction. 


The collapse of the wavefunction implies that immediately after the measurement that yielded q; 
a repeated measurement of Q will yield g; with no uncertainty. A small subtlety occurs if we have 
degenerate eigenstates. Suppose the wavefunction contains a piece 


VU = (av; + arr) Seteces (1.23) 


where yw; and y; happen to have the same eigenvalue g and the dots represent other terms. Then if 
we measure gq the state after the measurement collapses to the sum of those two terms 


yp = TORR (1.24) 


Vaal? + Tag?’ 


with the square root denominator included to provide the proper normalization to V. As a consistency 
check note that the probabilities p; to find the various eigenvalues as outcomes properly add to one: 


Dm = Delo? = 1, (1.25) 


) 


by the normalization condition for VY given in (1.20). The measurement postulate follows the Copen- 
hagen interpretation of quantum mechanics. 

Note that the measurement postulate uses the property that any vector in a vector space can be 
written as a sum of different vectors in an infinite number of ways. If we are to measure Q we expand 
the state in Qj eigenstates, if we are to measure Q2 we we expand the state in Q eigenstates, and so 
on and so forth. Each decomposition is suitable for a particular measurement. Each decomposition 
reveals the various probabilities for the outcomes of the specific observable. 


Exercise: Use the expansion V = 5°; a;7); to compute the expectation value (Q). We find 
Q) = fae Safer QS ajvy(o) 
a J 


= Naja; [ de vf(o)Qv(0) 
id (1.26) 
= Ss" OF a54q; fu Vi (@) oj (2) 
ig 


7 ye 050595015 = Ss" lai|?a: = Sidi. 
a,j a 4 


This matches our expectations: the expectation value of Q is the sum of the possible outcomes q; 
multiplied by the corresponding probabilities p;. This is a nice consistency check on our definition of 
expectation values. 


Example. Free particle on the circle x € |0, L}. 
We imagine the points x = 0 and x = L identified to form a circle of circumference L. A wavefunction 
W(x) on the circle must satisfy the periodicity condition 


U(x +L) = V(x), (1.27) 


Assume that at some fixed time we have the wavefunction 


Wie) = BQwe + [2 cos S82), (1.28) 


This wavefunction satisfies the periodicity condition, as you should check. We want to know what are 
the possible values of the momentum and their corresponding probabilities. 


Given our discussion, we must find the set of momentum eigenstates and rewrite the wavefunction 


as a superposition of such states. Momentum eigenstates are exponentials of the form e’**. Two 


things happen on the circle that do not happen in free space. First the momentum will be quantized 


as a consequence of the periodicity condition (1.27). Second, since space here is of finite length, the 
momentum wavefunctions will be normalizable. Consider the periodicity condition as applied to e’**. 
We need 


elke — eik(et+h) _, eth 4 _» kL=2nm, meZ. (1.29) 


Note that m can be any integer, positive, negative, or zero. We thus write for the momentum eigen- 
states, labeled by m 


2rima 


Wale) = Ne Es (1.30) 


with N areal normalization constant. The normalization condition gives 


L L 
ite ‘| eee aie = Nef ieee Ne Ss N=. (1.31) 
Therefore our momentum eigenstates are 
1 anime 
One) = VE’ oe (1.32) 
and these are states with momentum p,,, which is calculated as follows 
Bm = tom = Oe om = (1.33) 


Now that we are equipped with the momentum eigenstates we must simply rewrite the wavefunction 
(1.28) as a superposition of such states: 


1 1 6rixz 


eh a —eE nee — a e ? ‘ 
327 ./L Jf/32/7L 


We then recognize that we have 


Wea) = agile) gq he) + Tvs) + vale). (1.35) 


This is our key result: the original wavefunction written as a superposition of momentum eigenstates 
Um(x). We can now give the possible values p of the momentum and their corresponding probabili- 


ties P: 
2 2 
pa 2th pk ie L 
L 3 27 6 
Oh SL 1 
aT Pa ea) = 
: (1.36) 
6h 2) 1 
p= rar P= 5| SS 
L 3 3 
6rh 1/2 1 
SS. P= =| = -, 
P 7 it 3 


2 Uncertainty 


For random variables, the uncertainty is the standard deviation—the square root of the expected value 
of the square of deviations from the average value. Let Q be a random variable that takes on values 
Q1,---;Qn with probabilities p1,..., pn, respectively. The expected value is 


O-= S70. (2.37) 
and the variance (the square of the standard deviation) is 
(AQ) = > pi(Q:-Q). (2.38) 


This definition makes it clear that if AQ = 0, then the random variable is constant: each term in the 
above sum must vanish, making Q; = Q, for all 7. We find another useful expression by expanding 


the above definition = = 
(AQ)? = S> iQ? -25¢ viQiQ+ >— viQ 


(2.39) 
=@ -200+Q =@-@, 
where we used )°, p; = 1. Therefore 
ae 
(AQ? = @-Q. (2.40) 
Since, by definition (AQ)? > 0, we have the interesting inequality 
Q@>Q. (2.41) 


Now let us consider the quantum mechanical case. We have already defined expectation values of 
Hermitian operators, so we can now mimic the definition (2.40) and declare that the uncertainty AQw 
of an operator in a state V is a real number whose square is given by 


(AQ)G = (Q*)w — ((Q)w)?. (2.42) 


Sometimes for brevity, we omit the label for the state, 
(AQ)? = (Q*) — (Q)’. (2.43) 


Claim 1. The uncertainty can also be written as the expectation value of the square of the difference 
between the operator and its expectation value: 


(AQ)? = ((@- (Q))’). (2.44) 


Indeed, expanding the right hand side we have 


(Q? — 2Q(Q) + (Q)”) = (Q?) — 2(Q)(Q) + (Q)?_ = (Q) — (Q)’. (2.45) 


Claim 2. The uncertainty can be written as the integral of the square of the norm of (Q — (Q))¥: 
2 o z 2 
(ag? = f a|(@-@) we) (2.46) 


—Co 


Indeed to prove this we begin with (2.44) By a very similar proof, we can show this is equivalent to 
F i ae 2 
(AQ? = ((Q-(Q))*) = far ¥*(Q-(Q)) v. (2.47) 
Using the Hermiticity of Q and the reality of (Q) we can bring one of the two factors to act on W*: 


(AQ? = f ae [(O-(@)¥] @- (Qn = fae |(Q-(@)) ¥f". (2.48) 
Thiveompistes thetureat of chim: 
if AO 0sthen by cls re-muse have that for all 
(Q-(Q))¥@) =0, + Q¥(2) = (Q)¥@). (2.49) 


We see that W is an eigenstate of Q, which indeed means there is no uncertainty in the measurement. 
Of course if W is an eigenstate of Q then again QV = (Q)W and the uncertainty vanishes. All in all, 
we have established the equivalence 


AQv = 0 <> VWisan eigenstate of Q. (2.50) 
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1 Stationary States 


Consider the Schrédinger equation for the wavefunction V(x,t) with the assumption that the potential 
energy V is time independent: 


he oO 


2m Ox 


ih— = HW(x,t) = ( +V¥(o)] U(a,t), (1.1) 


where we displayed the form of the Hamiltonian operator H with the time independent potential 
V(x). Stationary states are a very useful class of solutions of this differential equation. The signature 
property of a stationary state is that the position and the time dependence of the wavefunction 
factorize. Namely, 


V(x,t) = g(t) p(x), (1.2) 


for some functions g and w. For such a separable solution to exist we need the potential to be time 
independent, as we will see below. The solution U(x, t) is time-dependent but it is called stationary 
because of a property of observables. The expectation value of observables with no explicit time 
dependence in arbitrary states has time dependence. On a stationary state they do not have time 
dependence, as we will demonstrate. 

Let us use the ansatz (1.2) for V in the Schrédinger equation. We then find 


(nS) wie) = 910) 40), (1.3) 


because g(t) can be moved across H. We can then divide this equation by W(x, t) = g(t)w(2), giving 


1 dg(t) le 2 3 
h = Aw(a). 1.4 
dt ~ Wey) oe 
The left side is a function of only t, while the right side is a function of only x (a time dependent 
potential would have spoiled this). The only way the two sides can equal each other for all values of t 
and zx is for both sides to be equal to a constant E with units of energy because H has units of energy. 
We therefore get two separate equations. The first reads 


in— = Eg. (1.5) 


This is solved by 
git) =e EY", (1.6) 


and the most general solution is simply a constant times the above right-hand side. From the z- 
dependent side of the equality we get 


Hy(x) = Ey(a). (1.7) 


This equation is an eigenvalue equation for the Hermitian operator H. We showed that the eigenvalues 
of Hermitian operators must be real, thus the constant E must be real. The equation above is called 
the time-independent Schrédinger equation. More explicitly it reads 


2 2. 
(--a < v(«)) (0) = Bv(a), (1.8) 


Note that this equation does not determine the overall normalization of ~. Therefore we can write 
the full solution without loss of generality using the g(t) given above: 


Stationary state: W(a,t) =e" y(r), with EER and Ai=Ev. (1.9) 


Note that not only is w(a) an eigenstate of the Hamiltonian operator H, the full stationary state is 
also an H eigenstate 
HW (az,t) = EV(z,t), (1.10) 


since the time dependent function in WV cancels out. 


We have noted that the energy E must be real. If it was not we would also have trouble normalizing 
the stationary state consistently. The normalization condition for WV, if & is not real, would give 


1 = f ae W* (x, t) U(z, t) = fae et t/he—tht/hay* (7 \ab(z) 
= ce BNn HF de w* (x)y(a) = eBImEHI f a w*(x)q)(a). 


The final expression has a time dependence due to the exponential. On the other hand the normal- 


(1.11) 


ization condition states that this expression must be equal to one. It follows that the exponent must 
be zero, i.e., EF is real. Given this, we also see that the normalization condition yields 


[- dou DoE. (1.12) 


How do we interpret the eigenvalue E? Using (1.10) we see that the expectation value of H on the 
state W is indeed the energy 


(A) = fu U* (x,t) H V(x, t) = [eve )EW(e.2) = Bf dx"(2,)U(o,t = £, (1.13) 


Since the stationary state is an eigenstate of H, the uncertainty AH of the Hamiltonian in a stationary 
state is zero. 


There are two important observations on stationary states: 


(1) The expectation value of any time-independent operator Q on a stationary state W is time- 
independent: 


(Q)w(2,t) = pu U* (x, t)QU(z, t) in et Bt/May* (2) Qe Fe hah(x) 


(1.14) 
-| dx ciBthe ABU My (2) Oub(a ey de ¥*(2)O(2) = (Q)ye), 


since the last expectation value is manifestly time independent. 


(2) The superposition of stationary states with different energies not stationary. This is clear because 
a stationary state requires a factorized solution of the Schrédinger equation: if we add two 
factorized solutions with different energies they will have different time dependence and the 
total state cannot be factorized. We now show that that a time-independent observable Q may 
have a time-dependent expectation values in such a state. Consider a superposition 


U(a2,t) = eye P14/ Pah (x) + ce *F24/ Maho (ax), (1.15) 


where 7, and w»2 are H eigenstates with energies E, and Ep, respectively. Consider a Hermitian 
operator Q. With the system in state (1.15), its expectation value is 


Oe i dx W*(«,t)QU(2,t) 
= / * dx (chetEV/Mapt (x) + cet®2"/ bs (a) (cre 2" "Qui (a) + ce ##""Qa(2)) 


= if dx (levies + |c2|?w3 Qube + clege! 1 Ba) Papt aby + chee CEB) hys3 Ou) 


(1.16) 
We now see the possible time dependence arising from the cross terms. The first two terms are 
simple time-independent expectation values. Using the hermitically of Q in the last term we 
then get 


(Q)v = lerl?(Q)y, + leal?(Q) ua 


ee Be 1.17 
+ ciege(B1— Ey da piQye + echo (1- | dx W1(Qy2)" ce 

The last two terms are complex conjugates of each other and therefore 
(Q)x = lel?(Q)u + leal’(Q)u_ + 2Re jens =m i da vides] : (118) 


We see that this expectation value is time-dependent if Fy 4 Ez and (1, Qw2) is nonzero. The 
full expectation value (Q)w is real, as it must be for any Hermitian operator. 


2 Solving for Energy Eigenstates 
We will now study solutions to the time-independent Schrddinger equation 


Hy(«) = Ey(2). (2.19) 


For a given Hamiltonian HI we are interested in finding the eigenstates ~ and the eigenvalues E, 
which happen to be the corresponding energies. Perhaps the most interesting feature of the above 
equation is that generally the value of EF cannot be arbitrary. Just like finite size matrices have a set 
of eigenvalues, the above, time-independent Schrodinger equation may have a discrete set of possible 
energies. A continuous set of possible energies is also allowed and sometimes important. There are 
indeed many solutions for any given potential. Assuming for convenience that the eigenstates and 
their energies can be counted we write 


wi(z), Ey 
wo(z), Er (2.20) 


Our earlier discussion of Hermitian operators applies here. The energy eigenstates can be organized 
to form a complete set of orthonormal functions: 


[i @eie) = 6, (2.21) 


Consider the time-independent Schrédinger equation written as 


aw _ 2m 
da? Rh 


(E-V(z))y. (2.22) 


The solutions 7(xz) depend on the properties of the potential V(x). It is hard to make general state- 
ments about the wavefunction unless we restrict the types of potentials. We will certainly consider 
continuous potentials. We also consider potentials that are not continuous but are piece-wise contin- 
uous, that is, they have a number of discontinuities. Our potentials can easily fail to be bounded. 
We allow delta functions in one-dimensional potentials but do not consider powers or derivatives of 
delta functions. We allow for potentials that become plus infinity beyond certain points. These points 
represent hard walls. 


We want to understand general properties of ~ and the behavior of 7 at points where the potential 
V(a) may have discontinuities or other singularities. We claim: we must have a continuous 
wavefunction. If 7 is discontinuous then 7’ contains delta-functions and w” in the above left-hand 
side contains derivatives of delta functions. This would require the right-hand side to have derivatives 
of delta functions, and those would have to appear in the potential. Since we have declared that our 
potentials contain no derivatives of delta functions we must indeed have a continuous w. 


Consider now four possibilities concerning the potential: 


(1) V(x) is continuous. In this case the continuity of w(x) and (2.22) imply ~” is also continuous. 
This requires 7’ continuous. 


(2) V(a) has finite discontinuities. In this case y” has finite discontinuities: it includes the product of 
a continuous w against a discontinuous V. But then 7’ must be continuous, with non-continuous 
derivative. 


(3) V(x) contains delta functions. In this case 7” also contains delta functions: it is proportional 
to the product of a continuous w and a delta function in V. Thus y’ has finite discontinuities. 


(4) V(x) contains a hard wall. A potential that is finite immediately to the left of z = a and becomes 
infinite for z > a is said to have a hard wall at x = a. In such a case, the wavefunction will 
vanish for x > a. The slope y will be finite as x — a from the left, and will vanish for x > a. 
Thus 7’ is discontinuous at the wall. 


In the first two cases =’ is continuous, and in the second two it can have a finite discontinuity. In 
conclusion 


Both ~ and y’ are continuous unless the potential has delta functions (2.23) 
or hard walls in which cases 7)’ may have finite discontinuities. 


Let us give an slightly different argument for the continuity of w and ae in the case of a potential 
with a finite discontinuity, such as the step shown in Fig. 1. 


V(x) 


Figure 1: A potential V(x) with a finite discontinuity at 7 = a. 


Integrate both sides of (2.22) a—« to a+ e, and then take « + 0. We find 


a+e d dw Im a+e 
= E- ; 2.24 
[owe 2 (2) = EL ae eve (2.24) 
The left-hand side integrand is a total derivative so we have 
dy dw Om. f7"" 
hati apes Ss eee —E# : 2.2 
El. Bl. 7 EL, BVO BVO (2.25) 
By definition, the discontinuity in the derivative of w at « = a is the limit as e€ > 0 of the left-hand 
side: ap yp a 
Ag {| —) = li — - — : DD 
(=) e90 ( dx a+e dix ) ( ®) 
Back in (2.25) we then have 
9) a+e 
hel tee dx (V(x) — E)W(2) . (2.27) 
dx e>0 A We 


The potential V is discontinuous but not infinite around x = a, nor is w infinite around x = a and, 
of course, EF is assumed finite. As the integral range becomes vanishingly small about x = a the 
integrand remains finite and the integral goes to zero. We thus have 


Ae (=) =0. (2.28) 


There is no discontinuity in This gives us one of our boundary conditions. 
To learn about the continuity of ~ we reconsider the first integral of the differential equation. The 
integration that led to (2.25) now applied to the range from 29 < a to x yields 


dip(x) 2m 
de ae (E-V(x (2.29) 


Note that the integral on the right is a bounded function of 2. We now integrate again from a — € to 
a+e. Since the first term on the right-hand side is a constant we find 


v(a+e) —V(a—€) = 2€ an 


m ate x 
Zl. - au | de] a. (2.30) 


Taking the ¢ — 0 limit, the first term on the right-hand side clearly vanishes and the second term goes 
to zero because ti dx’ (EF — V(2’)) is a bounded function of x. As a result we have 


Aww =0, (2.31) 


showing that the wavefunction is continuous at « =a. This is our second boundary condition. 


3. Free particle on a circle. 


Consider now the problem of a particle confined to a circle of circumference L. The coordinate along 
the circle is called x and we can view the circle as the interval x € [0, L] with the endpoints identified. 
It is perhaps clearer mathematically to think of the circle as the full real line x with the identification 


gvae+lL, (3.1) 


which means that two points whose coordinates are related in this way are to be considered the same 
point. If follows that we have the periodicity condition 


w(e +L) = Va). (3.2) 


From this it follows that not only w is periodic but all of its derivatives are also periodic. 


The particle is assumed to be free and therefore V(x) = 0. The time-independent Schrédinger 


equation is then 
h? dw 
-——— , =E8£ : : 
2m dx? v() ee) 


Before we solve this, let us show that any solution must have E > 0. For this multiply the above 
equation by *(a) and integrate over the circle x € [0,Z). Since w is normalized we get 


L 
we va) ae = = Ef vw =F: (3.4) 


The integrand on the left hand side can be rewritten as 


2 * 
ami Jy ae O"ae) ~ Ge al * 


and the total derivative can be integrated 


2 
zal (Yas) 


lea) 


Since w(x) and its derivatives are periodic, the contributions from x = L and x = 0 cancel out and we 
are left with 


| +E f |Sp ae =e (3.6) 


h2 Ly dap |2 
——— — > : 
p-# ['iefaee, e 


which establishes our claim. We also see that E = 0 requires w constant (and nonzero!). 
Having shown that all solutions must have EF’ > 0 let us go back to the Schrodinger equation, which 
can be rewritten as 


d?y) 2nE 
Bag Rp y. (3.8) 
We can then define k via ae 
ea So. (3.9) 
h 
Since EF > 0, the constant & is real. Note that this definition is very natural, since it makes 
hk? 
E = — 3.10 
a (3.10) 
which means that, as usual, p = hk. Using k? the differential equation becomes the familiar 
ay 
—— = —ky. Si 
a5 =— BY (3.11) 


We could write the general solution in terms of sines and cosines of kx, but let’s use complex expo- 
nentials: 


W(x) ~ et*?, (3.12) 


This solves the differential equation and, moreover, it is a momentum eigenstate. The periodicity 
condition (3.2) requires 


eth(e+h) — etka _,  gthb 4 _» kL =2nn, ne Z. (3.13) 


We see that momentum is quantized because the wavenumber is quantized! The wavenumber has 


discrete possible values 


27n 
kn = Z neZ. (3.14) 


All integers positive and negative are allowed and are in fact necessary because they all correspond to 
different values of the momentum p, = hk,. The solutions to the Schrédinger equation can then be 
indexed by the integer n: 


ia Nnem, (3.15) 
where N is a real normalization constant. Its value is determined from 
L L ; ; 1 
— ie x)da =f NPdx = NPL > N=, 3.16 
[ e@rvateyae = f = (3.16) 


so we have 


bala) = Saeln = eT | (3.17) 


The associated energies are 


2 hek2 = han? n? = Qn hn? 

Im amL? mL? 

There are infinitely many energy eigenstates. We have degenerate states because E,, is just a 

function of |n| and thus the same for n and —n. Indeed w,, and ~_, both have energy E,. The only 
1 


nondegenerate eigenstate is wo = TE which is a constant wavefunction with zero energy. 


Whenever we find degenerate energy eigenstates we must wonder what makes those states differ- 
ent, given that they have the same energy. To answer this one must find an observable that takes 
different values on the states. Happily, in our case we know the answer. Our degenerate states can be 
distinguished by their momentum: w, has momentum Qnn# and w_» has momentum (—27n4). 

Given two degenerate energy eigenstates, any linear combination of these states is an eigenstate 
with the same energy. Indeed if 


(3.18) 


n 


Hy, =Ey1, Hy: = Epp, (3.19) 
then 
H(ay, + b¥2) = aHy, +bHyy = aEy,+bEv2 = E(ay, + by). (3.20) 
We can therefore form two linear combinations of the degenerate eigenstates yw, and w_, to obtain 
another description of the energy eigenstates: 
Un +n ~ cos(knx) , 
Un — Wn ~ sin(kp2) . 
While these are real energy eigenstates, they are not momentum eigenstates. Only our exponentials 
are simultaneous eigenstates of both H and p. 


(3.21) 


The energy eigenstates wy, are automatically orthonormal since they are p eigenstates with no 
degeneracies (and as you recall eigenstates of a hermitian operator with different eigenvalues are 
automatically orthogonal) : 


L 1 L 2ri(m—n)x 
i i (eile de = al ei se: = “Oras (3.22) 
0 L Jo 


They are also complete: we can then construct a general wavefunction as a superposition that is in 
fact a Fourier series. For any U(zx,0) that satisfies the periodicity condition, we can write 


U(x,0) = So an vn(2), (3.23) 
neZ 
where, as you should check, the coefficients a, are determined by the integrals 


L 
a i da Ut (x) U(x, 0). (3.24) 


The initial state U(x,0) is then easily evolved in time: 
(x,t) = San Un(ae (3.25) 
neZ 
Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this 
document. 
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1 The Infinite Square Well 


In our last lecture we examined the quantum wavefunction of a particle moving in a circle. Here we 
introduce another instructive toy model, the infinite square well potential. This forces a particle 
to live on an interval of the real line, the interval conventionally chosen to be x € [0,a]. At the ends 
O and a of the interval there are hard walls that prevent the particle from going to x > a and x < 0. 
The potential is defined as follows and shown in figure 1. 


eye 0, O<aK<a, (1.1) 
Joo z<0,x2%>0 


It is reasonable to assume that the wavefunction must vanish in the region where the potential is 


A 
Z 


sx 


X=O X= 
Figure 1: The infinite square well potential 


infinite. Classically any region where the potential exceeds the energy of the particle is forbidden. Not 
so in quantum mechanics. But even in quantum mechanics a particle can’t be in a region of infinite 
potential. We will be able to justify these claims by studying the more complicated finite square well 
in the limit as the height of the potential goes to infinity. But for the meantime we simply state the 
fact: 

wa) =O) for 2 Dad fora Sia. (12) 


Since the wavefunction must be continuous we must have that it should vanish at 7 = 0 and at 7 =a: 


1. W(x =0) =0. 
2.0 =a) = 0-5 


These are our boundary conditions. You may wonder about the continuity of the first derivative 
w(x). This derivative vanishes outside the interval and continuity would say that 7’ should vanish at 
0 and at a. But this is impossible. A solution of Schrédinger’s equation (a second order differential 
equation) for which both the wavefunction and its derivative vanishes at a point is identically zero! If 
a solution exist we must accept that w’ can have discontinuities at an infinite wall. Therefore we do 
not impose any boundary condition on 7’. The two conditions above will suffice to find a solution. In 
that solution 7’ is discontinuous at the endpoints. 
In the region x € [0,a] the potential vanishes and the Schrédinger equation takes the form 

2 

xt 
and as we did before, one can show that the energy E must be positive (do it!). This allows us to 
define, as usual, a real quantity k such that 


2nE hk? 
a 
= ES: 1.4 
. h2 “a 2m (1.4) 
The differential equation is then 
dy 
—_ = -k’ 1.5 
at = —BY, (1.5) 


and the general solution can be written as 
w(x) = ccoskx+cosinkz, (1.6) 


with constants c; and cg to be determined. For this we use our boundary conditions. 
The condition ~(z = 0) = 0 implies that c, in Eq 1.6 must be zero. The coefficient of sin ka need 
not be, since this function vanishes automatically for « = 0. Therefore the solution so far reads 


w(x) = cgsinka. (1.7) 


Note that if we demanded continuity of 2’ we would have to ask for w~/(z = 0) = 0 and that would 
make cz equal to zero, and thus w identically zero. That is not a solution. There is no particle if 
wy =0. 


At this point we must impose the vanishing of w at © = a. 


esinka = 0 3 ka =nt > kn = (1.8) 
a 
Here n must be an integer and the solution would be 
n(e) = Nsin(*), (1.9) 
a 


with N a normalization constant. Which integers n are acceptable here? Well, n = 0 is not acceptable, 
because it would make the wavefunction zero. Moreover, n and —n give the same wavefunction, up 
to a sign. Since the sign of a wavefunction is irrelevant, it would thus be double counting to include 
both positive and negative n’s. We restrict ourselves to n being positive integers. 


To solve for the coefficient, we utilize the normalization condition; every y,,(a) must be normalized. 


a 
2 
1 =n? sin? ("ar = NP 4-0 + N=4/-. 
0 a a 


Therefore, all in all, our solutions are: 


Qn Qma2 ’ 


D) h2k2 h2 Dane 
tn = yf sin (2), Ey = = a n=1,2,--- 
a a 


(1.10) 


(1.11) 


Each value of n gives a different energy, implying that in the one-dimensional infinite square well there 
are no degeneracies in the energy spectrum! The ground state —the lowest energy state— corresponds 


to n = 1 and has nonzero energy. 


Figure 2: The four lowest energy eigenstates for the infinite square well potential. The n* wavefunction 
solution w, has n — 1 nodes. The solutions are alternately symmetric and antisymmetric about the midpoint 


a= a. 


Figure 2 shows the first four solutions to the 1-d infinite square well, labeled from n = 1 to n = 4. 


We note a few features: 


1. The ground state n = 1 has no nodes. A node is a zero of the wavefunction that is not at the 
ends of the domain of the wavefunction. The zeroes at « = 0 and x = a do not count as nodes. 
Clearly (a) does not vanish anywhere in the interior of [0, a] and therefore it has no nodes. It 


| 
| 


ve a) X=0 X= Vo x 


Figure 3: The infinite square well shifted to the left to make it symmetric about the origin. 


is in fact true that any normalizable ground state of a one-dimensional potential does not have 
nodes. 


2. The first excited state, n = 2 has one node. It is at x = a, the midpoint of the interval. The 
second excited state, n = 3 has two nodes. The pattern in fact continues. The n-th excited state 
will have n nodes. 


3. In the figure the dotted vertical line marks the interval midpoint x = $. We note that the ground 
state is symmetric under reflection about = 5. The first excited state is antisymmetric, indeed 
a 


its node is at « = $. The second excited state is again symmetric. Symmetry and antisymmetry 
2 


alternate forever. 

4. The symmetry just noted is not accidental. It holds, in general for potentials V(x) that are even 
functions of z: V(—x) = V(x). Our potential, does not satisfy this equation, but this could 
have been changed easily and with no consequence. We could shift the well over so that rather 
than having V(x) = 0 from 0 < xz < a, it extends from a5 <a< 5 and then it would be 
symmetric about the origin x = 0 (see figure 3). We will later prove that the bound states of 
a one-dimensional even potential are either even or odd! Here we are just seeing an example of 
such result. 


5. The wavefunctions y,(x) with n = 1,2,... form a complete set that can be used to expand any 
function in the interval x € [0,a] that vanishes at the endpoints. If the function does not vanish 
at the endpoints, the convergence of the expansion is delicate, and physically such wavefunction 
would be problematic as one can verify that the expectation value of the energy is infinite. 


2 The Finite Square Well 


We now examine the finite square well, defined as follows and shown in figure 4. 


es F <a, W>0 
ve)={ Mg, te: SERS, Vea (2.12) 


| 0, for |z|>a. 


Note that the potential energy is zero for |x| > a. The potential energy is negative and equal to 
—Vo in the well, because we defined Vo to be a positive number. The width of the well is 2a. Note 


Figure 4: The finite square well potential 


also that we have placed the bottom of the well differently than in the case of the infinite square well. 
The bottom of the infinite square well was at zero potential energy. If we wanted to obtain the infinite 
square well as a limit of the finite square well we would have to take Vo to infinity, but care is needed 
to compare energies. The ones in the infinite square well are measured with respect to a bottom at 
zero energy. The ones in the finite square well are measure with respect to a bottom at —Vo. 


We will be interested in bound states namely, energy eigenstates that are normalizable. For this 
the energy E of the states must be negative. This is readily understood. If EF > 0, any solutions in 
the region x > a where the potential vanishes would be a plane wave, extending all the way to infinity. 
Such a solution would not be normalizable. The energy F is shown as a dashed line in the figure. We 
have 

“Ve <0; (2.13) 


Note that since E is negative we have EF = —|E|. For a bound state of energy E, the energy E 
measured with respect to the bottom of the potential is 


E = E-(-W)=W-|E|>0. (2.14) 


Those E are the ones that can be compared with the energies of the infinite square well in the limit 
as Vo > oo. 


What are the bound state solutions to the Schrédinger equation with this potential? We have to 
examine how the equation looks in the various regions where the potential is constant and then use 
boundary conditions to match the solutions across the points where the potential is discontinuous. 
We have the equation Let’s examine the regions, where, for simplicity, we define A(x) by 


ay 2m 
En 


(E—V(x))~ =alayy, (2.15) 


where we have defined the factor a(x) that multiplies the wavefunction on the right-hand side of the 
Schrodinger equation. We then consider the two regions 


e region |x| > a: a(x) is a positive constant. The wavefunction in this region constructed with 
real exponentials. 


e region |z| < a: a(x) is a negative constant. The wavefunction in this region is constructed with 
trigonometric functions. 


The potential V(a) for the finite square well is an even function of a: V(—x) = V(a) We can 
therefore use the theorem cited earlier (and proven later!) that for an even potential the bound states 
are either symmetric or antisymmetric. We begin by looking for even solutions, that is, solutions w 
for which (—2x) = 7(z). 


Even solutions. Since the potential is piecewise continuous we must study the differential equation 
in two regions: 
e |x| <a 
aw 2m 2m 
ga ee Oe 


Vo — |E| is a positive constant thus define a real k > 0 by 


(Vo — |E|) (2.16) 


2m 


= sa (Vo — |BI) > 0, ee (2.17) 


It is interesting to note that this equation is not too different from the free-particle equation 


= 2 Indeed, Vo — |E| is the kinetic energy of the particle and thus k has the usual 


interpretation. The differential equation to be solved now reads 


i" = ky, (2.18) 


for which the only possible even solution is 


we) =coske, |al| <a. (2.19) 


We are not including a normalization constant because, at this state we do not aim for normalized 
eigenstates. We will get an eigenstate and while it will not be normalized, it will be normalizable, 
and that’s all that is essential. We are after is the possible energies. Normalized wavefunctions 
would be useful to compute expectation values. 


e |z| >a 
2m 2m|E| 
This time we define a real positive constant « with the relation 
2m|E| 
— a > 0. (2.21) 
The differential equation to be solved now reads 
wl = ew, (2.22) 


and the solutions are exponentials. In fact we need exponentials that decay as x — too, 
otherwise the wavefunction will not be normalizable. This should be physically intuitive, in a 


classically forbidden region the probability to be far away from the well must be vanishingly 
small. For x > a we choose the decaying exponential 


v(x) = Ae“, e>a, (2.23) 


where A is a normalization constant to be determined by the boundary conditions. More gener- 
ally, given that the solution is even, we have 


W(z) = Ae“, Ja] >a. (2.24) 


It is now useful to note that «? and k? satisfy a simple relation. Using their definitions above we 
see that the energy || drops out of their sum and we have 


ke +? = anne (2.25) 
At this point we make progress by introducing unit free constants €,7, and zo as follows: 
= ka>0O, 
= ets (2.26) 
2 = ees 


Clearly € is a proxy for « and 7 is a proxy for k. Both depend on the energy of the bound state. The 

parameter zo, unit-free, just depends on the data associated with the potential (the depth Vo and the 

width 2a) and the mass m of the particle. If you are given a potential, you know the number zo. A 

very deep and/or wide potential has very large zo, while a very shallow and/or narrow potential has 

small zo. As we will see the value of zg tells us how many bound states the square well has. 
Multiplying (2.25) by a? and using our definitions above we get 


ij eS Be (2:27) 


Let us make clear that solving for € is actually like solving for the energy. From Eq. (2.21), we 
can see 


2m|E|a? 2mVoa? | E| |E| 
7 ny oe eee = =, 32 
Sona SS =o = RMT 01 ; (2.28) 
and from this we get 
|E| &\e 


This is a nice equation, the left hand side gives the energy as a fraction of the depth Vo of the well 
and the right-hand side involves € and the constant zo of the potential. The quantity 7 also encodes 
the energy in a slightly different way. From (2.17) we have 


2 2 
2 = a = (Vy -|E)), (2.30) 


and using (2.14) we see that this provides the energy E, measured relative to the bottom of the 


potential 
2 


E=W-|E| = 7 (2.31) 


2ma? ~ 
This formula is convenient to understand how the infinite square energy levels appear in the limit as 
the depth of the finite well goes to infinity. Note that the above answer for the energies is given by 
the unit free number 7 multiplied by the characteristic energy of an infinite well of width a. 


Let us finally complete the construction. We must impose the continuity of the wavefunction and 
the continuity of ~’ at 2 =a. Using the expressions for w for x < a and for x > a these conditions 
give 

wv continuous at 7 =a => cos(ka) = Ae ** 


2.32 
—KAe *, ee) 


wv’ continuous at r =a = = —ksin(ka) 


Dividing the second equation by the first we eliminate the constant A and find a second relation 
between k and «! This is exactly what is needed. The result is 


ktanka =x — katanka =xa > £=ntann. (2.33) 


Our task of finding the bound states is now reduced to finding solutions to the simultaneous equations 


Even solutions: 77+€? = 2, €=ntann, €,n>0. (2.34) 


These equations can be solved numerically to find all solutions that exist for a given fixed value 
of zg. Each solution represents one bound state. We can understand the solution space by plotting 
these two equations in the first quadrant of an (n,€) plane, as shown in figure 5. 

The first equation in (2.34) is a piece of a circle of radius zo. The second equation, € = ntann, 
gives infinitely many curves as 7 grows from zero to infinity. The value of € goes to infinity for 7 
approaches each odd multiple of 7/2. The bound states are represented by the intersections in the 
plot (heavy dots). 

In the figure we see two intersections, which means two bound states. The first intersection takes 
place near 7 = 7/2 and with large € ~ zg. This is the ground state, or the most deeply bound 
bound-state. This can be seen from (2.29). Alternatively, it can be seen from equation (2.31), noting 
that this is the solution with smallest 7. The second solution occurs for 7 near 37/2. As the radius of 
the circle becomes bigger we get more and more intersections; z 9 controls the number of even bound 
states. Finally, note that there is always an even solution, no matter how small zo is, because the arc 
of the circle will always intersect the first curve of the € = 7 tan7 plot. Thus, at least one bound state 
exists however shallow the finite well is. 


Odd solutions. For odd solutions all of our definitions (k,«,z0,7,€) remain the same. The wave- 
function now is of the form 


sinkx, |z|<a 
¥(2) ee ~~ lel Sea 82) 
Matching 7 and w’ at x = a now gives € = —ncot7 (do it!). As a result the relevant simultaneous 
equations are now 
Odd solutions: 7?+€? =22, €=-—ncotn, £,n>0. (2.36) 


y 


Figure 5: Graphical representation of the simultaneous equations (2.34). The intersections of the circle with 
the 7 tan 7 function represent even bound state solutions in the finite square well potential. The deepest bound 
state is the one with lowest 7. 


Figure 6: Graphical representation of (2.36). The intersections of the circle with the curves € = —ncot7 are 
odd bound-state solutions in the finite square-well potential. In the case displayed there is just one bound state. 


In figure 6 the curve € = —n cot 7 does not appear for 1 < 7/2 because € is then negative. For zo < 
there are no odd bound-state solutions, but we still have the even bound state. 


We could have anticipated the quantization of the energy by the following argument. Suppose you 


Figure 7: Sketching eigenstates of a finite square well potential. The energies are Ey < Ey < E3. 


try to calculate energy eigenstates which, as far as solving the Schrédinger equation, are determined 
up to an overall normalization. Suppose you don’t know the energy is quantized and you fix some 
arbitrary fixed energy and calculate. Both in the even and in the odd case, we can set the coefficient of 
the sin ka or cos ka function inside the well equal to one. The coefficient of the decaying exponential 
outside the well was undetermined, we called it A. Therefore we just have one unknown, A. But we 
have two equations, because we impose continuity of = and of w’ at x =a. If we have one unknown 
and two equations, we have no reason to believe there is a solution. Indeed, generally there is none. 
But then, if we think of the energy EF as an unknown, that energy appears at various places in the 
equations (in k and «) and therefore having two unknowns A and F and two equations, we should 
expect a single solution! This is indeed what happened. 

In figure 7 we sketch the energy eigenstates of a square-well potential with three bound states of 
energies EL, < E2 < E3. A few features of the wavefunctions are manifest: they alternate as even, odd, 
and even. They have zero, one, and two nodes, respectively. The second derivative of w is negative 


for |z| < a and positive for |x| > a (it is in fact discontinuous at x = +a). The exponential decay in 
the region |x| > a is fastest for the ground state and slowest for the least bound state. 


Sarah Geller transcribed Zwiebach’s notes to create the first LaTeX version of this document. 
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1 General properties 


You will be proving the following facts in your homework: 


1. Given the Schrédinger equation with potential V () 


—-——54+(V(xz)- E)y =0 (1.1) 


there are no energy eigenstates with E < min, V(x). In other words, the situation indicated in 
Fig.1 cannot occur. 


VX 


__ E<Min (V(X) 


— —_— 
—_ =" 


Figure 1: There are no solutions with energy F less than the minimum of the potential V(z). 


2. For a one-dimensional potential over —co < x < oo, there are no degenerate bound states. Recall 
that a bound state is a normalizable energy eigenstate. Given that, lim; ... ~ = 0. 


We now show that the reality of V(x) allows us to work with real wavefunctions w(z). Even though 
there are complex solutions, we can choose real ones without loss of generality. 


Theorem 1. The energy eigenstates ~(a) can be chosen to be real. 


Proof. Consider our main equation for the complex wavefunction w(x): 


wo” + (E-V(2))p = 0. (1.2) 


Since (y")* = (~*)” and V(z) is real, the complex conjugation of the above equation gives 


(v")" + (EB - V(x)" = 0. (1.3) 


We see that ~*(x) is another solution of the Schrédinger equation with the same energy. The solution 
y*(x) is different from w(x) if there is no constant c such that ~* = cy. In that case y* and w 
represent two degenerate solutions and, by superposition, we can obtain two real degenerate solutions: 


br= 5+), Vim = yes (1.4) 


These are, of course, the real and imaginary parts of w. If W* = cw the real and imaginary parts yield 
the same real solution. In either case we can work with a real solution. 


If we are dealing with bound states of one-dimensional potentials more can be said: it is not that we 
can choose to work with real solutions but rather that any solution is, to begin with, essentially real. 


Corollary 1. Any bound state ~(x) of a one-dimensional potential is real, up to an overall constant 
phase. 

Proof. Recall that in one dimensional potentials there are no degenerate bound states. This means 
that the two real solutions w, and Wj, considered above must be equal up to a constant that can only 


be real: 
Vim(x) =cv,-(2), with cE R (1.5) 


It then follows that = = %, + iWim = (1+ ic)y,. Writing 1+ ic = V1 4+ c? e*? with real 8, shows that 
w is, up to a constant phase (, equal to a real solution. 

Our next result shows that for a potential that is a symmetric function of x, we can work with 
energy eigenstates that are either symmetric or antisymmetric functions of x. 


Theorem 2. If V(—x) = V(x), the energy eigenstates can be chosen to be even or odd under 7 > —2. 


Proof. Again, we begin with our main equation 
yp" + (E-V(z))p = 0. (1.6) 


Recall that primes denote here derivative with respect to the argument, so 7)”(x) means the function 
“second-derivative-of-y)” evaluated at x. Similarly ~”(—x) means the function “second-derivative-of- 
w” evaluated at —x. Thus we can change x for —x with impunity in the above equation getting 


b"(—2) + =7(E-V(a))o(-2) = 0, (1.7) 


where we used that V is even. We now want to make clear that the above equation implies that (—<) 
is another solution of the Schrédinger equation with the same energy. For this let us define a function 
v(x) and take two derivatives of it 


2 
ole) =¥(-2) + Syl) = ¥(-2)-(-), Bee) =v). 08) 
Using the last equation (1.7) becomes 
d? 2m 
qe Pl) + za (BE — V(2)) p(x) = 0. (1.9) 


showing that y(a) = v(—2) provides a degenerate solution to the Schrédinger equation: Equipped 
with the degenerate solutions w(x) and w(—«) we can now form symmetric (s) and antisymmetric (a) 
combinations that are, respectively, even and odd under x > — 2: 


Ys(x) = 5(¥(2) + ¥(-2)), Yale) = 5(¥(r)- ¥(-2)). (1.10) 


These are the solutions claimed to exist. 
Note that the above proof would not work for the case of odd potentials V(—x) = —V (a). We can’t 
say much in this case! 


Again, if we focus on bound states of one-dimensional even potentials the absence of degeneracy has a 
stronger implication: the solutions are automatically even or odd. 


Corollary 2. Any bound state of a one-dimensional even potential is either even or odd. 


Proof: The absence of degeneracy implies that the solutions ~(x) and v(—x) must be the same 
solution. Because of corollary 1, we can choose w(x) to be real and thus we must have 


w(—z) = cv(z), withceR. (1.11) 


Letting x + —z in the above equation we get w(x) = cy(—x) = c*y)(x) from which we learn that 
c? = 1. The only possibilities are c=+1. So (xz) is automatically even or odd under 2 > —2. 


We used the result of this theorem to find the bound states of the finite square well. Since that 
potential is even, we could restrict our work to search for even bound states and odd bound states. 
The potential cannot have bound states that are neither even nor odd! 


2 Bound states in slowly varying potentials 


We will now consider some insights that classical physics gives us about the behavior of energy eigen- 
states. This is sometimes called the semi-classical approximation because classical physics can 
sometime give an approximate description of the quantum physics. 

We begin with an example we already understand. We consider the total energy E of a particle 
that is the sum of a potential energy V and a kinetic energy K. When V is a function of position, 
K must also be a function of position in order for the sum EF to be conserved, as it must be. In our 
first example, shown in Figure 2, the potential V is constant and the energy FE is greater than V. A 
particle in such a potential will have a constant potential energy K and thus a constant momentum 
p=vV2mK. It is a fact that the wave representing the quantum particle has a de Broglie wavelength 
X equal to Planck’s constant h divided by the classical momentum. 

Indeed from the Schrédinger equation 


2m 2nk pe 
wy" = — =a ( —V)~ = 2 w = a (2.1) 
leading to real solutions of the form 
p 20 
a om = een 2.2 
wv cos (F.r) cos |G" | : (2.2) 


where we see that the wavelength of w is the de Broglie wavelength of a particle with momentum p. 
This wavefunction is real and thus is not a momentum eigenstate. It represents a superposition of a 


x 


Figure 2: On the left is an example of a constant potential, V = Vo < E where K is kinetic energy, and 
E = K+ Vo is the total energy. On the right is a sketch of the wavefunction for a particle in this potential. 
Constant V = > constant K, thus momentum is also constant. 


state with momentum p and a state of momentum p. Even in the classical case, the kinetic energy Kk 
only determines p? and not the sign of p. Our interest here is in real wavefunctions 7(a) as appropriate 
for energy eigenstates and what we have seen is that, for a constant potential, the wavelength of ¢(2) 
is the de Broglie wavelength associated to the classical momentum. 


Consider now the situation depicted in figure 3, where we imagine a classical particle moving in a 
linearly increasing potential V(x). This time the kinetic energy K(x) is also position dependent. As 
a result, the momentum of the classical particle p(x) = \/2mK (zx) is also position dependent. The 
insight now is that, if the potential is slowly varying, to a good approximation the wavefunction will 
have a position-dependent de Broglie wavelength A(a) given by 


(2.3) 


K, Pp decreases—> 
 increases—> 


Figure 3: A linearly varying potential, V(x) = ax. The energy FE = K + V(z) is fixed, and as x increases K (x) 
decreases, p(x) decreases, and (x) increases, resulting in a wavefunction of increasing wavelength. 


cos (57) : sin(S) ; (2.4) 


As we see in Figure 3, because of the linearly growing V(x) a particle with total energy FE will have 
decreasing kinetic energy K(x) as x increases. Thus the classical momentum will decrease, and we 
anticipate that the wavefunction will have an increasing local de Broglie wavelength A(x) as x increases. 
This is illustrated to the right of the figure. We will discuss below what we expect to happen to the 
amplitude of the wavefunction. 

The semiclassical approximation — setting the w(a) wavelength equal to the de Broglie wavelength 
of the classical particle — is accurate when the potential changes slowly. By slowly we mean that the 
change in the potential over a distance comparable to the local de-Broglie wavelength is very small 
compared to the potential: 


Ma) *| < |V(a)). (2.5) 


The left side of the inequality is an estimate for the change in V over a distance A(x), and therefore 
this quantity must be very small compared to the potential for the semiclassical approximation to be 
valid. This inequality above is the key condition for the semiclassical approximation. It is in fact the 
condition that allows one to set up the WKB analysis of the Schrédinger equation (8.06!). 


In figure 4 we show an arbitrary potential V(x) and consider the classical motion of a particle of 
total energy FE. For any point x9, we have V(x) + K(x) = E. The maximum kinetic energy occurs 
for the minimum value of the potential. Classically a particle cannot have negative kinetic energy, 
thus the particle cannot be found at points where V(z) is greater than the energy E. In the figure, 
this happens for x > xp and for x < xz, and these regions are called classically forbidden regions. A 
particle of energy EF will oscillate from x, to xz and back. As it moves its velocity changes, it is a 
function v(x) of position. The points xy and xp are called turning points, because at these points 
the particle motion reverses direction. 
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Figure 4: An arbitrary potential showing turning points x, and ag. The region to the right of xr and the 
region to the left of x, are classically forbidden regions. At any point x € [a,,2R] the sum of the potential 
energy V(x) and the kinetic energy K(x) is equal to the total energy EF. 


The oscillating classical particle spends more time in the regions where its velocity is small and 
less time in the regions where its velocity is large. This has quantum mechanical implications. Indeed, 
in the semiclassical approximation we find that the wavefunction amplitude is bigger at the places the 
particle spends more time and smaller at places it spends less time. This can be quantified. Consider 
the probability |(x)|?dx to find the particle within the infinitesimal region dx around x. This is set 
this proportional to the fraction of time the particle spends at dz: 


|ab(x)|?da ~ < : (2.6) 


where dt the time require to traverse dx and T is the half-period of oscillation, the time to go from 
the left turning point to the right turning point. With v(x) the local velocity of the particle, we have 
dx m 1 9 1 
w(x) |?dx ~ = dx > |p(a)|° ~ —~. (2.7) 
WON ae = er ~ Tre) Pe) 
This relation has to be interpreted with some care. Recall that (a) has very small wavelength for 
large p(x). So |w(zx)|? oscillates between zero and some peak value over very short distances along 2. 


On the other hand, the momentum p(x) appearing on the right has no such oscillations. Therefore, 
|2 


by |w(x)|? in the above relation we really mean the average of |2(x)|? over a few oscillations near 2, 
in other words the square of the amplitude of the wave 7)(xz). Writing the amplitude (a positive real 


number, of course) as Amp(¢()), we have 


Amp(¥(x)) ~ == ~ VA(2). (2.8) 


The amplitude of the wave is proportional to the square root of the de Broglie wavelength. Thus in 
figure 5, the momentum of the particle decreases and A(x) increases as x increases. So we expect the 
amplitude of the wave to increase, as sketched to the right. 


YX) increases —> 
Amplitude inceveases—? 


Figure 5: A potential V(x) and a kinetic energy that decreases for increasing x. Then the de Broglie wavelength 
increases with x. Equation 2.8 then implies that the amplitude of W(z) also increases with wz. 


As a simple illustration of the importance of averaging |~|? consider a state of large energy in the 
infinite well potential, as shown in Fig. 6. Within the box, the classical particle “bounces” between the 
walls with constant velocity. As a result the particle spends the same amount of time in every equal 
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Figure 6: Particle in a one dimensional box. The particle bounces with constant velocity. 


size interval dx within the box. The classical probability distribution P,, within the box is uniform: 
1 1 
Pala: = ; dx — Pa(x) = a (2.9) 


since the integral over the box x € [0,a] of Py(x)dx gives one. Now consider an energy eigenstate 


Un(x) with large n 
yin) = [2s (7 ) > |w (x)|? = sin? _ (2.10) 
es a a é et a )- : 


The associated quantum-mechanical probability density Pgm is 


Pon(z) = bale)? = = sin? (=) (2.11) 
a a 
and it is, for large n a rapidly oscillating function, seemingly quite different from the classical prob- 
ability Pq = 1/a. While the classical probability never vanishes, the quantum probability has many 
zeroes! Still, over arbitrary distances larger than the a/n (presumed small since n is very large) the 
average of the quantum probability density is approaches the classical probability density. Recall that 
the average of sin? over any integer number of oscillations is 1/2, so that its average over any interval 
that includes a large number of oscillations (not necessarily integer) is approximately equal to 1/2. 
We then have 
- 


Average, (Pgm(x)) ~ == Pas). (2.12) 


The semiclassical approximation is an explicit realization of what is loosely called the “Correspon- 
dence principle”, the idea that there are limits of quantum states in which their properties can be 
understood by a classical analysis of the system. 


3. Sketching wavefunction behavior in different regions 


Let’s examine the behavior of an energy eigenstate for the case of a general potential V(a). We rewrite 
the time-independent Schrodinger equation by dividing by ~w to get: 


ro) - A e-Vv(o), (3.1) 


which is convenient because there is no wavefunction on the right-hand side. We will consider the 
equation in two regions and at some special points. 


e E—V(zx) <0, classically forbidden region because the total energy is smaller than the potential. 
In this case the right hand side of eqn. (3.1) is positive. This means there are two possibilities: 
both w and w” are positive or both y and w” are negative. These possibilities are shown here 
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Figure 7: In the classically forbidden region either y > 0, wy’ > 0 or ~ < 0, w” <0. In both cases 
the wavefunction is convex towards the x-axis. 


Such possible behaviors are summarized by saying that the wavefunction is convex towards 
the axis. When classically forbidden regions reach to x = —oo or x = oo the behavior is of the 


type shown below: 


Figure 8: Wavefunctions approaching 7 = —oo or x = oo within classically forbidden regions. 


e E—V > 0, classically allowed region: the right hand side of eqn. 3.1 is negative. Thus either 
w >0,v” <0 or one has wy < 0,y~” > 0. Both options are shown below and are summarized by 
saying that the wavefunction is concave towards the axis: 
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Figure 9: w > 0, vw” <0 ory <0, ~” > 0. The wavefunction is concave towards the x-axis. 


A classically allowed region that is large would have the the top and bottom parts of the figure 
above alternating to form a function resembling a sin or cos. 


e V(xo) = E. Turning point: defined as a point xo at which the potential energy V (zo) is equal to 
the total energy E. Turning points separate classically allowed and classically forbidden regions. 
A turning point is an inflection point in the graph of (a), a point where the second derivative 


vanishes, since 


bla = = .0- (a0) = 20; (3.2) 


Not all inflection points are turning points. In fact nodes of the wavefunction are inflection points. 
This is clear from 7)” = 20 (E — V(x))w, since vanishing 7 implies vanishing 7”. 


Note, however, that having vanishing ~ and y’ at any point in the domain of the wavefunction is 
not allowed. This is because the Schrddinger equation is a second order linear differential equation. 
One can quickly show that 7 = y’ at some point xp implies that all higher derivatives of 7 vanish 
at that point. Assuming the wavefunction has a Taylor expansion we conclude that the wavefunction 
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must vanish identically. 
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Figure 10: It is impossible to have both w and w’ at any point xo: the Schrédinger equation would then force 
w(x) to vanish identically. 


We conclude this section by illustrating the quantization of the energy levels for bound states of 
an even potential. As shown in Figure 11, at the top, we have an even potential and we have marked 
four energies Ey < Eg < E3 < E,. Not all correspond to energy eigenstates. We imagine integrating 
the Schrédinger equation from + = oo down to x = 0. As shown in the small figure slightly to the 
right and below the potential, we assume w > 0 as x — oo. Since any bound state is automatically 
even or odd, the picture at x — —oo must either the one with y > 0 (the even extension) or the one 


Figure 11: An even potential (top figure) and the result of integrating the equation from x = oo towards zero 
for various values of the energy. We get a solution when w(x) for x > 0 can be matched with continuous w and 
w’ to an even or odd extension valid for x < 0. There is no solution for E, and £3. The ground state arises for 
Ey and the first excited state arises for E4. 


with w < 0 (the odd extension). To have a solution, the picture from the right must match properly 
at x = 0 with one of the two possible extensions at x < 0. 
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Consider, for example, the integration of the equation for E = E,. Coming in from x = oo the 
solution has an inflection point and becomes concave towards the axis. Beyond this point the first 
derivative uy’ decreases. Still, as shown in the picture ~’ does not vanish at x = 0, so that the even 
extension for x < 0 does not match properly with the x > 0 solution at 2 = 0, because y fails to 
be continuous. Matching with the odd extension is not a possibility because then ~ would fail to be 
continuous at x = 0. Thus EF; is not an allowed energy. 

As we increase the energy to E = Eo, the turning point occurs for larger x and 7’ manages to get 
exactly to zero at x = 0, matching w and w’ from the even extension. This time we have a solution. 
Of course this occurs for FE = E» precisely; a bit more energy or a bit less energy and the derivative 
is not continuous at x = 0. 

As we increase the energy to £3 the turning point moves to even larger x and y’ is negative by the 
time we get to x = 0, while yw is still positive. The solution can’t be matched with the even extension 
because of the discontinuous derivative. 

If we now increase the energy further to some value Ey, by the time we reach x = 0 the derivative 
y’ will be negative and wy will be exactly zero. As you can see in the figure, we can get a good match 
to the odd extension for x < 0. This is the first excited state. It is an odd wavefunction with a single 
node at x = 0. 


4 Shooting Method 


The shooting method is a numerical method to find the form of bound-state solution to the the 
Schrodinger equation 


(E—V(a))p =0. (4.1) 


The method is easily implemented for even potentials V(x) in which case we can search for even and 
for odd bound states. 

Consider first even bound states. We will attempt to integrate the differential equation from x = 0 
to x = oo. To start at x = 0 we need to fix = and yw’ at this point. Since the normalization of the 
bound state is not determined by the equation, we can choose 


ve=0) = 1, (4.2) 


The derivative must also be zero at this point: if the derivative is not zero the even wavefunction will 
have a discontinuous derivative at x = 0 (can you see that?). Therefore we must set 


(x4 =0) = 0. (4.3) 


To integrate the differential equation we now need to choose some energy. Pick some arbitrary value 
Eo. We know already that arbitrary values of the energy do not yield bound states. So what goes 
wrong if we now integrate the differential equation. What happens typically is that you get a solution 
that could not be normalized. As you work with your computer you notice that beyond certain point 
along x the solution diverges, perhaps by having yw — oo, that is ~ going up. 

You must then change the value of the energy until you find some value £ for which the solution 
also diverges, but this time with y — —oo, or w going down. This is a signal that there is some energy 
in the interval between Eo and FE, for which a normalizable solution exists. You then try to narrow 
the interval. If Ey < EF, you can do that by finding values larger than Eo for which the divergence is 
still up and values lower than FE, for which the divergence is still down. As you narrow the interval 
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you will still find these divergences, but they occur for larger and larger values of x. As you do so, 
you are getting a better and better approximation to the bound-state energy (see Fig.12). 
For odd solutions, the procedure is the same, but the boundary conditions at © = 0 are: 


(4.4) 


The first is needed for a continuous odd function. The second is arbitrary, up to normalization. 


If a potential is not even but has a hard wall, this provides a good starting point for integration of 
the Schrodinger equation. At this point the wavefunction is set to zero and the derivative set to one. 
If the potential has two hard walls, we can integrate starting from one wall and then demand that the 
solution reach exactly zero by the time we get to the second wall. 


A practical note: When working with computers to carry out these numerical integrations, it is 
first necessary to “clean out” the units from the differential equation. The first step in this process is 
to replace x for a unit free variable u. The relation between them is of the form + = bu where b is a 
quantity with units of length built from the parameters in the problem. 


Figure 12: Numerical integration of the Schrédinger equation leads to solutions that diverge up (Ww > 
oo) or down (7 — —oo) beyond some value of z. Bound states are found for values of energies at 
which the divergence changes from up to down, or from down to up. 


Sarah Geller transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this 
document. 
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1 The Delta Function Potential 


Consider a particle of mass m moving in a one-dimensional potential. The potential V(z) is is rather 
singular: it vanishes for all x except for x = 0 at which point it has infinite strength. More precisely, 
the potential is delta function localized at x = 0 and is written as 


V(z) =-ad(z), a>O, (1.1) 


Here a is a constant chosen to be positive. Because of the explicit minus sign, the potential is infinitely 
negative at « = 0; the potential is attractive. The potential is shown in Fig. 1, where we represent the 
delta function by an arrow pointing downwards. 


MG) 


| 
| 


Figure 1: A delta function well. 


We want to know if this potential admits bound states. For a bound state the energy & must 
be negative: this ensures that all of « 4 0 is classically forbidden, and the wavefunction will decay 
rapidly allowing a normalized solution. A bit of intuition comes by thinking of the delta function as 
approximated by a finite square well in the limit as the width of the well goes to zero and the depth 
goes to infinity in such a way that the product, representing the ” area” is finite (the delta function 
is a function with unit area, as it is clear from its integral). In Figure 2 we show two finite-well 
representations and sketch the wavefunction. We can see that the middle region provides the curving 


csi ve) ~ 4860) 


Figure 2: The delta function potential viewed as the limit as the finite square well becomes narrower and 
deeper simultaneously. We expect to get a wavefunction with discontinuous derivative. 


of the wavefunction needed to have a smooth derivative. In the limit as the width of the region goes 
to zero we would expect, if there is a bound state, to have discontinuous derivative. 

We can get further insight by considering units. The dimension-full constants in the problem are 
a,m, and h. Since a delta function has units of one over length, the constant a must have units of 
energy times length for the potential to have units of energy. Thus we have, as units 


a 
B= 12 
L ? ( ) 
but, as usual, the units of energy are 
h2 
BE = —,. 1.3 
mL? (1.3) 
From these two equations we find 
2 2 
ma h? 


The units of energy must be carried by the above combination of the constants of the problem. 
Therefore the energy E, of any bound state must be a number times that combination: 


ma? 


where # is a unit free positive number that we aim to determine. It is good to see a appearing in the 
numerator. This means that as the strength of the delta function increases, the depth of the bound 
state also increases, as we would naturally expect! 
Let us now turn to the relevant equations. We want to find an EF < 0 state. The wavefunction is 
constrained by the time independent Schrédinger equation 
h? d?a 
= (EF —V(x))w. 1.6 
a ge = (E-V (a) (1.6) 


For « #0, we have V(x) = 0, so this becomes 


d?w 2nE 2nE 
ao ( = Jv =Kp, 6 = —= > 0. Gaz) 
The solutions to this differential equation are of the form 
OG OR oe Alls (1.8) 


The potential is even: 6(—x) = d(x), so if we have a ground state it must be even and, of course, have 
no nodes. If there is an excited state, it must be odd and thus have a node at + = 0. The only odd 
solution we can build with the above exponentials is sinh «xx. But aw ~ sinh«ax cannot be normalized, 
it blows up at x = +oo. Therefore there cannot be an excited state in the delta function potential. If 
there are bound states there is just one of them! 

Let us use the above solutions to build the the ground state wavefunction. First, we can see that 
for x > 0 we must discard the solution e“”, because it diverges as x — oo. Similarly we must discard 
e *” for « <0. Since the wavefunction must be continuous at x = 0 the solution must be of the form 


Ae* x¢>0 
= ! 1.9 
ue) ie x<O0. ne) 


Is any value of « allowed for this solution? No, we will get another constraint by considering 
the derivative of the wavefunction and learning that, as anticipated, it is discontinuous. Indeed, the 
Schrodinger equation gives us a constraint for this discontinuity. Starting with 


~~ 4V(zr)p = Ey, (1.10) 


we integrate this equation from to €, with 0 < « 1, a range that includes the position of the 
delta function. This gives 


h? (dw 
2m \ dx 


dy) 
dx 


) a dee =B | dewey. (1.11) 


€ =. =€ 


The integral on the left-hand side returns a finite value due to the delta function. In the limit as 
€ — 0 the integral on the right-hand side vanishes because ~(x) is finite for all 2, while the region of 
integration is contracting away. This yields 
he, ($ dw 
] 
dx 


2m 30 \ dx 


) aw(0) = 0. (1.12) 


€ 


We define the discontinuity Ap of ~’ at x = 0 as 


dp\ _ ,. dy 
Ao (S) = ti (5 


) (1.13) 


€ 


We have therefore learned that 


Ao (=) és am (0). (1.14) 


Note that the discontinuity in 7’ at the position of the delta function is proportional to the value of 
the wavefunction at that point. At a node a delta function would have no effect; y’ would also be 
continuous. 

Applying the discontinuity equation to our solution (1.9), we have 


‘ dw dy as —KE —Ke\ _ _ 2ma 
is Niger a a a ee a 
This relation fixes the value of « 
ma 
and therefore the value Ey of the bound state energy 
hh? K? 1 ma? 
. 2m 2 he in) 


As we anticipated with the unit analysis, the answer takes the required form (1.5) and the undetermined 
constant # takes the value 1/2. 


2 The Node Theorem 


Recall the infinite well potential 


0 O 
V(x) = pani (2.18) 
oo elsewhere. 


The bound states take the form 


Vn(x) = i[2sin (=) (2.19) 


=1,2,... (2.20) 


and corresponding energies 


Note that %, has n — 1 nodes (zeroes). (The points x = 0 and x = a are not nodes, but rather 
endpoints. ) 


Figure 3: A smooth potential that goes to infinity as |x| > oo. 


This leads us to the node theorem. Consider a potential V(x) that is continuous and satisfies 
V (a) > oo as |z| — oo (Fig. 3). This potential has a number of bound states (energy eigenstates that 
satisfy w — 0 as |x| > oo), which we index y1, We, w3,.... Recall also that there are no degenerate 
bound states in one dimension. The node theorem states that w, has n — 1 nodes. We will give an 
intuitive, non-rigorous explanation of this phenomenon. 

For this argument we also recall that (xo) = ¢(ao) = 0 implies that ~(a) = 0 for all x. One can- 
not have vanishing derivative at a zero of the wavefunction. That applies to nodes or finite endpoints. 


ee | 
— — = 4 


ath O Q 


Figure 4: The screened potential Vq(x). 


First, we examine the potential and fix the location of = 0 at a minimum. We then define the 
screened potentials V,(x) as follows: 


Vile) = ie ‘ . . (2.21) 


As shown in Fig. 4, the screened potential V,(x) is an infinite well of width 2a whose bottom is the 
taken from V(x). The argument below is based on two plausible assumptions. First: As a — oo the 
bound states of V(x) become the bound states of V(x). Second: As a is increased the wavefunction 
and its derivative are continuously stretched and deformed. 

When a is very small, V,(2) is approximately a very narrow infinite well with a flat bottom — an 
infinite square well. This is because we chose x = 0 to be a minimum and any minimum is locally flat. 
On this infinite square well the node theorem holds. The ground state, for example, will vanish at the 
endpoints and will have no nodes. We will now argue that as the screen is enlarged we can’t generate 
a node. This applies to the ground state, as we explicitly discuss below, and to all other states too. If 
we can’t generate nodes by screen enlargement the node theorem applies to V(z). 

Why is this the case? Consider how we might develop an additional node while stretching the 
screen. To start with, consider the ground state in the top part of Figure 5. There is no node at this 
value of the screen and we have w/(—a) > 0 (left wall) and 7~’(a) < 0 (right wall). Suppose that as 
we increase a we produce a node, shown for the larger screen a’ below. For this to happen the sign of 
y’ at one of the endpoints must change. In the case shown it is the right endpoint that experiences a 
change in the sign of 7’. With the assumption of continuous stretching there would have to be some 
intermediate screen at which w’ = 0 at the right endpoint. But in that case, ~ = 7’ = 0 at this 
endpoint, and then 7(a) = 0 for all x, which is clearly impossible. 


Figure 5: Introducing a single node requires changing the sign of the derivative at the right endpoint: w’(a) < 0 
but ~’(a’) > 0. At some intermediate screen, the value of w’ at the right endpoint must become zero. But this 
is impossible. 


Figure 6: Introducing two nodes by having the wavefunction cross the z-axis in between the two boundaries 
(compare top and bottom). This is not possible as it would require in an intermediate screen (middle) in which 
w = wy’ =0 at some point. 


It is possible to introduce nodes without changing the sign of y’ at either endpoint. In this process, 
shown in Fig. 6, the wavefunction dips and produces two new nodes. This process can’t take place, 
however. Indeed, for some intermediate screen the wavefunction must be tangential to the x axis and 
at this point we will have w = y’ = 0, which is impossible. 

We conclude that we cannot change the number of nodes of any wavefunction as we stretch the 
screen. The n-th excited state of the tiny infinite square well, with n — 1 nodes will turn into the n-th 
excited state of V(x) with n — 1 nodes. In the tiny infinite square well the energy levels are ordered 


in increasing energy by the number of nodes. The same is true at all stages of the stretching screen 
and therefore for true for V(x). Any two consecutive energy levels cannot get permuted because, by 
continuity, this would require a situation where we have a degeneracy, which is not possible. 


3 Harmonic Oscillator 


The classical harmonic oscillator is a rich and interesting dynamical system. It allows us to understand 
many kinds of oscillations in complex systems. The total energy F of a particle of mass m moving in 
one dimension under the action of a restoring force Ff = —kzx, k > 0, is usually written as 


E = 3mv?+ $ka’. (3.22) 


The first term is the kinetic energy and the second term is the potential energy 


V(x) = Ska’. (3.23) 


The potential is quadratic in x. In such a system the particle performs oscillatory motion with angular 


w= oe +> k= mu". (3.24) 
m 


Trading & for w and using the momentum to express the kinetic energy, we can rewrite E’ as follows 


frequency w given by 


2 
Pp 1 2% 32 
f= — 5mMwW x. 3.25 
This is all for the classical harmonic oscillator. 
The quadratic potential is ubiquitous in physics, as it arises to first approximation when we expand 
an arbitrary potential around a minimum. To show this consider an arbitrary potential V(x) with a 


minimum at xo. For x ¥ xp, we can use a Taylor expansion to write 
V(x) = V(ao) + (x — 20)V' (x0) + 5 (a — 20)°V" (x0) + O((a — 20)*) . (3.26) 


Since 29 is a critical point V’(zo) = 0. Dropping the higher order terms, we then have that the 
potential is approximately quadratic 


V(x) ~ V (20) + 5V" (xo) (x _ x0). (3.27) 


This is a good approximation for x close to vp. Since xo is a minimum V"”(zq) > 0 and this is a 
harmonic oscillator centered at xo and with k = V" (ao). The additive constant V(aq) has no effect 
on the dynamics. 


Faced with the question of defining a quantum harmonic oscillator we are inspired by the above 
expression (3.25) for the energy and declare that @ and p will be operators with [%, p] = ih and that 
the Hamiltonian H is given by 


F Pp R eo : 
H= on smu &, [2p] =i: (3.28) 


V(x) = § mw? 2”. (3.29) 


Note that w has units of frequency: [w] = 1/T. We can use this to construct a characteristic energy 
hw. The quantum harmonic oscillator is a rather natural system directly inspired by the classical 


oscillator. 
Our first step is finding the energy eigenstates, the solutions of the time-independent Schrédinger 
equation: 
W? dp(x) 1 22 
Sr ae - 5 mux p(x) = Ey(z). (3.30) 


Here both E and y(a) are unknown. We expect that energy eigenstates only exist for certain quantized 
values of E. 

As a first step we will clean the equation of dimensionful constants. This helps us appreciate better 
the equation at hand. Moreover, it would allow us to easily put the equation on a computer. Each 
term in the equation must have units of energy times units of y, as we can see by looking at the 
right-hand side of the equation. Note that the units of y are not relevant to the consistency of the 
equation, as y appears on each term. We can therefore ignore the units of y. The units of energy 
on the left- hand side are constructed on the first term by a combination of constants and derivatives 
and in the second term by a combination of constants and powers of x. If we could work with a unit 
free coordinate u instead of x the units of energy would have to be produced just by the constants in 
the problem, and as we have seen, the only possibility is hw. A common hw factor will then simplify 
tremendously the structure of the equation as it will allow us to define a unit free energy! 

We therefore begin by introducing a unit-free coordinate u to replace the conventional coordinate x. 
We set 

oS atl, u unit free, [a] =L, (3.31) 
where a must be a constant with units of length. To determine a in terms of h,m, and w we equate a 


characteristic kinetic energy to a characteristic potential energy: 


h? h 
—, = in a a a (3.32) 
ma mw 


Now, plugging x = au into the time-independent Schrédinger equation yields 


he @ptu) 1 294 
Sma? du2 ~ gmur aru p(u) = Ey(u). (3.33) 


Here, we have used 
d dud 1d 


= = : 3.34 
dx dxdu adu ( ) 
Note that a = hw and mw?a? = hw, so we have 
2 
lis eu) t Mw ou) = Ey(u). (3.35) 
u 


We can see that things are working. As expected the energy units on the left-hand side are carried by 
hw! Multiplying by a we reach 


d?p(u) 
du? 


where we have defined a unit free energy €: 


t+u%y(u) = Ey(u), (3.36) 


£ Ea snwe. (3.37) 


tT 
ES 


If we know the pure number € we then know the energy &. Rearranging, we reach the cleaned-up 
version of the time-independent Schrodinger equation: 


dy 


fo = Ww _e)p. (3.38) 


This is our simplified, unit free version of the time-independent Schrddinger equation. It is clearly less 
cluttered than (3.30). 

The above differential equation must have solutions for all values of the energy parameter €, 
after all you could integrate it on a computer! Quantization must arise because solutions are not 
normalizable except for special values of €. To understand this issue as it relates to the equation we 
examine solutions for large values of |u|. In this limit, € can be ignored as compared to u?, and we 
have the approximate equation 

g"(u) © wp(u). (3.39) 
This equation cannot be solved by any polynomial! If y is a polynomial of degree n, the degree of the 
left hand side would be n — 2 and that of the right-hand side n + 2. This cannot work. Let’s try a 


solution of the form 
7. 
y(u) = uke /?, (3.40) 


The leading term in y” comes when we differentiate the exponential: 


y"(u) & a?uy(u) as |u| > oo. (3.41) 
Comparing with (3.39) we have solutions for a = +1, in which case we have 
plu) & Auke““/? 4+ Buke/? as |u| + 00. (3.42) 


The solution with coefficient B would not yield an energy eigenstate because it diverges as |u| — co 
and would not be normalizable. Note that the u* factor played no role in the analysis. This factor, 
-u’/2 could be a solution of the differential equation. 
This analysis suggests that, for our purposes, we should write 


however suggests that a polynomial multiplying e* 


p(u) = h(uje™/? (3.43) 


Note that there is no assumption or loss of generality in writing this expression. Indeed, any function 


~w*/2 as it is immediately clear (y(uje”/2)e—W/2, 


y(u) and be written some other function times e 
In writing (3.43) we are only hoping that the differential equation for h(w) is simpler. Clearly, if we 
find h(u) we have found y(u). We actually expect that h(u) may be a polynomial because the ansatz 
captures the large |u| dependence that prevents the solution for y(u) from being a polynomial. 
Plugging (3.43) into (3.38) and simplifying, we find a second-order linear differential equation for 


h(u): 


@h . dh 
2 lh=0O. 44 
sr US + (E-1)h=0 (3.44) 


It is actually possible at this point to see that getting a polynomial solution requires quantization of € 
Indeed, assume that h(w) is a polynomial of degree j: 


A(u) = wi tay t+agui7 +... (3.45) 


In the above equation the first term is then a polynomial of degree j — 2. Each of the other two terms 
are polynomials of degree 7. For the equation to make sense the coefficient of contributions to the 
coefficient of uJ and u/—! must vanish. The coefficient of u/ is 


Coefficient of uJ: -27+E€-1=0 > E€=2j4+1. (3.46) 


Thus, we get the quantization of energy: a polynomial solution h(w) of degree j requires € = 27 + 1. 
You may wonder about the subleading term of degree 7 — 1 whose coefficient must also vanish. 


Coefficient of uJ~': (—-2(j —1)+€—-1)a; =0. (3.47) 


Since the energy € has already been fixed, the only way to satisfy this condition is to set ay = 0. Thus 
the polynomial is actually of the form 


A(u) = wtagud2 +... (3.48) 


If this is supposed to lead to an energy eigenstate the vanishing of a; could have been anticipated. 
Since the harmonic oscillator potential is even we know that bound states must be either even or odd. 
~u*/2 is even, the solution y(u) will be either even or odd if h(w) is even or odd. If a1 had not 
vanished, h(u) would have two consecutive powers of u and could not be either even or odd. 


Since e 
We can analyze the equation more systematically using a series expansion: 
(oe) 
Ala) = De anu. (3.49) 
k=0 


A simple way to plug into the differential equation (3.44) is to select from each term the contribution 
to the coefficient of u?. For this one can imagine the terms aju? + aj41u!t! + aj42u!t? in h(u) and 
select the piece that contributes to the coefficient of u/: 


2 
Contribution from: _ > (9 +2)(9 + Iaj+e 
U 
Contribution from: — 2 S = 270; (3.50) 
u 


Contribution from: (E-—1)h : (E—1)a; 


The total coefficient of uJ in the left hand side of the differential equation must be set to zero, for all 
values of 7, for the differential equation to be satisfied. Therefore 


This can be written as a recursion relation: 


ae 


aj42 = G+DG+D aj, (3.52) 


This is a two-step recurrence relation. If you choose some ag you can construct a solution that contains 
only odd coefficients, az, a4,... as determined recursively by the above relation. That solution, of the 
form 

ay + agu? + agu*+--- (3.53) 
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would be even. Another solution is constructed by choosing some a, and then using the above recursion 
to find a3,a5,.... That solution, of the form 


ayu+ agu® Beets (3.54) 


would be odd. For any arbitrary value of € both solutions of the differential equation (3.44) exist, but 
neither one would be polynomial and neither one would be expected to be a good energy eigenstate. 
The general solution to (3.44) with arbitrary € is thus determined by the two constants (ag, a1) as 
they together determine all coefficients. This makes sense, because ag = h(0) and a, = h’(0) and 
the solution of a second order differential equation is be determined by knowing the function and its 
derivative at any point. 

Let us now demonstrate that if the series for h(w) never stops the corresponding y(u) is not an 
acceptable energy eigenstates. Let us see what would be the large u behavior of h(u) if it does not 
terminate. For large 7 the recursion relation (3.52) gives 


OES (3.55) 
aj j 
What kind of function grows this way? Note that 
ge eee 1 

9b st 2) 0 j 

eS S> ai a) = pe Gjai . (3.56) 
n=0 je even 
This series has coefficients cj = Ga for even 7, and so we see that 


cy (GH+2)/2)! f+27° 9’ 


for large j. This is just the behavior noted in (3.55) for h(w). So, if the series for h(u) does not 
terminate the wavefunction is behaves like 


p(u) = h(uje” /? a Mae Ear aon (3.58) 


which is the bad solution we identified in (3.42). This proves that h(u) must be a polynomial and the 
recursion relation must terminate for us to get an energy eigenstate! 

Now we discuss how to get a polynomial h(w), although the main conclusion was anticipated earlier 
in (3.46). If h(u) is to be of degree j it must have non vanishing a; and vanishing aj;+2, as determined 
from the recursion relation (3.52). The numerator in this recursion relation must vanish and we must 
choose € such that 

2j+1-E=0. (3.59) 


The solution will then take the form: 
h(u) = ajul At aj;_2w? stars (3.60) 


with powers decreasing in steps of two because this is what the recursion relation demands for having 
a solution. The solution will therefore be automatically even (if 7 is even) or odd (if 7 is odd). Say j 
is even and the solution is even with energy 27 + 1 as required. The second solution of the differential 
equation for that value of the energy would be odd, but the energy 27 +1 that made the even solution 
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terminate will not make the odd solution terminate. This means that the second solution of the 
differential equation is not an energy eigenstate. 
We usually call the degree j of the solution using the letter n. Then, 


E=2n+1, n=0,1,2,... (3.61) 


corresponds to the polynomial solution 


Aly =e ag = eee. es OO, (3.62) 
The energy of the solution yp(u) = hn (ue /? is 
h h 
pS so = eae: (3.63) 
2 2 
We have ; 
Ey, = hw(n+5), ee eo eae (3.64) 


We see that the energies are quantized and the energy levels are evenly spaced. The ground state 
energy is Ey = fw/2. The corresponding power series solutions h,,(u) are the Hermite polynomials, 
usually denoted as H,,(w) 

S000) eee ye creer (3.65) 


The factor of 2” here is a convention choice. The Hermite polynomials are solutions of (3.44) with 


€ = 2n+ 1 therefore they satisfy the differential equation 


a Hy, dH, 
Tab 2u aA + 2nH, = 0. (3.66) 
The first several Hermite polynomials are 
Ao(u) = 1 
Ay(u) = 2u 
u(t) ; (3.67) 
Ho(u) = 4u* — 2 
H3(u) = 8u? — 12u 


The generating function for the Hermite polynomials is an exponential, with formal parameter z: 


oe) n 


ent 42e4 — H,,(u). (3.68) 


u=0 
It is not too hard to show that the polynomials defined by this expansion satisfy the requisite differential 
equation (3.66) and are normalized as claimed in (3.65). 


: : 2 0129 2_ oh 
Let us write the energy eigenstates in terms of x. Recalling that u“ = 2*/a*, where a“ = —- the 
relation 
2 
gn(u) ~ Hn(ule™”?, (3.69) 


then gives us 
vn(z) = Nn Ha (« me) eH” =n =0,1,2,..., (3.70) 
where N,, is a normalization constant. 


Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this 
document. 
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1 Algebraic Solution of the Oscillator 


We have already seen how to calculate energy eigenstates for the simple harmonic oscillator by solving 
a second-order differential equation, the time-independent Schrédinger equation. 


Let us now try to factorize the harmonic oscillator Hamiltonian. By this we mean, roughly, writing 
the Hamiltonian as the product of an operator times its Hermitian conjugate. As a first step we rewrite 
the Hamiltonian as 63 

a= smu” (« + ss) . (1.1) 
Motivated by the identity a? + 6? = (a — ib)(a + 7b), holding for numbers a and b, we examine if the 
expression in parenthesis can be written as a product 


nA nA a2 . 
ae ets . caper 
(s-=) (s+) a + P 5 + — (ép - £2), 
mw mw mw? mw 


(1.2) 


where the extra terms arise because < and , as opposed to numbers, do not commute. We now define 
the right-most factor in the above product to be V: 


ve“a4—2, (1.3) 
Mw 


Since @ and p are Hermitian operators, we then have 
i 
yi oe (1.4) 
and this is the left-most factor in the product! We can therefore rewrite (1.2) as 


ad h 
44 = viv+—1, (1.5) 
MAU) mw 


and therefore back in the Hamiltonian (1.1) we find, 


H = smu VIV + $hw1. (1.6) 


This is a factorized form of the Hamiltonian: up to an additive constant Ep, H is the product of a 
positive constant times the operator product V'V. We note that the commutator of V and V? is 
simple 
2 ip. ip Geo 98 are 2h 
m mw mw mw 


ev. ev = ie (1.8) 


This suggests the definition of unit-free operator operators @ and @!: 


This implies that 


é = V 
an. * 
(1.9) 
ai eA 
= Fs Vl 
Due to the scaling we have 
[a@,a7] =1. (1.10) 


The operator @ is called annihilation operator and 4! is called a creation operator. The justification 
for these names will be seen below. From the above definitions we read the relations between (4, a!) 


and (4 ,p): 


(1.11) 


(1.12) 


While neither G nor a! is hermitian (they are hermitian conjugates of each other), the above equations 
are consistent with the hermiticity of @ and ». We can now write the Hamiltonian in terms of the a 
and a operators. Using (1.9) we have 


2h 
ViV = —ala, (1.13) 
Mw 


and therefore back in (1.6) we get 


(1.14) 


The above form of the Hamiltonian is factorized: up to an additive constant H is the product of a 
positive constant times the operator product @'@. In here we have dropped the identity operator, which 


is usually understood. We have also introduced the number operator N. This is, by construction, a 
hermitian operator and it is, up to a scale and an additive constant, equal to the Hamiltonian. An 
eigenstate of H is also an eigenstate of N and it follows from the above relation that the respective 
eigenvalues E' and N are related by 


EB = htw(N +4). (1.15) 


Let us now show the powerful conclusions that arise from the factorized Hamiltonian. On any 
state w that is normalized we have 


(H)y = (Hy) = hub, ala) + phw(},), (1.16) 
and moving the @! to the first input, we get 
(H)y = hw(av,ap) + thw > Zhw. (1.17) 


The inequality follows because any expression of the form (y, y) is greater than or equal to zero. This 
shows that for any energy eigenstate with energy FE: Hw = Ew we have 


Energy eigenstates: E > shw . (1.18) 


This important result about the spectrum followed directly from the factorization of the Hamiltonian. 
But we also get the information required to find the ground state wavefunction. The minimum energy 
shi will be realized for a state w if the term (aw ; aw) in (1.17) vanishes. For this to vanish ay must 
vanish. Therefore, the ground state wavefunction yo must satisfy 


&yo = 0. (1.19) 


The operator @ annihilates the ground state and this why @ is called the annihilation operator. Using 
the definition of @ in (1.11) and the position space representation of p, this becomes 


(s+ a4) yo(z) =0 A (+45) yo(z) =0. (1.20) 
mw i dx 

Remarkably, this is a first order differential equation for the ground state. Not a second order 
equation, like the Schrodinger equation that determines the general energy eigenstates. This is a 
dramatic simplification afforded by the factorization of the Hamiltonian into a product of first-order 
differential operators. The above equation is rearranged as 


dyo mw 
— SS SS ; 1.21 
dx ae ( ) 


Solving this differential equation yields 


(1.22) 


where we included a normalization constant to guarantee that (Yo, 0) = 1. Note that yo is indeed 
an energy eigenstate with energy Ep: 


Heo =hw(alat+4)yo =thwy > Eo = 


Ne 


hw. (1.23) 


Before proceeding with the analysis of excited states, let us view the properties of factorization 
more generally. Factorizing a Hamiltonian means finding an operator A such that we can rewrite the 
Hamiltonian as AtA up to an additive constant. Here At is the Hermitian conjugate of A, an operator 
that is defined by 

(,Aty) = (Av,¢). (1.24) 


We say that we have factorized a Hamiltonian H if we can find a A for which 


(1.25) 


where £o is a constant with units of energy that multiplies the identity operator. This constant 
does not complicate our task of finding the eigenstates of the Hamiltonian, nor their energies: any 
eigenfunction of At A is an eigenfunction of H. Two key properties follow from the factorization (1.25). 


1. Any energy eigenstate must have energy greater than or equal to Eo. First note that for an 
arbitrary normalized (x) we have 


(b, Hb) = (p,ATAY) + E0(v,¥) = (Av, Av) + Eo, (1.26) 


Since the overlap (Ar, Aw) is greater than or equal to zero, we have shown that 
(b,Hp) > Bo. (1.27) 


If we take y to be an energy eigenstate of energy E: H w = Ew, the above relation gives 
E> Ep. (1.28) 
This shows, as claimed, that all possible energies are greater than or equal to Ep. 


2. A wavefunction wo that satisfies 
Avo = 0, (1.29) 


is an energy eigenstate that saturates the inequality (1.28). Indeed, 
Hyp = AtAyo+ Eovo = Al(Ado) + Env = Eno. (1.30) 
The state wo satisfying Aw) = 0 is the ground state. For conventional Hamiltonians this is a 


first order differential equation for wp) and much easier to solve than the Schrédinger equation. 


2 Operator manipulation and the spectrum 


2 


We have seen that all energy eigenstates are eigenstates of the Hermitian number operator N = at 
This is because H = hw(N + 3). Note that since 4yo = 0 we also have 


Ngo = 0. (2.1) 


We can quickly check that 


[N,4] = [ata,a] = [al,aja = -a, ax 
[Na] = [ala,a'] = aaa] = al, | 
which we summarize as 
[N,a] = -4, 
(2.3) 
[N, alc al, 
Using these identities and induction you should be able to show that: 
N,(4)*] = —k(a)*, 
L | (a) (2.4) 


[Ny te") = "ka", 


These relations suggest why N is called the number operator. Acting on powers of creation or anni- 
hilation operators by commutation it gives the same object multiplied by (plus or minus) the number 
of creation or annihilation operators, k in the above. Closely related commutators are also useful: 


yall =a (arr? 


(2.5) 
[@, (atyF] = k(aty 


5)*]. We will also make use of the following 


These commutators are analogous to [p,(#)*] and [#, (p) 
Lemma which helps in evaluations where we have an operator A that kills a state ~ and we aim to 


simplify the action of AB , where B is another operator, acting on ~. Here is the result 


(2.6) 
AB = [A,B]+BA, (2.7) 

as can be quickly checked expanding the right-hand side. It then follows that 
ABy = ((A,B])+BA)y = [A, Bly, (2.8) 


because BA w = B(Av) = 0. This is what we wanted to show. This is all we need to know about 
commutators and we can now proceed to construct the states of the harmonic oscillator. 


Since @ annihilates yo consider acting on the ground state with a!. It is clear that a’ cannot also 
annihilate yo. If that would happen acting with both sides of the commutator identity [4 ,at] = 1 on 
yo would lead to a contradiction: the left-hand side would vanish but the right-hand side would not. 
Thus consider the wavefunction 

gy = ayo. (2.9) 


We are going to show that this is an energy eigenstate. For this purpose we act on it with the number 
operator: 
Neg, = Nalyo = [N,a!]p0, (2.10) 


where we noted that Nyo = 0 and used Lemma (2.6). Given that [N, at] = at, we get 
Ny = al yo = Yi. (2.11) 


Thus y is an eigenstate of the operator N with eigenvalue N = 1. Since yo has N eigenvalue zero, 
the effect of acting on yo with a! was to increase the eigenvalue of the number operator by one unit. 
The operator @! is called the creation operator because it creates a state out of the ground state. 
Alternatively, it is called the raising operator, because it raises (by one unit) the eigenvalue of N. 
Since N = 1 for vy, it follows that vy, is an energy eigenstate with energy E, given by 


Ey, =hw(1+4) = hw. (2.12) 
It also turns out that vy, is properly normalized: 


t 


(91, 91) = (@' yo, a' yo) = (Yo, 44" yp) , (2.13) 


where we used the Hermitian conjugation property to move the 4! acting on the left input into the 


right input, where it goes as (a')' = @. We then have 


(¥1, 91) = (Yo, 44" p0) = (Yo, [a, "] yo) = (Yo, Yo) = 1, (2.14) 


where we used (2.6) in the evaluation of Ga‘). Indeed the state y, is correctly normalized. 


Next consider the state 
vy = aalyo. (2.15) 


This has 


Ny, =Nalatyo = N, ata yo = 2atatyo = 26%, (2.16) 


SO (2 is a state with number N = 2 and energy Ey = hw. Is it properly normalized? We find 


(2,92) = (aa! 0, ala! yo) = (yo, aG4'A' G0) = (v0, 4[4, @'"] yo) ea 
= (Yo, 244" yp) = (yo, 0)' = 2. 
The properly normalized wavefunction is therefore 
a eee (2.18) 
J2 
We now claim that the n-th excited state of the simple harmonic oscillator is 
(2.19) 
Exercise: Verify that this state has N eigenvalue n. 
Exercise: Verify that the state y,, is properly normalized. 
Since the N eigenvalue of y is n, its energy Ey, is given by 
E, = hw(n+4). (2.20) 


Since the various states yp are eigenstates of a Hermitian operator (the Hamiltonian H) with different 
eigenvalues, they are orthonormal 


(Yr: Ym) = Smn- (2.21) 


We now note that ay, is a state with n—1 operators a! acting on yo because the G eliminates one 
of the creation operators in yy. Thus we expect G~n ~ Yn_1. We can make this precise 


1 n =I 
_ at n _ A at n _ at n 
a = a— (a = —=(|4, (a = —=(4 : 2.22 
At this point we use (2.19) with n set equal to n — 1 and thus get 
apn = Val (n-D!lyn-1 = Vnyp-1- (2.23) 


By the action of a! on y, we get 


, 1 
(aly ry = Vint Dents = Vn FT ents. (2.24) 


(2.25) 


These relations make it clear that @ lowers the number of any energy eigenstate by one unit, except 
for the vacuum yo which it kills. The raising operator @' increases the number of any eigenstate by 
one unit. 


Exercise: Calculate the uncertainty Az of position in the n-th energy eigenstate. 


Solution: By definition, 
(Az), = (8 )on — (25, - (2.26) 
The expectation value (i) vanishes for any energy eigenstate since we are integrating x, which is odd, 


against |y~,(x)|?, which is always even. Still, it is instructive to see how this happens explicitly: 


(2) on = (Yn Bn) — —— (Yn »(@+ a!) pn) ) (2.27) 
using the formula for @ in terms of @ and @!. The above overlap vanishes because G~, ~ Yn—1 and 
al yn ~ Yn41 and both yp_1 and yy41 are orthogonal to y,,. Now we compute the expectation value 
of £7 


h 


(Pon = (Pn s#en) = so (ens (@+@)(4 + A!) pn) 
Mad 
i (2.28) 


Since G4y~n ~ Yn—2 and alalon ~ Ynt+2 and both yyp_2 and Yn+2 are orthogonal to yp, the aa and 
atat terms do not contribute. We are left with 
A 


(Pon = ga (¥n, (Gal + aa) pn). (2.29) 


At this point we recognize that ata = N and that aat = [a,at]+a1@ =1+N. Asa result 


= so (en, (1+2N)y,) = soll +2n). (2.30) 


We therefore have ; 
(Az)? = (n+4). (2.31) 


mw 
The position uncertainty grows linearly with the number. 
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1 The Step Potential 


X*0 
Figure 1: The step potential. 


We now begin our detailed study of scattering states. These are un-normalizable energy eigenstates. 
They simply cannot be normalized, just like momentum eigenstates. These energy eigenstates are 
not states of particles, one must superpose scattering states to produce normalizable states that can 
represent a particle undergoing scattering in some potential. Here we examine the step potential 


(Figure 1), defined by 
0, «— U; 
V(x) = - (1.1) 
Yo, 220. 
Our solutions to the Schrédinger equation with this potential will be scattering states of definite 


energy EF. We can consider two cases: E > Vo and E < Vo. In both cases the wavefunction extends 
infinitely to the left and is non-normalizable. Let us begin with the case E > Vo. 


2 Step Potential with EF > VY 


ers e w=” 


Figure 2: The energy F of the stationary state is greater than the step Vo. The full « axis is classically allowed. 


The stationary state with energy E is of the form 
VaD=v@je”, (2.2) 


and we will focus first on the unknown w(x). In order to write a proper ansatz for ~(a) we visualize a 
physical process in which we have a wave incident on the step barrier from the left. Given such a wave 
traveling in the direction of increasing x, we would expect a reflected wave and a transmitted wave. 
The reflected wave, moving in the direction of decreasing x, would exist for x < 0. The transmitted 
wave, moving in the direction of increasing x, would exist for « > 0. The ansatz for the energy 
eigenstate must therefore contain all three pieces: 


Ae? + Be“*® <0, 
v(x) = (2.3) 


Ceike x>0. 


Recall that e’**, with k > 0, represents a wave moving in the direction of increasing x, given the 
universal time dependence above. Therefore A is the coefficient of the incident wave, B is the coefficient 
of the reflected wave, and C is the coefficient of the transmitted wave. The waves for x < 0 have 


wavenumber & and the wave for x > 0 has wavenumber k. These wavenumbers are fixed by the 


Schrodinger equation 
2nE = 2m(E — Vo) 
Qe. DS os 

—— WD? k = ee (2.4) 
There are two equations that constrain our coefficients A,B, and C: both the wavefunction and 
its derivative must be continuous at 7 = 0. With these two conditions we can solve for B and C in 
terms of A. This is all we could expect to do: because of linearity the overall scale of these three 
coefficients must remain undetermined. In fact, we can think of A as the input value and B and C as 


output values. Let us begin: 


e v(x) must be continuous at x = 0. Thus 
A+B=C. (2.5) 


e w(x) must be continuous at 2 = 0. Thus 


ikA-—ikB =ikC > A-Be= 


| sl 
Q 
— 
NO 
a 
we 


Solving for B and C in terms of A, we get 


Bo k- 
A k+ 


a] a1 
Q 
to 
>= 


= _ 2; 


If A is real B and C are real. For E = Vo, we have k = 0 and equations (2.7) give B = A and C = 2A. 
Therefore, for E = Vo the energy eigenstate is 


B=y2 t= ie 3 (2.8) 


and looks like this: 


Figure 3: Energy eigenstate for E = Vo. 


We get further insight into the solution by evaluating the probability current to the left and to the 
right of the x = 0 step. Recall the form of the probability current for a wavefunction ~ is 


h Ow 
J=—I + 2.9 
un eae (2.9) 
A short calculation shows that the current Jz to the left of the step is 
hk hk hk 
J, = —(|A? -|B)?) =Ja-—Je, Ja=—|A?, Jp =—|BP. (2.10) 
m m m 


There is no interference arising from the incident and reflected waves. The total current to the left of 
the step is simply the current J,4 associated with the incident wave minus the current Jg associated 
with the reflected wave. The current Jr to the right of the step is 


_ _ hk 
jn=— |CP = Jo. (2.11) 


In any stationary solution there cannot be accumulation of probability at any region of space because 
the probability density p is manifestly time-independent. While probability is continuously flowing 
in scattering solutions, it must be conserved. From the conservation equation gl + 8p = 0, the time 
independence of p implies that the current J must be x-independent. In particular, our solution (2.7) 
must imply that J; = Jr. Let us verify this: 


hk hk bk 
— —(\A\? — |B?) = 1 —} }|Al? 
dt | |° — |BI?) =( (=) ) | 


_ “*( Akk ) 2 hk 4k? 


(2.12) 


= = — ——__|A/? = AkIC|? = Jr, 
(b+ bP 5|Al IC| R 


om (k+R) 


Ic/? 


m 


as expected. The equality of Jz and Jr implies that 

Jp Jeo 
=—+—. 2.13 
JA 7 JA ( ) 


We now define the reflection coefficient R to be the ratio of the probability flux in the reflected wave 


Ja—JSJp=Jo 7 JnzJptdo - 1 


to the probability flux in the incoming wave: 


Jpg |BI? ( . 
R= = = a els 2.14 
Ja |Al? k+k} — i) 


This ratio happens to be the norm squared of the ratio B/A, and it is manifestly less than one, as it 
should be. We also define the transmission coefficient T to be the ratio of the probability flux in the 
transmitted wave to the probability flux in the incoming wave: 


a = ; z 
pa de _kiecP _ kb 4k Akh (2.18) 
Ja k\Al? k(k+k)? (k+k)? 


The above definitions are sensible because R and T’, given in terms of current ratios, add up to one: 
R+T =1, (2.16) 


as follows by inspection of (2.13). Note that T # Ce because of wavenumbers to the right and to the 
left of the step are not the same. 

Recall that for E = Vo we found k = 0. In that case we have full reflection: R = 1 and T = 0. 
Indeed the probability current associated with the constant wavefunction that exists for x > 0 (see 
(2.8)) is zero. Additionally we can give an argument from continuity. The coefficients R and T must 
be continuous functions of the energy E. For E < Vo we expect T = 0 since the forbidden region is all 
of « > 0 and an exponentially decaying wavefuntion cannot carry probability flux. If T= 0 for any 
E < Vo it must still be zero for EF = Vo, by continuity. 


3 Step Potential with EF < VY 


Figure 4: The step potential barrier. 


When E < Vo the region x > 0 is a classically forbidden region. Let us try to solve for the energy 
eigenstate without re-doing all the work involved in solving for B and C' in terms of A. For this 
purpose we first note that the ansatz (2.3) for x < 0 can be left unchanged. On the other hand, for 


x > 0 the earlier solution 
42 = 2m(E — Vo) 


V(x) = Ge, he ) 


(3.17) 


should become a decaying exponential 


_ 2n(Vo — E 
en amy Os ae = amnVo" 2) : (3.18) 
We note that the former becomes the latter upon the replacement 


ko ik. (3.19) 


This means that we can simply perform this replacement in our earlier expressions for B/A and C/A 
and we obtain the new expressions. In particular from (2.7) we get 


Bk-it 
cies 3.20 
A k+ikt oe) 

Therefore B i(k—in) : 

Uk — 1K Kt 2i6(E) 
= = = 3.21 
A i(k +ik) Kk — ik ‘i : ee} 
with 


6(E) = tan~*(=) = tan" (/ — =) : (3.22) 


Since the magnitude of A is equal to the magnitude of B, we have J4 = Jp and Jo = 0. Thus T = 0 
and R = 1. As noted above, the ratio B/A is a pure phase. The phase of the numerator « + 7k is 6(E£) 
and the phase of the denominator « — ik is —é(E), thus the total phase 26(£) for the ratio. We did 
not absorb the minus sign into the phase; in this way 6(£) + 0 as E > 0. Note that 6(£) is positive 
and does not exceed 7/2. In fact a sketch of 6(£) is given in Figure 5. 


Figure 5: The phase 6(£) as a function of energy E < Vo. 


The total wavefunction for x < 0 is interesting 
w(a) = Ae 4 (— Aci) e-ike 
= AcE) (eW5(E) git = iS {E) g-ike) (3.23) 
= 2iAe®™) sin(ka — 5(E)) 
This means that the probability density is 


|p|? = 4A? sin?(ka — 5(E)). (3.24) 


The point x9 > 0 determined by the condition kao = 6(£) is the point in the forbidden region where 
the extrapolation of the allowed-region solution would vanish. Of course in the forbidden region x > 0, 
the probability density ||? is a decaying exponential. 


yy |> 
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Figure 6: Norm squared for the energy eigenstate when EF < Vo. For x > 0 the probability density decays 
exponentially with «. The point xo is the point where the extrapolation of the x < 0 probability density would 
have vanished. 


For future use we record the derivative of the phase 6(/) with respect to the energy 


fe ay “A 1 
§(B) = = a B=) (3.25) 


Note that this derivative becomes infinite both for EF — 0 and for E > Vo. 


= 


Figure 7: The derivative 6’(E£) as a function of energy E < VY. 


4 Wavepackets in the step potential 


Now we examine the more physical scenario. As we’ve seen with the free particle, the stationary states 
are not normalizable, and physical particles are actually represented by wavepackets built with an 
infinite superposition of momentum eigenstates. We can do similarly thing with our energy eigenstates. 
We will consider energy eigenstates with E > Vo or equivalently with k? > k2, where 


_— 2mME _ %mV 
a ae 
and we will superpose them. To begin we write the energy eigenstates in a slightly different form, 


including the time dependence. Setting A = 1 and using the values for the ratios B/A and C/A we 
find the solution 


k? =), (4.26) 


. (cit a ae aa) CE Oe 4.27 
0s 2k pike ¢—iE(k)t/h z>0 = 
kk . 


We can form a superposition of these solutions by multiplying by a function f(k) and integrating 
over k 


fee dk f(k) s te bok ee Pet, “<0, 


West) = i 
he dk f(k as elke e tB(R)E/h a>. 


(4.28) 


Here f(k) is a real function of & that is ok zero except for a narrow peak at k = ko. Note that 
we have only included momentum components with energy greater than Vo by having the integral’s 
lower limit set equal to k. The integral only runs over positive k because only in that case the e’** waves 
are moving towards positive x, and are therefore genuine incident waves. The above is guaranteed to 
be a solution of the Schrédinger equation. 


We can split the solution into incident, reflected and transmitted waves, as follows. 


Weel Dit) eer et) L< 0, 
U(z,t) = (4.29) 
Wianales t) ; xz>0. 


Naturally both Vj,-(x,t) and W;,.f(x,t) exist for « < 0 and Wyans(x,t) exists for x > 0. We then 
have, explicitly 


Wine(x < 0, t) 


i’ dk f(b te ee. 
k 


Vrep(a <0,t) = is dk f(k (FF p erttre tm (4.30) 


k+k 


Vaal > 04) = i, dk f (k) (7, Jetirersetnn, 
i 


How does the peak of Vj,-(a,t) move? For this we look for the main contribution to the associated 
integral which occurs when the total phase in the integrand is stationary for k ~ ko. We therefore 


require 
d hk? ¢ 
ka 
dk 2m h/ |, 


This is the relation between t and x satisfied by the peak of W;,,. It describes a peak moving with 
constant velocity hkg/m > 0. Since V;,-(a,t) requires that x < 0, the above condition shows that we 
get the peak only for t < 0. The peak of the packet gets to x = 0 at t = 0. For t > 0, Vine (a,t) is not 
zero, but it must be rather small, since the stationary phase condition cannot be satisfied for any x 
in the domain x < 0. 


h 
=0 > a“ Fk, _ 9 = pa (4.31) 
m m 


Consider now W,.f(x,t). This time the stationary phase condition is 
d wk 
FT alae 
m he 
The relation represents a peak moving with constant negative velocity —hko/m. Since W;e¢(x, t) 
requires that x < 0, the above condition shows that we get the peak only for t > 0, as it befits a 


reflected wave. For t > 0, Vpef(x,t) is not zero, but it must be rather small, since the stationary phase 
condition cannot be satisfied for any x in the domain x < 0. 


hk hk 
=0 > e+—t=0 = r= -—t. (4.32) 
m m 


Finally, let us consider W,gn;. The stationary phase condition reads: 


a (2 Ark? t dk hko 
dk (ix 2m =) és "Tk lko” om 
Using 7 
- 2mVo dk ek 
ie == 
fe” dk Ok’ 
back to the earlier equation we quickly find that 
hk 
Transmitted wave peak: x = = is 


(4.33) 


(4.34) 


(4.35) 


with k evaluated for k = kp. Since x > 0 is the domain of Virans this describes a peak moving to the 
right with velocity hk/m for t > 0. For t < 0, Wtrans(z,t) is not zero, but it must be rather small, 


since the stationary phase condition cannot be satisfied for any x in the domain x > 0. 


In summary, for large negative time V;,,, dominates and both V,¢f and Vtrans are very small. For 
large positive time, both V,¢f and WVt;-ans dominate and V;,, becomes very small. These situations 
are sketched in figures 8 and 9. Of course for small times, positive or negative, all three waves exist 
and together they describe the complex process of collision with the step in which a reflected and a 


transmitted wave are generated. 


Figure 8: At large negative times an incoming wavepacket is traveling in the + direction. 


~"S-- se - > 
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Figure 9: At large positive times we have a reflected wavepacket traveling in the —% direction and the trans- 


mitted wavepacket traveling in the + direction. 


Let us now examine a wavepacket built with energies EF < Vo. Recall that in this situation 
B/A= —e%(E) Therefore for an incident wave, all of whose momentum components have energy less 


than Vo, 


i 
taiesoe) = i dk f(k) ett eiBt/ (4.36) 
0 
the associated reflected wavefunction is 
i 
Vrep(a <0,t) = -{ de fhe ee EN (4.37) 
0 
Using the method of stationary phase again to find the evolution of the peak, 
d Et h2ko hkot 
— | 26(E) — kx — — = 26'(E = 0. 4. 
al (E) x ain => 0(E)— x m™ 0 (4.38) 
From this we quickly find 
ayia (¢ -— 2h6'(E)), (4.39) 
m 


where the derivative is evaluated at E(ko). The reflected wave packet is moving towards more negative 
x as time grows positive. This is as it should. But there is a time delay associated with the reflected 
packet, evident when we compare the above equation with x = — fko¢, The time delay is given by 


time delay = 2h6'(E). (4.40) 


The derivative 6’(E) was evaluated in (3.25) and it is positive. We see that the delay is particularly 
large for wave packets of little energy or those with energies just below Vo. 


We conclude the analysis of the step potential by discussing what it means to observe the particle 
in the forbidden region. It would be contradictory if the observer could make the following two 
statements: 


1. The particle is located in the forbidden region. 
2. The particle has energy less than Vo. 


Both statements taken to hold simultaneously would imply that the particle has negative kinetic 
energy, something that is inconsistent. In particular with E < Vo we would have a negative kinetic of 
magnitude Vo — E. 


E<V, 


X=0 
Figure 10: The step potential with potential energy Vo. If we could observe a particle in the forbidden region 
with energy F then the kinetic energy would be negative. 


First note that in the solution the particle penetrates into the forbidden region a distance of about 


1/K, where, you will recall that 
2m(V — E) 
n= = 9 (4.41) 
To be sure the particle is in the forbidden region its position uncertainty Ax must be smaller than the 


penetration depth: 


Kee >, (4.42) 
K 


The particle acquires some momentum p due to the position measurement: 


> — > Ar. 4.4 
a a = 7 ee) 


Due to this momentum induced by the position measurement there is some additional contribution E’ 


to the kinetic energy 
2 F262 
fe SoS a Ee (4.44) 


E' = 
2m 2m 


where we used (4.41). From this inequality we find that the total energy will exceed Vo 
Eig = E+E’ > E+(Vo-E)=W. (4.45) 


While the argument is heuristic, it gives some evidence that no negative kinetic energy will be detected 
for a particle that is found in the forbidden region. 


Sarah Geller transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this docu- 
ment. 
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1 Resonant transmission in a square well 


Consider the finite square well 


V(z) = 


f 
0, or |x| >a, (1.1) 


—Vo, for |z| <a. 


Here Vo > 0 has units of energy. We consider an energy eigenstate, a scattering solution with E > 0 
that represents an incoming wavefunction approaching the well from the left. An ansatz for the 
eigenstate will be of the form 


Ae*? +. Bethe L<-a 
w(x) =< Cel? 4 Det, |zl<a , (1.2) 
Fe. L>a. 


Ael™,,. pol | ad ar aes 


— 


Figure 1: The square well with all the relevant waves as one, with coefficient A is incident from the left. 


Here A is the coefficient of the incoming wave that exists for « < a, and B is the coefficient of the 
reflected wave (see Figure 1). Both of these waves have wavenumber k. In the well region |x| < a we 
have a wave moving to the right, with coefficient C and a wave moving to the left, with coefficient D. 
The wavenumber in this region is called ko. To the right of the well we have just one wave, moving 


to the right. Its coefficient is F' and its wavenumber k. Note that even though the potential is an 
even function of x, a non-normalizable energy eigenstate need not be even nor odd. The symmetry is 
broken by the condition that a wave is incident from the left. The values of k and k’, both positive, 
are determined by the Schrodinger equation and are 

2mE 9 2m(E + Vo) 

There are four boundary conditions: continuity of ~ and yy’ at « = —a and at x = a. These four 
equations can be used to fix the coefficients B,C, D, and F' in terms of A. We define reflection and 


= 


transmission coefficients R and T as follow 
_ |BP _ |FP 


— |Al?’ — | Al? ° (1.4) 


From probability current conservation we know that the currents to the left and to the right of the 
well must be equal so that 


JA? — |B? = |FP. (1.5) 
This is not an independent equation; it must follow from the boundary conditions. It implies that 
BP FP 
R+T= =); 1.6 
lap * Tap a 


showing that our definition of R and TJ’ makes sense. 
Solving for R and T is straightforward but a little laborious. Let us just quote the answer that 

one gets. The transmission coefficient is the following function of the energy F of the eigenstate: 
1 1 Ve 


en 


Tr 1 E(B +Ve) sin?(2k2a). (1.7) 


Since the second term in the right-hand side is manifestly positive, we have T < 1. As EF — 0, we 


have z —+ 1+ 00, which means T > 0. As E > ov, we have T > 1. 


We can remove all units from this result by defining 


B 9 2ma?Vo 
aie eg = ; 
Vo he? 
Then, 
2ma?(E + Vo) 2ma?Vo 


(kya)? = 7 = (lt+e) > 2k,a = 2/1 +e, (1.9) 


so we have 


1 
= 14 in? (22/1 : 1.10 
4 4e(1 +e) smal 20y Lae) ( ) 


We now can see that the well becomes transparent, making 7’ = 1, for certain values of the energy. 
All we need is that the argument of the sin function be a multiple of 7: 


2eVl+e = na, neZ. (1.11) 
Not all integers are allowed. Because e > 0, the left-hand side is bigger than or equal to 2z) and 
therefore ; 
ras) (1.12) 
T 


Call the associated energies E;, = €,Vo. Then 


nn? nen? Kh 
— = 1.13 
€n + 422 2m(2a)2Vo ’ ( ) 
so that 2272 
vn 
E VY = —~. 1.14 


Note that E, + Vo is the energy of the scattering state measured with respect to the bottom of the 
square well. The right-hand side is the energy of the n-th bound state of the infinite square well of 
width 2a. We thus have a rather surprising result: We get full transmission for those energies E, > 0 
that are in the spectrum of the infinite-square-well extension of our finite square well. The inequality 
n> a guarantees that E,, > 0. Since infinite square well bound states are characterized by fitting an 
integer number of half wavelengths, we have a resonance type situation in which perfect transmission 
is happening when the scattering waves fit perfectly inside the finite square well. The phenomenon we 
have observed is called resonant transmission! 


Figure 2: We get resonant transmission across the finite well at the positive bound-state energies of a would-be 
infinite well. 


The fitting of an exact number of half wavelengths can also be seen directly from the vanishing of 
the sin function in (1.7) giving k2(2a) = nz which implies 
2 
(2a) =nt + 2% = nd. (1.15) 


We show in Figure 3 the transmission coefficient T as a function of e = E’/Vo for a square well with 
zq = 137/4. In this case we must have n > iu or n > 7. 


2 The Ramsauer—Townsend Effect 


Carl Ramsauer and John Sealy Townsend published separately their investigations in 1921. They 
were studying the elastic scattering of low-energy electrons off of rare-gas atoms, such as noble gases. 
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Figure 3: The transmission coefficient T as a function of e for a square well with zp = Bn At the energies for 
which T = 1 we have resonant scattering. The values of E7, Eg,... are indicated in the figure. There are three 


cases of resonant transmission for 0 < E < Vo. Note that the spacing between the points where T = 1 grows as 
€ grows. 


These gases have their electronic shells fully filled and are both very un-reactive and have high ion- 
ization energies. The potential is created by the nucleus, and becomes visible as the incident electron 
penetrates the electron cloud. This potential is a spherically symmetric attractive potential for the 
electrons — some kind of finite spherical well. In the experiment some electrons collide with the atoms 
and scatter, mostly bouncing back. We can thus view the reflection coefficient R as a proxy for the 
scattering cross section! 

Ramsauer and Townsend reported a very unusual phenomenon. At very low energies, the scattering 
cross section was high. But the energy dependence was surprising. As the energy was increased the 
scattering went down towards zero, to go up again as the energy was increased further. Such mysterious 
behavior had no sensible classical explanation. What one has at play is quantum-mechanical resonant 
scattering. The cross section going to zero means the reflection coefficient going to zero, and the 
transmission coefficient going to one! The first resonant transmission happens for electrons of around 
one electron-volt (such electrons have a speed of about 600 Km/s). Figure 4 provides a sketch of 
both R and T as a function of energy. Our one-dimensional square well potential does not provide a 
good quantitative match for the data, but it illustrates the physical phenomenon. A three-dimensional 
spherical square well is needed for a quantitative analysis. 


Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this 
document. 


Figure 4: The reflection and transmission coefficients as a function of energy for the Ramsauer—Townsend 
effect. Note that R+ T= 1. At the arrow, we have R = 0, so there is no scattering. All of the electrons go 
straight through the noble gas atoms! They experience resonant transmission. This first occurs around leV. 
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1 Scattering in One Dimension 


Physicists learn a lot from scattering experiments. Rutherford learned about the structure of the atom 
by scattering alpha particles off thin layers of gold. We have elastic scattering if particles do not change 
type; this typically requires low energies. At high energies, scattering becomes quite complicated due 
to the creation of new particles. 

The scattering of a particle off a fixed target is studied by working with a particle moving in a 
potential, the potential created by the target. Even in the case of particle collisions it is generally 
possible to study the problem using center of mass coordinates and, again, considering scattering in a 
potential. 


Figure 1: A potential of range R. The potential vanishes for 7 > R and is infinite for x < 0. 
We consider (elastic) scattering in the general potential shown in Fig. 3. This potential is given by 


Vix) O<a<R, 
Vie =< 0 z>R, (1.1) 


OO xa<O. 


We call this a finite range potential, because the nontrivial part V(x) of the potential does not extend 
beyond a distance R from the origin. Moreover, we have an infinite potential wall at « = 0. Thus 
all the physics happens in x > 0 and incoming waves from x — oo will be eventually reflected back, 
giving the physicist information about the potential. The restriction to x > 0 will also be useful when 


we later study the more physical case of scattering in three dimensions. In that case, when using 
spherical coordinates, we have r > 0 and most of the concepts to be learned here will apply. 


Consider first the case where V(x) = 0, so that 


V(x) = ° . : (1.2) 


as shown in figure 2. We have a free particle, except for the wall at x = 0. The solutions can be 


constructed using linear combination of momentum eigenstates e+“** where 
2nE 
2 
— 72 (1.3) 


Figure 2: The zero potential and its energy eigenstates ¢(a) = sin ka are needed to compare with the V(x) 4 0 
problem. 


The solution has an incoming wave e~*** and an outgoing wave e’** combined in such a way that 


(0) = 0 as required by the presence of the wall: 
b(a) ~ et? — ete, (1.4) 


Better yet, we can divide this solution by 27 to find 


The incoming wave is the first term to the right of the first equal sign and the reflected wave is the 
second term. Both carry the same amount of probability flux, but in opposite directions. 


Now consider the case V(x) 4 0. This potential always acts over the finite range 0 < x < R and 
we will eventually be interested in computing the energy eigenstate wavefunction 7(x) in this region. 
For the moment, however, let us consider 7(x) in the region x > R. We will take the incoming wave 
to be the same one we had for the zero-potential solution ¢(x): 


eke 


24 


Incoming wave: = (1.6) 
The outgoing wave to be added to the above requires an e’*, to have the same energy as the incoming 
wave solution. We now claim that the most general solution includes a phase factor, so we have 


ka 


of 
Outgoing wave: ee = , OER. (1.7) 


We note that the phase 6 cannot be a function of x: the free Schrédinger equation for x > R only 
allows phases linear in x, but that would change the value of the momentum associated with the 
outgoing wave, which we have already argued must be k. Furthermore, 6 cannot be complex because 
then we would fail to have equality between incident and reflected probability fluxes. This condition 
requires that the norm-squared of the numbers multiplying the exponentials e+’** be the same. Thus, 
6 is a real function that depends on the energy FE of the eigenstate and, of course, of the potential 
V(x). A natural range of 6 is from zero to 27, but it will be easier to let it range through all of the 
reals R in order to achieve a 6 that is a continuous function of the energy. Assembling together the 
incident and reflected components of the x > R solution we get 


4 (-e + genre | = esin(kx +6), for «>R. (1.8) 


This is called the canonical solution for x > R. For any particular feature of the solution ¢(x) that 
we might find at kx = ao, giving x = ag/k, we would find the same feature in y)(x) at ki + 6 = ao, 
giving & = ao/k — 6/k. For small 6 > 0, the wave is pulled in by an amount 6/k, and the potential is 
exerting attraction. For small 6 < 0, the wave is pushed out by an amount |6|/k, and the potential 
is exerting repulsion. Note also that 6 and 6 +7 give exactly the same 7(x). This is simplest to see 
from the first form in (1.8). 


P(x 


Figure 3: The solution ¢(a) for zero potential is shown with a dashed line. For x > R, the solution w()is 
shown as a continuous line. As compared to ¢(x), it is spatially shifted towards the origin by a distance 6/k. 


We define the scattered wave 1,(x) as the extra wave in the solution ~(x) that would vanish for 
the case of zero potential, i.e. 
P(x) = o(x) + Ys(2). (1.9) 
Note that both @ and w have the same incident wave, so wy, must be an outgoing wave. We find, 
indeed 


etka eika+2i6 e tke eikx eikax+2i6 elke 
= — ——_ - = —__ — 1.10 
U(0) = ¥e)-6@) = -S- +$ S45 - GS eS) 
and therefore we get 
w(x) = e*sinde™™ = A,e**, with A, = e sind. (Lat) 


A, is called the scattering amplitude, being the amplitude for the scattered wave. While we can’t 
normalize our states yet (for that we need wave packets) the probability to scatter is captured by 


|A,|? = sin? 6. (1.12) 


1.1 Time Delay 


The quantity 6(£) determines the time delay of a reflected wave packet, as compared to a wave packet 
that would encounter a zero potential. Indeed we claim that the delay is given by 


A= 2nd (2). (i133) 


where prime denotes derivative with respect to the argument, and the evaluation is done for the 
central energy of the superposition that builds the wave packet. If At < 0, then the particle spends 
less time near x = 0, either because the potential is attractive and the particle speeds up, or because 
the potential is repulsive and the particle bounces before reaching « = 0. If At > 0, then the particle 
spends more time near x = 0, typically because it slows down or gets temporarily trapped in the 
potential. 

We write the incoming wave in the form 


Vine(ast) = [die g(kje eI, a> R, (1.14) 
0 


where g(k) is a real function peaked around k = kg. We write the associated reflected wave by noticing 
that the above is a superposition of waves as in (1.6) and the reflected wave must be the associated 
superposition of waves as in (1.7). We must therefore change the sign of the momentum, change the 
overall sign and multiply by the phase e?”: 


Weanle,f) = -{ dk g(k)e** e25() ett og > R, (1.15) 
0 


We now use the stationary phase approximation to figure out the motion of the wave packet peak. We 
must have a stationary phase when k = ko: 


d E(k)t 
ae 9 papas iar a 
0 Tk (ie + 0(k) = ) a 
ee ce 
ET Ele, AR |g B 
pee 2 ge) Se 
dk |, FE | gre) OE [gy (1.16) 
1 dE di 
2642 ae rs 
h dk | 71 re ) 
hko /», 46 
= aL Fy pea = 
: m ( dE | ekg) i), 
oo hk dé 
= —(t—2h basal (1.17) 
w dE | 13(40) 


Here kg/m is the familiar group velocity vp of the wavepacket. If there had been no phase shift 6, 
say because V(a) = 0, there would be no time delay and the peak of the reflected wavepacket would 
follow the line x = vot. Therefore the delay, as claimed, is given by 


At = 2h6'(E(ko)). (1.18) 


We can compare the time delay to a natural time in the problem. We first rewrite 


ee 0.9 
We know that 1ldE dw 
p= oe = tis; (1.20) 
so we can then write dé5 At 
= (1.21) 


Multiplying this by # we have 


(1.22) 


The left-hand side is unit free, and the right hand side is the ratio of the time delay to the time the 
free particle would take to travel in and out of the range R. 


1.2. An Example 


Consider the following example in which we have an attractive potential: 


—-Vo, for 0<a2<a, 
VQS=s 2, tor 2a, (1.23) 


coo, for x<O0. 
The potential is shown in Figure 4 where the energy F > 0 of the energy eigenstate is indicated with 
a dashed line. The solution is of the form 


e*sin(ke +6) 2>a, 


V= Shes UA 2 Gs ee 


Figure 4: An attractive potential. We are interested in calculating the energy eigenstates for all E > 0. 


The solution in the region x > a is simply the canonical solution, and the solution in the region x < a 
follows because we must have w(0) = 0. The constants k and k’ are given by 


2nE p2 — 2m(E + Vo) 


Matching 7 and yw’ at x = a we find the conditions that will eventually give us the unknown phase 
shift 6: 


Asin(k’a) = e® sin(ka + 6) (1.26) 
k' Acos(k’a) = ke® cos(ka + 6). (1.27) 


Dividing the latter equation by the former, we reach 


kcot(ka + 6) = k' cot k’a. (1.28) 
We now make use of the identity 
cot Acot B—1 
t(A + B) = ———————_ 1.29 
conan) cot A+ cot B’ nz) 
from which we find re ‘ a 
cot ka cot 0 — 
— cot k'a = cot(ka + 6) = —————. 1.30 
hae cone te) cot ka + cot 6 ( ) 
Solving this equation for cot 6 yields 
tan ka + © cot k’ 
(ie ee (1.31) 


1- F cot k'atan ka 


While this is a complicated formula to analyze directly, we can always plot it with a computer for 
different values of Vo. As usual we characterize the well by the constant zo, defined by 


2mVoa? 
es ie (1.32) 
Note also that 
ka = /224+(ka)?. (1.33) 


Figure 5 shows the phase factor 5, the quantity sin? 6, the time delay 16 and the amplitude |A| 
inside the well as functions of ka, for ze = 3.40. 


Note that the phase 6 begins at 0 for zero energy and reaches —7 for infinite energy. The excursion 
of the phase is thus 7 and as we will see, this happens because for this value of zo the potential would 
have one bound state. 

The value of |A,|? = sin? 6 represents the probability of scattering and it peaks for the value of ka 
for which the phase 6 has absolute value 7/2. 

Next in the plot is the unit-free delay 14 
explained: as the particle moves over the well, its kinetic energy is increased by Vo, the particle speeds 


up as it reaches and bounces off the wall. 


Notice that the delay is negative. This is easily 


The last plot shows the magnitude |A| of the constant that gives the amplitude of the wavefunction 
in the region 0 < x <a. 


At very large energy FE >> Vo, the particle hardly notices the potential. Indeed we see that 


6 approaches —z which, as noted below (1.8), is equivalent to 6 = 0 and means no phase shift. 


Accordingly sin? 5 + 0, meaning no scattering. We also hav 4 dp — 0, meaning no time delay and, 


finally, |A] — 1 as in the free solution. 


Figure 5: Various quantities plotted as functions of ka, for 22 = 3.40. Top: the phase shift, going from zero to 
—n. Second: the scattering amplitude |A,|? = sin? 6. Third: the delay relative to the free transit time. Last: 
the norm |A| of the amplitude of the wave inside the well. 
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1 Levinson’s Theorem 


Levinson’s theorem relates the number Ny of bound states of a given potential to the excursion of the 
phase shift 6(£) as the energy goes from zero to infinity: 


Ny = =(6(0) ~6(00)), (1.1) 


To prove this result consider an arbitrary potential V(x) of range R, with a wall at x = 0. This 
potential, shown to the left in Figure 1, has a number of bound states, all of which are non-degenerate 
and can be counted. There is also a set of positive energy eigenstates: the scattering states that, 
belonging to a continuum, cannot be counted. Our proof requires the possibility of counting states, 
so we will introduce a second in nite wall, placed at x = L for large L. Of course, this will change the 
spectrum, but as LZ becomes larger and larger the changes will become smaller and smaller. We think 
of L as a regulator for the potential that discretizes the spectrum and thus enables us to enumerate 
the states. It does so because with two walls, the potential becomes a wide infinite well and all states 
become bound states. The potential with the regulator wall is shown to the right in Figure 1. 


V&) 4\ 60) 


Figure 1: Left: An arbitrary one-dimensional potential V(a) of range R. Right: The same potential with a 
regulator wall placed at x = L. 


The key of the proof will be to compare the counting of states in the regulated V 4 0 potential to 
the counting of states in the V = 0 potential, also regulated with a second wall at « = L. Consider 
therefore the regulated V = 0 potential and the positive energy eigenstates. These correspond to the 
wavefunction $(x) = sinkx, with the second wall requiring ¢(a = L) = 0. We thus have 


kL=nnr, with n=1,2,... (1.2) 


The values of & are now quantized. Let dk be an infinitesimal interval in wavenumber with dn the 
number of states in dk when V = 0. Thus, 


daL=a7c > d= 2 a. (1.3) 
7 


OO 


dh x 


Figure 2: With the regulator wall the wavenumber k takes discrete values. dk is an infinitesimal interval in k 
space. 


When V(x) 4 0, the solutions for x > R, all of which are positive energy solutions, have the form 
w(x) = e® sin(ka +6). (1.4) 

The boundary condition w(L) = 0 implies a quantization 
kL +6(k) = nn, (1.5) 


with n’ integer. We can again differentiate to determine the number of positive energy states dn’ in 
the interval dk, with V # 0: 
dé — ee! 1 (dé 
dkh + 7 dks =dnx —> dn = Zak + =(S)ae. (1.6) 


The number of positive energy solutions lost in the interval dk as we turn on the potential V is given 
by dn — dn’, which can be evaluated using (1.3) and (1.6): 


dn—dn! = -+ (Gat). (7) 


The total number of positive energy solutions lost as the potential V is turned on is given by integrating 
the above over the full range of k: 


a | 
# of positive energy solutions lost as V turns on = -{ ao = ——(d(co) —4(0)). (1.8) 
0 Tv 


Figure 3: The positive energy states of the V = 0 setup shift as the potential is turned on and some can become 
bound states. 


Although we lose a number of positive energy solutions as the potential V is turned on, states do 
not disappear. As one turns on the potential from zero to V continuously, we can track each energy 
eigenstate and no state can disappear! If we lose some positive energy states those states must now 
appear as negative energy states, or bound states! Letting N, denote the number of bound states in 
the V 4 0 potential, the result in (1.8) implies that 

N= OCs, (1.9) 


T 


This is what we wanted to prove! 


2 Resonances 


We have calculated the time delay At = 2h6’(E) associated with the reflected wavepacket that emerges 
from the range R potentials we have considered. If the time delay is negative, the reflected wavepacket 
emerges ahead of time. We can ask: Can we get an arbitrarily large negative time delay? The answer 
is no. A very large time delay would be a violation of causality. It would mean that the incoming 
packet is reflected even before it reaches « = R, which is impossible. In fact, the largest negative time 
delay would be realized (at least classically) if we had perfect reflection when the incoming packet hits 
x = R. If this happens, the time delay would be Se where vo is the velocity of the packet. Indeed, 
2R is the time saved by the packet that did not have to go in and out the range. Thus we expect 


vO 


dé 2k 
i ] = 2h— > —-—. 2.1 
time delay dE =~ % (2.1) 
This can be simplified a bit by using k derivatives 
dé 1 dé 2 dé 2k 
2h = 2h = > 2.2 
dE a dk vodk ~ v9 ; ( ) 
which then gives the constraint 
dé 
as ee 2.3 
dk ~ 22) 


The argument was not rigorous but the result is rather accurate, receiving corrections that vanish for 
packets of large energy. 

Alternatively, we can ask: Can we get an arbitrarily large positive time delay? The answer is 
yes. This can happen if the wave packet gets temporarily trapped in the potential. In that case we 
would expect the probability amplitude to become large in the 0 < x < R region. If the wavepacket 
is trapped for a long time we have a resonance. The state is a bit like a bound state in that it gets 
localized in the potential, at least for a while. In order to get a resonance it helps to have an attractive 
potential and a positive energy barrier. We can achieve that with the potential 


ee for «<0 

Vey —Vo for 0<a<a (2.4) 
Vi for a<a< 2a 
0 for x > 2a. 


The potential, with Vo,Vi > 0, is shown in figure 4. In order to have a resonance we take explore 
energies in the range zero to Vj. In such range of energies we can expect to find some particular values 


that lead to resonant behavior, namely, large time delay and large amplitude for the wavefunction in 
the well. 
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Figure 4: We search for resonances with energy F in the range (0,V;). In this range the V; barrier produces a 
classically forbidden region x € (a, 2a), that can help localize the amplitude around the well. 


Given the three relevant regions in the potential we define 


2 —E 2nE 
yr = E+ VO) cea mF) i2 = 7 (2.5) 


In the region 0 < x < a we must use trigonometric functions of k’x. In the region a < x < 2a we use 
hyperbolic functions of ka and in the region x > 2a we use the canonical solution with phase shift and 
wavenumber k. In the middle region a < x < 2a we could use a combination of solutions 


{e**, e **}, or {coshaxx, sinhxx}, or {coshx(x—a), sinhk(x—a)}. (2.6) 


The last pair is most suitable to implement directly the continuity of the wavefunction at x = a. So 
we can write for the wavefunction w(x): 


Asin(k’2) 0<a<a 
w(x) = 4 Asin(k’a) cosh«(z —a) + Bsinhk(x—a) a<a< 2a (2.7) 
e sin(ka + 6) ze > 2a 


After implementing the remaining boundary conditions we can solve for the phase shift 6. After a 
modest amount of work one finds: 


. iF . 
ka sink’a coshka + E cos k’a sinh ka 


tan(2ka+ 6) = z . 2.8 
( ) Ka sink’a sinhka + E cos k’a cosh Ka ve) 
This expression is fairly intricate so it is best to do numerical work. For this we define 
2mVoa? WMV) a? 
a. = 35.3 a = 7 oe eS ka. (2.9) 
which allow us to express both k’a and «a as functions of u 
(kia)? = 2+u?, (xa)? = 22-u?. (2.10) 


At this point (2.8) can be used to determine 6 as a function of u = ka and the constants zo, 21. Suppose 


we pick values for our parameter controlling equations. In Figure 5 we show results for z{ = 1 and 
2 
zy = 5. 
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Figure 5: Plot of various quantities as a function of u = ka, with the potential characterized by 22 = 1 and 
2 =5. (a) 6(£) increases quickly around us = 1.85, or equivalently E = 0.69V;, crossing —7/2 and signaling 
resonant behavior. (b) Plot of |A,|? = sin? 5, showing peaks each time |6| = 1/2. (c) The coefficient |A| of the 
wavefunction at the well peaks at the resonance, showing high probability of finding the particle at the well. 
(d) The time delay is positive and peaks at resonance. 


Consider part (a) of the figure, showing 6(ka). At the beginning 6 decreases linearly, a sign of a 
negative time delay, as the low energy waves reflect at the edge x = 2a of the V; barrier. As 6 crosses 
—r/2 there is no resonance, even though |A,|? = sin? 6 is equal to one. Indeed we see no bump in the 
amplitude |A|. As the energy is increased and u = ux = 1.8523 we get a resonance. This time 6 is 
increasing rapidly and 6 crosses —7/2 again, making |A,|? = 1. The signal of resonance is the very 


high |A| the peak in the time delay. This time delay reaches the value of about 14, meaning the delay 
is fourteen times the free transit time 4a/vo! 


3 Modeling the Resonance 


We would like to have further insight into the nature of resonances. In particular we want to appreciate 
the general features of the phenomenon. Additionally, so far we can identify resonances by looking at 
the behavior of 6 but, can we find an equation that defines resonances? 

As a first step, we model the behavior of a phase near resonance. Recalling that a resonance 
requires |6| cross the value 7/2 and that 6, physically, is the same as 6 increased or decreased by 
multiples of 7 we can choose to have 6 vary from nearly zero to nearly 7. We can achieve this with 
the following simple function. 


6=tan! (4) : with B>0, a>O. (3.1) 
a-—k 

here a and £ are positive constants with the same units as k. To see what this function does, we first 

plot the argument of the arc tangent at the top of Figure 6. Note that the argument varies quickly in 

the region (a — 8,a+ 3). The variation of the associated phase 6 is shown in the figure below. To 

have a sharp increase in the phase we must have small 6 compared to a. 


Figure 6: The constant 6 must be small compared to a to obtain a sharp variation. A resonance, as shown 
here requires 6 increasing with energy. 


Two relatively short calculations give us further insight: 


dé 1 oo ae 
dklk=a 8’ al" = sins = 


pe 


The first one informs us that, all things being equal, the delay is large if G is small. The second gives 


(3.2) 


the norm-squared of the scattering amplitude as a function of k, with a peak at k = a. This equation 


is most famously expressed in terms of the energy. For this we note that 
h2 2 h2 
k+a)(k-—a) & Dm LON — &) (3:3) 
m 


a 


when working with k = a. It thus follows that 


and therefore ' a 
B “Tr 
es? = 2 We. 2 = 1 2 ipo? (3.5) 
B* + 7a (E — Ea) (#— Eq)? + 4r 
Where we have defined the constant I with units of energy: 


ip or 
2 


_ 2aBh? 
r= ae (3.6) 


1yp2 
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The energy dependence of |w,|* follows the so-called Breit-Wigner distribution, 


1p2 
al 


2 ~~ 
IPs ~ (E- Ey)? + 41? 


(3.7) 


The distribution is shown in Figure 8. The peak value for |2),|? is attained for E = Ey and is one. 
We call T’ the width at half-maximum because the value of |7)5|? at E = Eq + 5I is one-half. Small T 


corresponds to a narrow width, or a narrow resonance. 


patie ee $- ee 
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Figure 7: The Breit-Wigner distribution. T is the width of the distribution at half maximum. 


To understand better the significance of I we define the associated time 7, called the lifetime of the 


resonance: h 
m 
= = . 3.8 
"= T ~ dash 8) 


As you probably would expect, the lifetime is closely related to the time delay associated with a 
wavepacket of mean energy equal to the resonant energy. Indeed, we can evaluate the time delay At 
for k = a to get 


dé dk dé 2h (1 2h 2m 
At = 2h— = 2h— — = = = =A4r. : 
iB = "a6 Wa ~ E(B) = Ea ~ aah a 


We therefore conclude that the lifetime and the time delay are the same quantity, up to a factor of 
four. 


= ;At. (3.10) 


Unstable particles are sometimes called resonances. The Higgs boson, discovered in 2012, is an 
unstable particle with mass 125 GeV. It can decay into two photons, or into two tau’s, or into a bb 
pair, among few possibilities. The width [ associated to the particle is 4.07 Mev (44%). Its lifetime 
T is about 1.62 x 10~?? seconds! 


We now try to understand resonances more mathematically. We saw that, at resonance, the norm 
of A, reaches a maximum value of one. Let us explore when A, is large. We have 
sin 6 sin 6 tan é 


A, = sinde® = —_— = = — 3.11 
i ds etd cos 6 — isin d 1—itand ( ) 


At resonance 6 = 7/2 and A, = i, using the first equality. On the other hand, while we usually think 
of 6 as a real number, the final expression above indicates that A, becomes infinite for 


tand = -i, (3:12) 


whatever that means! If we recall that taniz = 7itanh z we deduce that the above condition requires 
6 — —ioo, a rather strange result. At any rate, A, becomes in nite, or has a pole, at tand = —i. We 
will see that the large value |A;| = 1 at resonance can be viewed as the “shadow” of the infinite value 
A, reaches nearby in the complex plane. 
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Figure 8: In the complex k plane, resonances are identified as poles of the scattering amplitude A, located 
slightly below the real axis. Bound states appear as poles on positive imaginary axis. 


Indeed, we can see how A, behaves near resonance by inserting the near-resonance behavior (3.1) 
of 6 into (3.11): 
or B 
As = s = ‘ 3.13 
"Tig Gs sy 


When k = a, meaning at the resonant energy, we get A, = 7, as expected. If we now think of the 
wavenumber k as a complex variable, we see that the pole of A, is a pole atk = k, = a—if. The 
real part of k, is the resonant energy, and the imaginary part 6 encodes the lifetime. For small 6 
the resonance is a pole near the real axis, as illustrated in Figure 8. The smaller 6 the sharper the 
resonance. As we can see, the value of |A,| on the real line becomes large for k = a because it is 
actually infinite a little below the axis. 

The lesson in all of this is that we can indeed take (3.12) seriously and look for resonances by 
solving for the complex k values for which 


Resonance condition: tand(k) = —i. (3.14) 


The real part of those k’s are the resonant energies. The imaginary parts give us the lifetime. 
The idea of a complex k plane is very powerful. Suppose we consider purely imaginary & values of 
the form k = ik, with « > 0. Then the energy takes the form 
hae 
E = —-— <0, (3.15) 
2m 
which is suitable for bound states. Indeed one can show that bound states appear as poles of A, along 
the positive imaginary axis, as shown in Figure 8. The complex k-plane has room to fit scattering 
states, resonances, and bound states! 


Sarah Geller transcribed Zwiebach’s handwritten notes to create the first LaTeX version of this docu- 
ment. 
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1 Schrodinger Equation in 3D and Angular Momentum 


We have so far considered a number of Hermitian operators: the position operator, the momentum 
operator, and the energy operator, or the Hamiltonian. These operators are observables and their 
eigenvalues are the possible results of measuring them on states. We will be discussing here another 
operator: angular momentum. It is a vector operator, just like momentum. It will lead to three 
components, each of which is a Hermitian operator, and thus a measurable quantity. The definition 
of the angular momentum operator, as you will see, arises from the classical mechanics counterpart. 
The properties of the operator, however, will be rather new and surprising. 

You may have noticed that the momentum operator has something to do with translations. Indeed 
the momentum operator is a derivative in coordinate space and derivatives are related to translations. 
The precise way in which this happens is through exponentiation. Consider a suitable exponential of 
the momentum operator: 7 

en , (1.1) 
where a is a constant with units of length, making the argument of the exponential unit free. Consider 
now letting this operator act on a wavefunction w(x) 


ipa 


er p(x) = etaey(a), (1.2) 


where we simplified the exponent. Expanding the exponential gives 


ipa d ad a® 

aa = (14 bee , 

eve) = (tase + ast gat +) YO) as 
_ dy | a? dy a? Bw = 
SEE Ge aaa Oe 


since we recognize the familiar Taylor expansion. This result means that the operator eh moves 
the wavefunction. In fact it moves it a distance —a, since (a + a) is the displacement of w(x) by a 
distance —a. We say that the momentum operator generates translations. Similarly, we will be able 
to show that the angular momentum operator generates rotations. Again, this means that suitable 
exponentials of the angular momentum operator acting on wavefunctions will rotate them in space. 


Angular momentum can be of the orbital type, this is the familiar case that occurs when a particle 
rotates around some fixed point. But is can also be spin angular momentum. This is a rather different 
kind of angular momentum and can be carried by point particles. Much of the mathematics of angular 
momentum is valid both for orbital and spin angular momentum. 

Let us begin our analysis of angular momentum by recalling that in three dimensions the usual # 
and operators are vector operators: 


f arr h hfo 0d 0 
P = (Bes PyPs) = a ~ (Sapa): 


Z (1.4) 
& = (4,9, 2). 
The commutation relations are as follows: 
[eae De | = th , 
[9, By] = ih, (1.5) 
[2, Dz | = ih. 


All other commutators are involving the three coordinates and the three momenta are zero! 


Consider a particle represented by a three-dimensional wavefunction ~(, y, z) moving in a three- 
dimensional potential V(r). The Schrédinger equation takes the form 


2 
5 Vv) + Vee) = Ev(r). (1.6) 


We have a central potential if V(r) = V(r). A central potential has no angular dependence, the 
value of the potential depends only on the distance r from the origin. A central potential is spherically 
symmetric; the surfaces of constant potential are spheres centered at the origin and it is therefore 
rotationally invariant. The equation above for a central potential is 


2 
VE) FEV PvE) = Bv(r). (1.7) 


This equation will be the main subject of our study. Note that the wavefunction is a full function of r, 
it will only be rotational invariant for the simplest kinds of solutions. Given the rotational symme- 
try of the potential we are led to express the Schroddinger equation and energy eigenfunctions using 
spherical coordinates. 


In spherical coordinates, the Laplacian is 


a ec ee ae ees, 1.04 i 
Viea (VV) aah) ha anéoe° 60 ante oe) ne) 


Therefore the Schrédinger equation for a particle in a central potential becomes 


r- sin 0 + 
r Or? r2 


sin 6 06 00 sin? 6 0¢? 


2m 


h? E Or. a 1 0 0 1 o>) e+ Verve ee (1.9) 


In our work that follows we will aim to establish two facts: 


1. The angular dependent piece of the V? operator can be identified as the magnitude squared of 
the angular momentum operator 


a) 1 oe L? 
ind = 1.10 
and 00° 0.’ sin? @ Og? ie (E10) 
where at 4,  % a 
L? = Ly Le + Lyly + Lhe. (iid) 


This will imply that the Schrédinger equation becomes 


h? [1 0? 1 L? 
r Or2 . r2 2 


Jurvine = Ey (1.12) 
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or expanding out 


h2 2 L2 
; +o (ry) tb + V(ry = Ey. (1.13) 
mr or 
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2. Eq. (1.7) is the relevant equation for the two-body problem when the potential satisfies 
V(ri,r2) =V (ri — ral), (1.14) 


namely, if the potential energy is just a function of the distance between the particles. This 
is true for the electrostatic potential energy between the proton and the electron forming a 
hydrogen atom. Therefore, we will be able to treat the hydrogen atom as a central potential 
problem. 


2 The angular momentum operator 


Classically, we are familiar with the angular momentum, defined as the cross product of r and p: 
L=r xp. We therefore have 
L = (Lz, Ly,L,)= rxp, 
Ly = YPz — 2Py; 
Ly = 2px — LPz , 
Lz = XPy — YPx - 


(2.1) 


We use the above relations to define the quantum angular momentum operator L and its components, 
the operators (Ly, Ly, Lz): 


(2.2) 


In crafting this definition we saw no ordering ambiguities. Each angular momentum operator is the 
difference of two terms, each term consisting of a product of a coordinate and a momentum. But note 
that in all cases it is a coordinate and a momentum along different axes, so they commute. Had we 
written Ly = pj — py2, it would have not mattered, it is the same as the Lz above. It is simple to 


check that the angular momentum operators are Hermitian. Take L,, for example. Recalling that for 
any two operators (AB)! = BI AT we have 


(Ex)t = (G6. — 26,)' = (G6.)t — (2by)t = ptgt — plat. (2.3) 


Since all coordinates and momenta are Hermitian operators, we have 


(Lz)! = $.9=pye = Gbe= by = Le; (2.4) 


where we moved the momenta to the right of the coordinates by virtue of vanishing commutators. 
The other two angular momentum operators are also Hermitian, so we have 


Dies Tes (2.5) 


All the angular momentum operators are observables. 

Given a set of Hermitian operators, it is natural to ask what are their commutators. This compu- 
tation enables us to see if we can measure them simultaneously. Let us compute the commutator of 
Ds with Le 

[Boas Ly] = [Wz = Zpy, 2px = De | (2.6) 
We now see that these terms fail to commute only because 2 and p, fail to commute. In fact the first 
term of L, only fails to commute with the first term of Lig Similarly, the second term of L, only fails 
to commute with the second term of i. Therefore 


ak, = [YPz, ZDz | + |Z 
2 | 
= 9 [pes Z] Po t+ & Les Dz| Dy (2.7) 


= ih(@py — GPx) - 
We now recognize that the operator on the final right hand side is L, and therefore, 
(Bar|, Sa: (2.8) 


The basic commutation relations are completely cyclic, as illustrated in Figure 1. In any commutation 
relation we can cycle the position operators as in > y7 > Z — # and the momentum operators as 
in py + Py + pz — py and we will obtain another consistent commutation relation. You can also see 
that such cycling takes L, > iy + L, — Lz, by looking at (2.2). We therefore claim that we do not 
have to calculate additional angular momentum commutators, and (2.8) leads to 


Baan oe = ih Lz , (2.9) 


This is the full set of commutators of angular momentum operators. The set is referred to as the 
algebra of angular momentum. Notice that while the operators L were defined in terms of coor- 
dinates and momenta, the final answer for the commutators do not involve coordinates nor momenta: 
commutators of angular momenta give angular momenta! The L operators are sometimes referred to 


as orbital angular momentum, to distinguish them from spin angular momentum operators. The spin 
angular momentum operators 5, oy and §, cannot be written in terms of coordinates and momenta. 
They are more abstract entities, in fact their simplest representation is as two-by-two matrices! Still, 
being angular momenta they satisfy exactly the same algebra as their orbital cousins. We have 


[ nim) al = ih S., 
[ Sis 5. | = ih Sy , (2.10) 
[Sz,52] = th Sy. 
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Figure 1: The commutation relations for angular momentum satisfy cyclicity). 


We have seen that the commutator [#,p] = ih is associated with the fact that we cannot have 
simultaneous eigenstates of position and of momentum. Let us now see what the commutators of L 
operators tell us. In particular: can we have simultaneous eigenstates of L, and oe. As it turns 
out, the answer is no, we cannot. We demonstrate this as follows. Let’s assume that there exists a 
wavefunction ¢9 which is simultaneously an eigenstate of Ly and Ly 


Lz bo = Axb0 5 
Ly ¢0 = Aydo- 


Letting the first commutator identity of (2.9) act on ¢9 we have 


(2.11) 


ihLzd0o = [Bess Ly ]¢o = Ly Ly do _ Ly Lego 
= LyAy Go _ LyAxo (2.12) 
= (cry = AyAx) Po = 0; 


showing that L.¢) = 0. But this is not all, looking at the other commutators in the angular momentum 
algebra we see that they also vanish acting on ¢p and as a result Az and , must be zero: 


[Ly Le]éo = thLedp = ihrAydo —> Xe = 0, 


MS is 2 2.13 
[ Lz, Ly |¢o0 = thLy do = th Ayoo => Ay = 0. ( ) 


All in all, assuming that ¢o is a simultaneous eigenstate of L, and Ly has led to [n¢) = Ly bo = 
Lz¢9 = 0. The state is annihilated by all angular momentum operators. This trivial situation is not 


very interesting. We have learned that it is impossible to find states that are nontrivial simultaneous 
eigenstates of any two of the angular momentum operators. 

For commuting Hermitian operators, there is no problem finding simultaneous eigenstates. In fact, 
commuting Hermitian operators always have a complete set of simultaneous eigenstates. Suppose we 
select L, as one of the operators we want to measure. Can we now find a second Hermitian operator 
that commutes with it? The answer is yes. As it turns out, L?, defined in (1.11) commutes with L, 
and is an interesting choice for a second operator. Indeed, we quickly check 


(2.14) 


So we should be able to find simultancous eigenstates of both L, and L?. We will do this shortly. 
The operator L? is Casimir operator, which means that it commutes with all angular momentum 
operators. Just like it commutes with L,, it commutes also with Lz and Ly. 


To understand the angular momentum operators a little better, let’s write them in spherical coor- 
dinates. For this we need the relation between (r,9,¢) and the cartesian coordinates (2, y, z): 


x =rsin@cosd¢, r =vVax*+y24 22, 
y =rsinésing, 0 = cos”! (2) ‘ (2.15) 
a =reos?e , ro) = tan”! (4) é 
We have hinted at the fact that angular momentum operators generate rotations. In spherical co- 
ordinates rotations about the z axis are the simplest: they change ¢ but leave @ invariant. Both 


rotations about the x and y axes change @ and ¢. We can therefore hope that L, is simple in spherical 
coordinates. Using the definition L, = %p, — yp, we have 


ie *(05 = ve) (2.16) 
Notice that this is related to ¥ since, by the chain rule 
0 
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where we used (2.15) to evaluate the partial derivatives. Using the last two equations we can identify 


Digs Ste (2.18) 


This is a very simple and useful representation. It confirms the interpretation that L, generates 
rotations about the z axis, as it has to do with changes of ¢. Note that L, is like a momentum along 
the “circle” defined by the ¢ coordinate (¢ = @+ 27). The other angular momentum operators are a 
bit more complicated. A longer calculation shows what we suggested earlier, that 


L? LO Seb 1 @& 
he sano nag | + agaae en) 


3 Ejigenstates of Angular Momentum 


We demonstrated before that the Hermitian operators L, and L? commute. We now aim to construct 
the simultaneous eigenfunctions of these operators. They will be functions of 9 and ¢ and we will call 
them wpn(0,¢). The conditions that they be eigenfunctions are 


Letom = hime, meR 


s. 3.1 
Bem = WEE+1) Wm, ER. es 


As befits Hermitian operators, the eigenvalues are real. Both m and / are unit free; there is an A in 
the Lz eigenvalue because angular momentum has units of h. For the eigenvalue of L? we have an h?. 
Note that we have written the eigenvalue of L? as ¢(€+1) and for @ real this is always greater than or 
equal to —1/4. In fact @(¢+ 1) ranges from zero to infinity as @ ranges from zero to infinity. We can 
show that the eigenvalues of L? can’t be negative. For this we first claim that 


(b, Lb) > 0, (3.2) 


and taking 7) to be a normalized eigenfunction with L? eigenvalue \ we immediately see that the above 
gives (w, rv) = \ > 0, as desired. To prove the above equation we simply expand and use Hermiticity 


(b, L?y) = (o, L2y) + (¥, Ly) + (¥, ey) 
= (Lyd, Leb) + (Lyb, Lyb) + (£2), Lp) > 0, 


because each of the three summands is greater than or equal to zero. 


(3.3) 


Let us now solve the first eigenvalue equation in (3.1) using the coordinate representation (2.18) 
for the L, operator: 


h OwWem Ove : 
= = Amen = ae A 
86 mpm 4% > b ime (3.4) 
This determines the ¢ dependence of the solution and we write 
Vem(8,¢) = e'™? P78), (3.5) 


where the function P;"(@) captures the still undetermined @ dependence of the eigenfunction Wem. We 
will require that wen be uniquely defined as a function of the angles and this requires that 


Wem (8, ot 27) = Wem (9, ?) . (3.6) 


There is no similar condition for 6. The above condition requires that 


eim(o+2r) = eime -» e2rim _ 1. (3.7) 


This equation implies that m must be an integer: 


meZ. (3.8) 


‘One may have tried to require that after ¢ increases by 27 the wavefunction changes sign, but this does not lead to 
a consistent set of Wem’s. 


This completes our analysis of the first eigenvalue equation. The second eigenvalue equation in (3.1), 
using our expression (2.19) for L?, gives 


Le EO Yes O 1 
= (S55 (sine | sagan) Yim = WL(L+1)dem- (3.9) 


We multiply through by sin? 6 and cancel the h? to get 


2 
(sin g sin 0 g 2 


— in? 
a6 a9 * aga) vem = —U(0 + 1)sin® Oem . C20) 


Using Wyn, = e”"? P!”(6) we can evaluate the action of cae on Wem, and then cancel the overall e’”? to 
é ad 
arrive at the differential equation 


sind (sing 76 ) m°P;” = —0(€4+ 1)P;”" sin’ 6, (3.11) 
or, equivalently, 

, d ‘ dP;” ; - 2 2 mi __ 

sind (sind 7 ) + (£(€ + 1)sin* 6 — m*) Pi” = 0. (3.12) 


We now want to make it clear that we can view P;” as a function of cos@ by writing the differential 
equation in terms of « = cos@. Indeed, this gives 


Oi nar > sine Sara 


dO dOdx dx dé dx ote) 
The differential equation becomes 
d dP}" 
2 2 £ \ \ 2 2) pm oo 
(1 2”)—|(1 x?) | + [£(€+ 1)(1 — 2?) — m2) P(x) = 0, (3.14) 
and dividing by 1 — x? we get the final form: 
d dP") | m? “2 _ 
=|(a oe ! lee +1) - 3|P G=o: (3.15) 


The P;"(x) are called the associated Legendre functions. They are not polynomials. All we know at 
this point is that m is an integer. We will discover soon that @ is a non-negative integer and that for 
a given value of @ there is a range of possible values of m. 


To find out about @ we consider the above equation for m = 0. In that case we write P)(x) = P?(z) 
and the P,(x) must satisfy 


d [a 22) ah 


= A 490 PIG) =O: (3.16) 


This is the Legendre differential equation. We try finding a series solution by writing 


RiG)= Sap, (3.17) 
k=0 


assuming that P;(x) is regular at x = 0, as it better be. Plugging this into the DE yields we find that 
the vanishing of the coefficient of «* requires: 


(K+ 1)(K + 2)anyo + (C0 +1)-—k(K+1)]a, = 0. (3.18) 
Equivalently, we have 
Ak42 eé+1) —k(k +1) 
= : 3.19 
aK (k + 1)(K+4+ 2) a 
The large & behavior of the coefficients is such that unless the series terminates P; diverges at x = +1 


(since « = cos@ this corresponds to 0 = 0,7). In order for the series to terminate, we must have 
(€+1) =k(k+1) for some integer k > 0. We can simply pick ¢ = k so that ag42 = 0, making P, (2) 
a degree k polynomial. We have thus learned that the possible values of @ are 


£= 0,1,2,3,.... (3.20) 


This is quantization! Just like the m values are quantized, so are the ¢ values. The Legendre polyno- 
mials Po(xz) are given by the Rodriguez formula: 


P(x) = wala,) (2? i) (3.21) 


The Legendre polynomials have a nice generating function 


= 1 
Pigiss == — 2, 3.22 
y (a) V1—2¢s+4+ 8? ( ) 


A few examples are 
Poe) = 1; Pilz) = 2, P(e) = ¢(Ge*—-1). (3.23) 


P;(x) is a degree ¢ polynomial of definite parity. 


Having solved the m = 0 equation we now have to discuss the general equation for P;"(x). The 
differential equation involves m? and not m, so we can take the solutions for m and —m to be the same. 
One can show that taking |m| derivatives of the Legendre polynomials gives a solution for Pj"(x): 


PI"(x) = (1 —2?)Iml/2 (=) P,(a). (3.24) 


Since P) is a polynomial of degree ¢, the above gives a non-zero answer only for |m| < @. We thus have 
solutions for 


See: (3.25) 


It is possible to prove that no other solutions exist. One can think of the Wg, eigenfunctions as first 
determined by the integer @ and, for a fixed @, there are 20+ 1 choices of m: —@,—£+1,...,2. 


Our Wem eigenfunctions, with suitable normalization, are called the spherical harmonics Y;,,,(9, @). 
The properly normalized spherical harmonics for m > 0 are 


Yem(0,¢) = | : 4 es mt 1)"e'™? P™ (cos 8). (3.26) 


For m < 0, we use 
Yom (9, ¢) = (-1)™ [Ye-m(6, O)]" - (3.27) 
We thus have 


Beye = himYim, 


‘ 3.28 
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The first few spherical harmonics are 
1 
Yo,o(9, 6) = = (3.29) 
4dr 
3 4 3 os } 
¥i,41(0,9) = ¥4/ —e** sind = F4/ —-—* (3.30) 
T 87 Or 


| 3 3 2 
Yi 0(8, om) = ae cos 6 = re rv (3.31) 


Being eigenstates of Hermitian operators with different eigenvalues, spherical harmonics with different 
£ and m subscripts are automatically orthogonal. The complicated normalization factor is needed to 
make them have unit normalization. The spherical harmonics form an orthonormal set with respect 
to integration over the solid angle. This integration can be written in many forms: 


[eo = [ae i sind = [a6 [aos ae (3.32) 


The statement that the spherical harmonics form an orthonormal set with respect to this integration 
means that 


/ dQ YZ sy(8,4)Vemn(0,6) = 50,0°Sins (3.33) 


4 ‘The Radial Wave Equation 


Let us now write an ansatz for the solution of the Schrodinger equation. For this we take the product 
of a purely radial function Rge(r) and a spherical harmonic 


wr, 0, ¢) = Ree(r)Yem(9, ?). (4.1) 


We have put subscripts F and ¢@ for the radial function. We did not include m, because, as we will 
see the equation for R does not depend on m. We can now insert this into the Schrédinger equation 
(1.13) 
h2 1 62 f 2 
——-— > (rREeeyYi 
2m r Or? ("ReeYem) + 


Since the spherical harmonics are L? eigenstates we can simplify the equation to give 


Smit RewYom + V(r) ReeYim = EREYm- (4.2) 


h? 1 d?(rRE) Vie h£(0 + T} 
2nr dr? Ane Qmr? 


RrYm +V(r)ReYim = ERgYm.- (4.3) 


Canceling the common spherical harmonic and multiplying by r we get a purely radial equation 


h? d?(rRuv) 
2m dr? 


h2e(€ +1) 
Qmnr? 


(rRe¢) + V(r)(rRee) =H, (rRee), (4.4) 
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It is now convenient to define 
uge(r) = rRe(r). (4.5) 


This allows us to rewrite the entire DE as 


h? dure 
2m dr? 


“ls (V(r) =f aa = Eure. (4.6) 


This is called the radial equation. It looks like the familiar time-independent Schrédinger equation 
in one dimension, but with an effective potential 


A20(€ +1) 


Vorr(r) > V(r) + Oran 


. (4.7) 
that features the original potential supplemented by a centrifugal term, a repulsive potential propor- 
tional to the angular momentum squared. Because of this term, the radial equation is slightly different 
for each value of £. As anticipated, the quantum number m does not appear in the differential equation. 
The same radial solution uge(r) must be used for all allowed values of m. 

Recall our decomposition of the wavefunction: 


U r 
(70,8) = Rpelr ¥en (0,8) = MY, (8,6). (4.8) 
The normalization condition requires 
2 
ee pes b|? = | Pavan une VAY as (4.9) 


The angular integral gives one, the explicit factors of r cancel and we get 
(oe) 
i dr |uge|” = 1. (4.10) 
0 


Indeed uge(r) plays the role of a one-dimensional wavefunction for a particle moving in the effec- 
tive potential along r. Since only r > 0 is allowed, we must consider the possible behavior of the 
wavefunction for r = 0. 


We can learn about the behavior of the radial solution at the origin under the reasonable assumption 
that the centrifugal barrier dominates the potential as r + 0. In this case the most singular terms 
of the radial differential equation must cancel each other out, leaving less singular terms that we can 
ignore in this leading order calculation. So we set: 


h? d? ure h2 ee + if) 
2m dr? Qmr? 


upe = 0, asr—-0. (4.11) 


or equivalently 
ure — &¢+1) 
dr? r2 


The solutions of this can be taken to be uge = r* with s a constant to be determined. We then find 


UE? - (4.12) 


s(s—1) = &@€4+1) > s=4+1, s=-42, (4.13) 
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thus leading to two possible behaviors near r = 0: 


1 

une ~ rt, or une ye (4.14) 
r 

For € > 0 the second behavior is not consistent with normalization, the wavefunction diverges at r > 0 

far too quickly. For ¢ = 0 the second behavior, leading to R ~ 1/r, is in fact not a solution of the 

Schrodinger equation. Therefore we have established that for all 2 > 0 we must have 


upe ~wer't!, asr— 0. (4.15) 


Note that wze vanishes at r = 0. Even for @ = 0, we have u ~ r and wu vanishes at r = 0. Effectively 
there is an infinite wall at r = 0 consistent with the impossibility of extending r to negative values. 


Recall that the full radial dependence of the wavefunction is obtained by dividing uge by r, so that 
Ree w~ cr’. (4.16) 


This allows for a constant non-zero wavefunction at the origin only for = 0. Only for ¢ = 0a particle 
can be at the origin. For ¢ 4 0 the angular momentum “barrier” prevents the particle from reaching 
the origin. 


Sarah Geller and Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX 
version of this document. 
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1 The Hydrogen Atom 


Our goal here is to show that the two-body quantum mechanical problem of the hydrogen atom can 
be recast as one in which we have center-of-mass degrees of freedom that behave like a free particle 
and relative-motion degrees of freedom for which we have dynamics controlled by a central potential. 

The hydrogen atom consists of a proton and an electron moving in three dimensions. We label the 
position and momentum operators of the proton as Xp, Ppp, and those of the electron as Xe, pe. These 
are canonical variables, meaning they satisfy the canonical commutation relations: 


[(Xp)i, (Pp) j] = *hdiz,  [(Xe)i, (De) sj] = thdi;. (1.1) 


Here the subscripts 7,7 = 1, 2,3 denote the various components of the vector operators. Furthermore, 
the proton variables commute with the electron variables. We have two pairs of independent canonical 
variables. 

The wavefunction for the system is a function of the positions of both particles: 


U (Xp, Xe) , (1.2) 


and the quantity 
| Wx, xe) |? dag xe, (1.3) 


is the probability to find the proton within a window d°®x, of x, and the electron within a window 
d?x, of x. The Hamiltonian of the system is given by 


a) a) 
Pp Pe 

= ——+4+2-+-+4+V — . 1.4 
2Mp 2Me ([e xp) ( ) 


Note that the kinetic energy is simply the sum of the kinetic energy of the proton and kinetic energy 
of the electron. The potential only depends on the magnitude of the separation between the two 
particles, not on their individual positions. 

In order to simplify the problem, we will introduce two new pairs of independent canonical variables. 
The first pair is associated with the center-of-mass (CM) motion. We introduce the total momentum 
operator P and the CM position operator X, given by 


MeXe + MpXp 


(1.5) 


Met ing 


The operator X is given by the typical expression for the center-of-mass of the system, but with the 
positions replaced by position operators. Using the commutation relations (1.1), we can show that X 
and P are canonical conjugates: 


——1 
“~~ 
»~ 
ad 
— 
—~ 
> 
— 
& 
Lal, 
| 


be Teal 6 2 6 


Me + Mp 
= fie 1 tit [ey ) (De) j] + viene [(Xp)i ? Pp) | (1.6) 
= a thd; ; me thoi; 5 
Me = Mp Me + Mp 
resulting in the expected 
(&. (P);| = thd;; . (1.7) 


For the second pair of canonical variables we will define relative position and momentum operators. 
The relative position operator is the natural variable implied by the form of the potential: 


X = Xe—Xp. (1.8) 


Since the second pair of canonical variables must commute with the first pair, we must check that 
x defined above, commutes with X and with P. The commutation with X is automatic and the 
commutation with P works thanks to the minus sign in the above definition. We must now construct a 
relative momentum operator p that is canonically conjugate to x. It must be built from the momentum 
operators of the two particles, so we write 


P= QPe— SPp, (1.9) 


with a@ and £ coefficients to be determined. To be canonically conjugate, the relative operators must 
satisfy 
[(%)i,(B)j] = 155 > atB =1, (1.10) 


using the above definitions of x and p and the proton and electron commutators. Finally, the relative 
momentum must commute with the CM coordinate 


(Xs, (p)| =0 + mea- mB =0. (1.11) 


The two equations for a and § can be solved to find 


Mp Me 


Me + Mp B Me + Mp ( ) 
We define the total mass M and the reduced mass yz as follows 
MeMp 
M= Se 1.13 
Me + Mp, fh Me + Mp ( ) 


The reduced mass of a pair of particles with very different masses is approximately equal to the mass 
of the lower-mass particle. Using these definitions 


a=, pH. (1.14) 


Thus, collecting the relative variables we have 


p = u(2 Pe) = —? Be =e Bp R= kK. (1.15) 
Note that the relative momentum p can be written in terms of velocities as follows: p = [u(ve — vp). 
The relative momentum vanishes if the motion is only CM motion, in which case the velocities of the 
two parties are the same. 

We can now rewrite the Hamiltonian in terms of the new variables. Solving for the original 
momentum operators in terms of P and p, we find 


an Mp * a Me x 
BS a Pie Pe S 5, Pte: (1.16) 


We can then rewrite the kinetic terms of the Hamiltonian in the form 


~ 2 oo) 2 
1 Ms» 2 R 
Pp de Pe = FP? Pp 6 p? 
2Mpy =2Me 2Mp \ M M 
Vi fied Dilip | 25 
P? 4 P-p4 1.17 
Ome (Fs Ne tee iG 
p2 p? 
~ 2M" Qu 


Happily the term coupling the two momenta vanishes. Thus the center of mass degrees of freedom and 
the relative degrees of freedom give independent contributions to the kinetic energy. The Hamiltonian 
can then be written as 


ot oh (1.18) 


In position space, the total and relative momentum operators can be expressed as gradients 


‘ h h 
i i 
Each V has a subscript indicating the type of coordinate we use to take the derivatives. Just like we had 
a wavefunction U(x,-,x,) the new canonical variables require that we now think of the wavefunction 
as a function U(X, x) of the new coordinates. 

We solve the time-independent Schrédinger equation by using separation of variables 


W(X, x) = Vom (X)Vye1(x). (1.20) 
Plugging these into the time-independent Schrédinger equation HU = EW, we reach 


A 2 


p 
Wre Ta 
oa Fe 


P2 
ea (CX 
ag You) 


Drei(x) + V (la) oa) Wo (X) = EVoy(X)Vye(x). (1.21) 


Dividing by the total wavefunction Voy (X)WVyre1(x), this becomes 


P2 
er p< 
aug Yom (%) 


2 
Wom (X) 


Qu 


| E | a) Ui ae: (1.22) 


The first term on the left-hand side is a function of X only and the second term on the left-hand side 
is a function of x only. Their sum is equal to the constant FE and since x and X are independent 
variables, each term must individually be constant. We thus set the first term equal to the constant 
Eom and the second term equal to the constant F,.|, resulting in the following equations: 


Feu X) <3) Baek). (1.23) 
= 2 
aaa) Dre(x) = EverVrei(X), (1.24) 
Lb 


E= Eom + Fre. (1.25) 


We get two Schrodinger equations. The first equation tells us that the center of mass moves as a free 
particle of mass M. Thus, the CM energy is not quantized and we get plane wave solutions. The 
second equation is for the relative motion, and as we wanted to show, it is described as motion in a 
central potential. The third equation tells us that the total energy is the sum of the center-of-mass 
energy and the energy from the relative motion. 


2 Hydrogen atom spectrum 


We now have the tools to study the hydrogen atom, which has a central potential given by 


Ze? 
Vee, (2.1) 
r 
where Z is the number of protons in the nucleus. For hydrogen we have Z = 1. But it is worth 
considering Z > 1 in which case we are describing the motion of an electron around the nucleus of 


some other atom. We will also define following physical constants. 


2 
e The fine structure constanta: a= < ~ a 


e The Bohr radius ag. This is the characteristic length scale in the problem. It can be calculated 
by equating kinetic and potential energies expressed in terms of ap and ignoring all numerical 
constants: 

h2 e2 


eS (2.2) 
Med ao 


Here the mass should be the reduced mass, which in this case can be taken rather accurately to 
be the mass of the electron. We then have explicitly, 


_ ie = ee = he ihe 
10 me ema (<)me? ame? 
197 MeV fm 1970 eV Angstrom (2.3) 
= = = - x 137 
0.51 x 10 eV(z347) 0.51 x 10% eV 
= 0.529 Angstroms ~ 53pm. 
For the energy scale estimate we have 
2 2 4 
eof me*\ _ e€ a ee oh 


There are other characteristic lengths that are interesting: 


aaj = Compton wavelength of electron = A. ~ 390fm, 


classicalelectronradius ~ 2.8fm. 


2) 
g 
S 

II 


Let us do some work now! The radial Schrodinger equation for bound states E < 0 reads 


u= Eu. 


( 2 d? | Ree +1) 2 


Imdr2 '  2mr2 r 


(2.6) 


We could label the wavefunction u as ug as the solutions will certainly depend on ¢ and the energy E. 


As usual, we like to work with a unit free coordinate. This could be achieved by writing r = agx, with 


x unit free and ag carrying the length units of r. It will be more convenient to use a slight variation 


to eliminate Z from the equation and some factors of two. We will take the new unit-free coordinate 


xz to be defined from 


The Schrédinger equation then becomes 


( Rh? 4Z? d? AZ WI +1) 22%? -) 
u= Eu 


2m. .ap dx? an 2m. 2? ag 2x 
oh? Z? @ i(i+1 277671 
> ( 5 ( a ( BE ’) : u = Eu. 
mag dx x ag « 


Note that 
Dn? 7? - 2h? Z me? Ze? 


maz mag hi? fig.” 


which reduces our differential equation to 


a  @€+1) 1\  #23&£ 
dz2 x? ry (222°) we 
ao 


We now define the unit-free parameter « that encodes the energy: 


Ko = eerie > 0. 
(ae) 


& is a unit-free version of the bound state energy. The differential equation is then 


(2.7) 


(2.10) 


(2.11) 


(2.12) 


We can further simplify this equation prior to solving it by examining the limiting cases. In the limit 


x — oo, the dominant terms are the second derivative and the term on the right-hand side, giving 


du 9 + 
a HKU ve 
dx 


(2.13) 


Since « is unit free, we can make the above exponent equal to a new unit-free coordinate p: 


p= Ke = —r. (2.14) 


This time we get 


da ee+1) 1 
= -u. 2.1 
( dp? p za) = oa 


Note that we did not get « to disappear from the equation. This is good news: the equation should 
fix the possible values of « (or possible energies). The equation above is not quite ready for a series 
solution: we would find a three term recursion relation, which is rather complicated. To make progress 
we discuss the behavior for small and large p. 


For p > co we now get u ~ u*? and of course we hope for u = e~? for normalizability. As we 
discussed before, for p —> 0 the radial solution must be of the form u ~ p+), This information about 
the behavior for small and for large p suggests a good ansatz for u(p) 


u(p) = pt W(p)e?. (2.16) 


where W(p) is a yet to be determined function that we hope satisfies a simpler differential equation. 
To derive this differential equation for W(p), we plug our ansatz into Eq. (2.15). As a little help on 
the calculation, we give an intermediate result: 


e( +1) 
2 


2(€ +1) 
p 


gh Ww" 4 2(€+ 1) 


utu = ( w" W 4 2W") pftlen?. (2.17) 


With a little more work we finally get the differential equation for W: 


PW dw yl 
page +24) 


Oa 2(L 4 )| W = 0. (2.18) 


This looks a bit more complicated than the differential equation we started with but it leads to a very 
nice one-step recursion relation. As usual we write W as a series expansion 


W= 5S axp*, (2.19) 
k=0 


and plugging back into (2.18), we group terms of order p* to derive a recursion relation 


1 
ansik(k +1) +2(€4+1)(k + Langs — Qkox 4 E 204 1) a, err 


(2.20) 
1 
> aggilk(e +1) +2(2+1)(6+)) = a (2(¢+ k+1)- -) 
which gives 
ak, (K+ 1)(k + 2€ + 2) : 


Detailed examination shows that for normalizable wave functions, the series must terminate. To see 
this we examine the large & behavior of the above ratio: 


QAk+1 2 2k 2 


~ = ; 2.22 
ak k?2 k ( ) 
Note that mi < 2; thus if the ratio ma leads to a divergence so will the ratio 2. Taking 
Gk+1 2 2 
— + = —  @Q 7 2.23 
te ka +1 = 7 41% (2.23) 
and this is solved by 
9k 
Therefore, the sum 
W= Ne axp* ~ ao Ss" Sr age”? . (2.25) 
k=0 k=0 


This is precisely sufficient to make the ansatz in (2.16) un-normalizable. 
In order to get a normalizable solution the series for W must terminate. Suppose W is a polynomial 
of degree N so the coefficients satisfy 


an #O and anii=O0. (2.26) 


From Eq.(2.21) this implies 
1 
— = 2N+¢+1). (27) 
kK 
Quantization has happened! The energy-encoding parameter « is now related to integers! Note that @ 
can take values @ = 0,1,2,... as it befits it being an angular momentum quantum number. Moreover 
N can take values N = 0,1,2,..., since a polynomial of degree zero exists, being equal to a constant. 
Define the principal quantum number n as follows: 


1 
n= N+é+1= —, with £=0,1,2,..., N=0,1,2,..., andn=1,2,3,....) (2.28) 
K 


Importantly, note that for a fixed n we must have 
0 <2 <n-1, and 0<N <n-l. (2.29) 


If n and @ are known, N is determined from N + @+1=n. So the independent quantum numbers so 
far are n and &. Interestingly, the energies depend only on n, since « depends only on n. Using Eq. 
(2.11), the energy dependence on principle quantum number is given by 


972 2 
pea. (2.30) 
ao 
and using & = + we get 
wre 
E = ——->5 2.31 
2a9 n? 28)) 


These are the energy levels of the hydrogen atom! Since at any fixed value of n > 1 there are various 
possible @ values, the spectrum is highly degenerate. Even more, each value of @ amounts to 20+ 1 
states, given the possible values of m. One way to visualize the spectrum is shown in Figure 1. All 
integer points in the (N, 2) positive quadrant represent states. The states with common value of n lie 
on the dashed lines. 


Figure 1: All points with integer N > 0 and integer @ > 0 represent hydrogen atom states. The figure helps 
us to count the number of possible states for given value of n. Each dot along the diagonal line for a given n 
represents a possible state. 


Figure 1 helps us to count the number of bound states for a given value of n. Recall that for each 
n, € can take values from 0,...,2.—1 and for each value of @, m takes values from —£ to @. The following 
table counts the states for the first few values of the principal quantum number n. A given state is 
specified by its values for (n,1,m), all of which are referred to as thequantum numbers for hydrogen 
states. Each number has a very important physical meaning: n tells us about the energy eigenvalue, 
hl(¢ + 1) is the eigenvalue of the square of angular momentum and fim is the eigenvalue of the z 
component of angular momentum. 


Number of States 


n value | £ values m, values total states 
n=1 £=0 m=0 1 state 
n=2 £=0 m=0 1 
l=1 m=-—1,0,1 + 3 

= 4 states 

n=3 l=0 m= 1 
l=1 m=-—1,0,1 +3 

(=2 | m= —2,...,2 +5 

= 9 states 


The total number of states for arbitrary principal quantum number n can now be calculated: 
n-1 

2in-—1 
# of states forn = S°(2¢+1) = letra er eee eee (2.32) 


2 
£=0 


This is in agreement with the partial results in the table. A more familiar representation of the states 
of hydrogen is given in Figure 2. The different columns indicate the different values of @. We have 
also indicated in the figure the values of N, the degree of the polynomial entering the radial solution. 
Note that for a given @, that is, for a fixed radial equation, the value of N increases as we go up the 
column. The number N corresponds to the number of nodes in the solution. 
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Figure 2: Plot of energy levels E ~ —1/n? indicating also the angular quantum number @ and the degree N of 
the polynomial. The spectrum is highly degenerate. 


Recall that we defined p = ae Together with « = a this gives 


Z 
pee; (2.33) 
nag 


The eigenstates are labelled by quantum numbers (n, ¢,m) and the wavefunctions are 


+1 
nim = NM) (0,8) = NM ——Waelp}en? Yin( 6,8) 


= Np Wrelp) €? Yim(9, 4), (2.34) 
ee 
polynomial of degree 
N=n—(é+1) 


where N is a normalization constant. Therefore, using the expression for p and absorbing constants 
into NV we have 


] i r Zr 
a ee (=) ( polynomial in = ) Sor Ue (2.35) 


ao of degree N=n—(€+1) 


For the ground state of hydrogen (Z = 1), we have (n, @,m) = (1,0,0). Having zero angular momentum 
the associated wavefunction has no angular dependence. The normalized wavefunction is 
1 —r/a 
broo(r,8,¢) = qe! (2.36) 
For normalized hydrogen wavefunctions at n = 2 and n = 3 see 
http: //hyperphysics.phy-astr.gsu.edu/hbase/quantum/hydwf.html 


Sarah Geller and Andrew Turner transcribed Zwiebach’s handwritten notes to create the first LaTeX 
version of this document. 
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